Real Time Dynamics of Symmetry Breaking 



Zsolt Szep 
Ph.D. Thesis 

Supervisor: Andras Patkos 



physics program 
program leader: Horvath Zalan 
subprogram leader: George Pocsik 



Department of Atomic Physics 
Eotvos Lorand University, Budapest 



September 2001 



Contents 



1 Introduction and motivation 



2 General formalism for non-equilibrium QFU] 

77L Schwinger- Dyson approach 

Kinetic equations for the two-point functions 

^.1.2 Iterative solution of the fieisenberg equations 

\>,.1.'S Mode-function expansion 

2.2 Classical linear response theory 

3 One-component <l>"^ theory 

3.1 Classical kinetic theory of Landau Damping] 

p. 1.1 Kinetic theory of particles with background held dependent mass . . . 

p. 1.2 The wave equation of the slow modes 

p. 1.3 Connection with the formalism of Danielewicz and Mrowczyhski . . . 

p. 1.4 The induced source 

p. 1.5 Discussion] 

J72 Ettective order parameter dynamics and the decay of a metastable vacuum state 

far from equilibrium 

p. 2.1 Classical cut-ott" <1> "^-theory on lattice 

p. 2. 2 Time-history of the order parameter 

F2.3 Motion near the (meta)stable point] 

p. 2. 4 Jump from the metastable to the stable vacuum 

p. 2. 5 Connection with the nucleation picture 

p. 2. 6 Conclusions 

4 The damping of the Goldstone modes 

1.1 Etfective theory for the soft modes in 0(JN) model 

1.1.1 I'he ettective equations of motion for the slow modes 

1.1.2 'I'he two-point function of the fast modes 



i 



CONTENTS 



4.1.3 Non-equilibrium linear dynamics of the slow modes | 69 

4.1.4 Dispersion relations for on-shell waves] 70 

4.1.5 Solution of the initial value problem ot" the slow tields| 75 



4.1.6 Non-local dynamics of the slow modes | 77 



4.1.7 Conclusions 82 



4.2 Nonlinear relaxation of classical 0(2) symmetric fields in 2 -|- 1 dimensions 



4.2.1 The model and its numerical solution 



83 
83 



4.2.2 Relaxation to the h ^ equilibrium] 84 

4.2.3 The h ^ equilibrium] 89 

Diffusion to the h = equilibrium] 91 



4.2.4 



4.2.5 Conclusions] 91 



5 Real time dynamics of the Higgs effect 



93 



5.1 Mean Field Equations in the Abelian Higgs model] 94 



5.2 High-temperature expansion beyond the HTL approximation 



5.2.1 Even terms 



96 
97 



5.2.2 Odd contributions] 98 

99 

100 



5.2.3 Linear response and physical modes 

5.2.4 The infrared separation sc"ale 

573 Discussion] 101 



A The introduction of the one-particle distribution function and the calculation 



o ^ the damping rate 



105 



B Equivalence of the one-loop perturbation theory and the once iterated solu- 
ti jjn of the generalised Boltzmann-equations ] 109 



U The dynamics of the classical 0{J\) model 

D Evaluation of some integrals appearing in Subsection 4.1.3] 



111 
117 



u 



Chapter 1 

Introduction and motivation 



In the last few years remarkable progress has been made in the investigation of the non- 
equilibrium behaviour of finite temperature systems, both for scalar and gauge fields. Actually, 
after understanding equilibrium features of these systems many of the workers in the field have 
eagerly blazed away at investigating non-equilibrium aspects of finite temperature Quantum 
Field Theories, with a special emphasis on the real time simulation of its high temperature 
limit. There are a few good reasons for doing this. 

First of all, time dependent phenomena in hot relativistic quantum field theory play an 
important role in applications to cosmology, in the context of baryogenesis and inflation, and 
in the heavy ion collisions. Not all of these phenomena are non-perturbative (for example 
relaxation phenomena to equilibrium from states slightly out of equilibrium can be treated 
perturbatively) but as the departure from equilibrium increases non-perturbative tools are 
needed. Non-perturbative calculations under general circumstances are very difficult. Real time 
numerical simulations in the QFT would face the problem of complex weights appearing in the 
quantum expectation value of time dependent quantities. Fortunately a classical description 
can be accurate for physical quantities which are determined by low momentum modes of 
the theory, because these modes are highly occupied at high enough temperature. For high 
momentum modes, the classical treatment is incorrect, but in general the influence of these 
modes can be taken into account perturbatively. Initially proposed in for the computation 
of the rate of sphaleron transitions, the numerical solution of the classical approximation proved 
to be eventually a powerful tool in evaluating quantities for systems developing in time. 

Preheating and the possibility of Electroweak baryogenesis 

An intensively investigated field is the transition from an inflationary to the radiation dom- 
inated epoch in the early Universe. In the context of chaotic inflation it was realised that 
the initially homogeneous inflaton field can decay very rapidly into narrow bands of its own 
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low momentum modes or into modes of other scalar fields through parametric resonance. This 
phenomenon goes by the name of preheating. 

The large amplitudes of the parametrically amplified modes make this problem non-pertur- 
bative, but at the same time the large occupation number of these modes opens the possibility 
to study this problem in the classical field approximation on lattice (for example and refer- 
ences therein). It was then realised that in addition to the standard high-temperature phase 
transitions that is assumed to occur in the state of thermal equilibrium |Q, there exists a new 
class of phase transitions which may occur at the intermediate stage between the end of infla- 
tion and the establishing of thermal equilibrium These new types of phase transitions 
are referred as to cosmological non-thermal phase transitions after preheating, and can restore 
symmetry on scales up to 10^^ GeV, due to large amplitude non-thermal fluctuations. An 
important feature of the non-thermal phase transition based on preheating is that the classi- 
cal non-linear dynamics of infiaton field coupled with other boson fields can induce symmetry 
breaking for the latter. 

Symmetry breaking is a fundamental component of theoretical particle physics because 
through the Higgs effect it accounts for the masses of the elementary particles. These particles 
have got their masses during the electroweak phase transition. The importance of this phase 
transition is revealed by its impact on the observed baryon asymmetry in our Universe. A 
successful theory of baryogenesis requires beyond the existence of baryon number violating 
processes C and CP violating processes and departure from thermal equilibrium. These three 
requirements are known as Sakharov conditions. In the standard model (SM) baryon number 
is not conserved because of the non-perturbative processes in the gauge sector that involve the 
quantum anomaly. The usual scenario assumes that the universe was in thermal equilibrium 
before and after the electroweak phase transition, and far from it during the phase transition. 
In order to significantly drive the primordial plasma out of equilibrium a strongly first order 
phase transition is needed. This is also a condition for the sphaleron processes to stop quickly 
after the transition in order to prevent the washing out of the created baryon asymmetry. The 
order of the EW phase transition was thoroughly analysed, both perturbatively [0 and with 
MC simulations p|, and the conclusion was, that there is an end point on the line of the first 
order phase transition separating the symmetric and broken phases at a Higgs mass of about ^ 
72 GeV. For greater Higgs mass values there is no first order PT that makes the things hopeless 
in view of the LEP experiments that exclude a Higgs mass below 113 GeV. So, the possibility 
of baryogenesis in the usual scenario is ruled out in the SM. This also means that if the EWPT 
occurred in thermal equilibrium then the onset of the Higgs effect was a smooth process. 

In the usual scenario one can go to a minimal supersymmetric extension of the SM. The 
phase structure of MSSM is much more difficult to explore because of its large parameter 
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space. According to perturbation theory and 3- dimensional lattice simulations |T0[ as well as 



4-dimensional lattice simulation |]Tl|, there could be an EWPT in MSSM that is strong enough 
for baryogenesis up to a value of ~ 105 GeV for the lightest Higgs mass. 

Recently it was shown that the baryogenesis can take place below the EW scale in a very 
economical extension of the SM in which only one more scalar SU{2) x [/(l)-singlet field is 
required, to be taken the inflaton |13|. The novel scenario make use of the possibilities that 
reside in the phenomena of preheating. The preheating precedes the establishment of thermal 
equilibrium. With an appropriate choice of the parameters of the inflation one can achieve 
that the reheating temperature remains below the EW scale, therefore no intensive sphaleron 
processes are allowed. There is an important baryon asymmetry production during rescattering 
after preheating due to the fact that the low modes that carry a large amount of the energy 
quickly "thermalize" at a high effective temperature (a few times the value of the EW scale) 
where the rate of sphaleron transitions is very effective. 

Thermalization and relaxation 

The approach to equilibrium of initially out-of-equilibrium states is a highly important issue 
in many branches of physics ranging from inflationary cosmology (the spectrum of density fluc- 
tuation in the early universe) through particle physics (the problem of baryogenesis, formation 
of DCC in heavy ion collisions) to statistical physics (dynamics of phase transitions, realiza- 
tion of Boltzmann's conjecture that is an ensemble of isolated interacting systems eventually 
approaches thermal equilibrium at large times.) 

In the recent field theoretical studies of thermalization and relaxation classical fields are 



intensively investigated The interest stems from the fact that in order to explore the 

reliability of truncation and expansion schemes they should be benchmarked against this exact 
solution. Recently, evolution of equal-time IPI correlation functions derived from the effective 
time dependent action [|15] were confronted with the results of the exact time evolution. By 
solving equations non-local in time, obtained from 2PI effective action thermalization of quan- 



tum fields was demonstrated ||T6|. These approximate methods can be formulated both for 



classical and quantum cases. In the quantum case the analogue of the classical ensemble aver- 
aging over the initial conditions is the quantum expectation value. With this correspondence 



the formal derivations, and hence the results, are quite similar |T7[] 



Heavy Ion collision 

Two phenomena of the QCD chiral phase transition have received much attention in the 
past decade due to the fact that both may be within reach of discovery by the current or 
forthcoming relativistic heavy-ion collision experiments: the formation of disoriented chiral 
condensate (DCC) and the end point (the critical point E) of the first order phase transition 
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line for decreasing chemical potential /i in the jj, — T QCD phase diagram. 

When both areas are studied in the context of heavy-ion collision experiments a knowledge 
of the evolution of far from equilibrium dynamical system of strongly interacting quantum fields 
of bosonic and fermionic degrees of freedom is needed. During the evolution the system rapidly 
expands and cools creating and emitting baryons and mesons. 

As investigating such systems and in particular phase transitions directly in QCD is out of 
question at the moment, a model that encodes relevant aspects of QCD in a faithful manner is 
needed. 

It was argued in [|T^ that the 0(4)-symmetric linear cr-model is representative of the physics 
relevant for the formation of DCC. This model is in the same static universality class as the 
QCD with two massless quarks, that is a good approximation to the world at temperatures and 
energies below Xqcd- The basic idea of DCC formation is that regions of misaligned vacuum 
might occur. These are regions where the field = (cr, vr) instead of taking the true ground 
state value (f , 0) is partially aligned with the vr direction. It is expected that the relaxation 
of the misaligned vacuum region to the true vacuum would proceed through coherent pion 
emission with a charge distribution violating the isospin symmetry. For energetical reasons 
(see Ref. [0 ) the emission of a large number of pions that produce an easily detectable signal 
could only be envisaged in a highly non-thermal chiral phase transition, that is when the system 
is far out of equilibrium. This situation is approximated with a sudden quench from high to 
low temperature during which the long wavelength modes of the pion fields become unstable 
and grow relative to the short wavelength ones. 

At the critical point E the phase transition is second order and belongs to the Ising univer- 
sality class. The pions remain massive but the correlation length of the a field, that plays the 
role of the order parameter, diverges due to growing long wavelength fluctuations. An increase 
in the correlation length compared with the one corresponding to the equilibrium energy den- 
sity of the QCD plasma can be observed by analysing the low transversal momentum pions in 
which the a fields decay. 

At present it is not possible to construct the exact mapping between the axes of QCD phase 
diagrams and the axes of the Ising phase diagrams. With the new possibility of simulating 
the QCD phase diagram at non-zero chemical potential proposed in this question may be 
solved in the near future. 

Critical slowing down that occurs near the critical point E requires long equilibration times 
meaning that the plasma will inevitably slow out of equilibrium if the plasma is cooled. Since 
the rate of cooling is faster than the rate at which the system can adapt to this change, a 
non-equilibrium evolution is guaranteed. 

The phenomenon we are interested in both cases cited above is related to the growth of 
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the correlation length basically produced by the long wavelength, classical modes. This is 
an argument favouring classical real time numerical simulations. The inclusion of quantum 
correction is an important question and has been implemented in mean-field approximation. 

Recently the effect of dynamical behaviour of the a correlation length around the critical 
point E was studied in Ref. |]20|. An intuitive mapping between QCD and the Ising model 



was proposed identifying the magnetic field and the reduced temperature of the Ising model 
with the temperature of the QCD and the chemical potential, respectively. The result of 
the investigation, a factor of 2-3 increase in the correlation length proved to be insensitive to 
moderate tilt in this mapping. In Ref. |2l| basically the same investigation was performed with 



a second order dynamics in the context of $ theory. 

The classical approximation and its problems. 

It is well known that at high temperature the thermodynamics of a QFT in (3 + 1) dimen- 
sions is successfully described through dimensional reduction (DR) by an effective classical 3 
dimensional theory derived for the static modes of the original theory. Its couplings are deter- 
mined by the integration over the non-static modes [^, |2^. Corrections to the dimensional 



reduction are small when the thermal mass and external momenta are small compared to the 
temperature, a requirement that is satisfied when the theory is weakly coupled. Unlike the 
case of QFT at finite temperature where the Bose-Einstein distribution effectively introduces 
an UV-cutoff, the formulation of a classical field theory is meaningful only if a A UV-cutoff 
is present, otherwise one encounters the Ray leigh- Jeans divergences. The introduction of A 
means that the parameters of the effective theory must depend on it in such a way as to cancel 
the A-dependence of the regularized loop integrals in the effective theory. 

Encouraged by the success of DR one could naively expect that the dynamics of soft modes 
is described by classical Hamiltonian dynamics i.e. calculation of time dependent correla- 
tion functions implies solving classical Hamiltonian equations of motion for arbitrary initial 
conditions, over which a Boltzmann weighted averaging with the 3 dimensional DR action is 
performed. This was the essence of the proposal of Grigoriev and Rubakov for calculating 
time dependent correlation functions. The expectation above works only for scalar $^ theory 
and for Abelian-gauge theory in (1 + 1) dimension. In the former case it was shown that the 
ultraviolet divergences of the time dependent correlation functions can be absorbed by the same 
mass counterterms as needed in the static theory ||24|| . 

For non-Abelian gauge theory in higher spatial dimensions the dynamics of infrared fields is 
classical, but they are not described by classical Hamiltonian dynamics. Here the problem of UV 
divergences is more involved. Actually, the existence of a classical theory poses two questions: 
1. what is the divergence structure of the hot real time classical Yang- Mills theory 2. if there 
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are well-separated scales in the theory (for example in a gauge theory) what is the interplay 
between them. The answer to the first question is, that at one loop the linear divergences of the 
classical theory are related to the quantum hard thermal loops (HTL) discovered by Braaten 
and Pisarski. For example, the divergent part of the classical self-energy can be obtained as 
the classical limit of HTL self energy ||2^, |2^, It was realised in that there are UV 
divergences in the classical thermal Yang-Mills theory which are non-local in space and time (on 
the lattice they are sensitive to the lattice geometry), they do not occur in equal time correlation 
functions and so they cannot be absorbed by local counterterms. The second question led in 
the end to the serious challenge of constructing an effective theory for the low-momentum 
modes by incorporating as accurately as possible the effect of the quantum but perturbative 
high-momentum modes. This step has to be done if the classical theory exhibits strong cut-off 
dependence, or equivalently lattice-spacing dependence, because if this occurs it shows that the 
classical dynamics is sensitive to quantum modes with momenta of the order of the cutoff. 

In non-Abelian gauge theory at temperatures where the weak coupling limit is valid, there 
is a hierarchy of three momentum scales that played an important role already in the static case 
p8| , p9[] . There is the "hard" scale T of the typical momentum of a particle, there is the "soft" 
scale gT associated with colour-electric screening and finally there is the "ultrasoft" scale g^T 
associated with colour-magnetic screening and where non-perturbative effects appear. 

The first numerical simulations were done in a purely Yang-Mills classical theory with a 
cutoff A of order T, but it raised the suspicion of a strong cutoff dependence of the measured 
quantity, the sphaleron rate. This indeed was the shown for example in ||30|| . 

One can incorporate the effect of the "hard" modes of order ~ T on the dynamics of softer 
modes using a classical theory based on the HTL effective action obtained after integrating out 
these modes. There are difficulties with this approach due to the fact that the HTL equations 
of motions are non-local in space and time and a local Hamiltonian form exists only in the 
continuum. 

The classical effective HTL theory leads to UV divergences and has to be defined with 
an UV-cutoff gT <^ A <^ T. In order to ensure the UV insensitivity, the coefficients of the 
effective theory (effective HTL Hamiltonian) must be calculated with an infrared cutoff A, so 
that the cutoff dependence of the parameters in the effective theory cancels against the cutoff 
dependence of the classical thermal loops. The problem is, that on the lattice it is very difficult 
to derive the expression of the A (IR cutoff )-dependent piece in the HTL Hamiltonian. 

In order to reduce the sensitivity to the UV scale A one can go further and integrating 
out the soft degrees of freedom down to scale g^T <^ A ^ gT starting from the classical 
effective theory for the soft fields. This was first performed by Bodeker, who showed that the 
resulting theory at the scale g'^T takes the form of a Boltzmann-Langevin equation |3^. It was 
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demonstrated later that this equation is insensitive to the ultraviolet fluctuations 



|3l 



Presentation of the thesis 

This thesis reflects the route followed by its author towards the derivation of the dynam- 
ical equations allowing the study of the real time onset of the Higgs-effect. Since this study 
would involve both analytical and numerical investigations we have developed and tested both 
techniques in simpler models. 

First, I have investigated the dynamics of the one-component scalar theory both in the 
linear-response approximation and then for large deviations from equilibrium. 

Next, I proceeded to the real time characterisation of the Goldstone effect when a condensate 
breaks the 0(N) symmetry of the fleld theory. Here again both quantum and classical systems 
were studied. 

Finally, I have studied the corrections to the Hard Thermal Loop dynamics in the Abelian 
Higgs model, which reflect the presence of the scalar condensate. 
I have flve publications related to my thesis: 

• A. Patkos, Zs. Szep, Phys. Lett. B446 (1999) 272-277 

• A. Jakovac, A. Patkos, P. Petreczky, Zs. Szep, Phys. Rev. D61 (2000) 025006 

• Sz. Borsanyi, A. Patkos, Zs. Szep, Phys. Lett. B469 (1999) 188-192 

• Sz. Borsanyi, A. Patkos, A. Polonyi, Zs. Szep, Phys. Rev. D62 (2000) 085013 

• Sz. Borsanyi, Zs. Szep, Phys. Lett. B508 (2001) 109-116 
The structure of the thesis is as follows. 

In Chapter 2 I briefly review the formalism we use in the following Chapters: the Green- 
function approach to transport theory of quantum scalar flelds developed by Danielewicz and 
Mrowczyhski and the classical linear response theory of Jakovac and Buchmiiller. I show the 
connection between the former method with others existing in the literature. 

In Chapter 3 I present the work done within the one-component scalar $^ model. An 
effective theory of low frequency fluctuations of selflnteracting scalar flelds is constructed in 
the broken symmetry phase coupled to the particles of a relativistic scalar gas via their locally 
variable mass. The non-local dynamics of Landau damping is investigated in a kinetic gas 
model approach. 

Next, I present a numerical investigation of the thermalization in the broken phase of 
classical $^ theory in 2 -|- 1 dimensions. The $ fleld is coupled with a homogeneous external 
"magnetic" fleld that induces a transition from a metastable state to the stable ground state. 
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The dynamics of the system is described using an effective equation of the order parameter. 
This description is consistent with the miclcation theory in a first order phase transition. 

Chapter 4 also consists of two Sections. In the first one an effective theory of the soft 
modes in the broken phase of 0{N) symmetric model is presented to linear approximation 
in the background when the effect of high-frequency fluctuations is taken into account at one- 
loop level. The damping of Higgs and Goldstone modes is studied together with large time 
asymptotic decay of an arbitrary initial configuration. 

In the second Section the real time thermalization and relaxation phenomena are numerically 
studied in a 2 -|- 1 dimensional classical 0(2) symmetric scalar theory. The near-equilibrium 
decay rate of on-shell waves and the power law governing the large time asymptotics of the 
off-shell relaxation is checked to agree with the analytic results based on linear response theory. 
The realisation of the Mermin- Wagner theorem is also studied in the final equilibrium ensemble. 

Chapter 5 deals with the real time dynamics of the Higgs effect. The effective equations of 
motion for low-frequency mean gauge fields in the Abelian Higgs model are investigated in the 
presence of a scalar condensate, near the high temperature equilibrium. 

More details on the presentation can be found at the beginning of each Section. 
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Chapter 2 

General formalism for non-equilibrium 
QFT 

We describe in this Chapter different techniques apphed in the hterature to the investigation 
of the real time behaviour of quantum and classical fields. We begin with the derivation of the 
exact Schwinger- Dyson equation in a general framework. We introduce approximation schemes 
and transform the SD equations into simpler kinetic equations. Next, we demonstrate to leading 
order of the perturbative expansion with respect to the powers of the coupling the equivalence 
of the Schwinger- Dyson equations with the iterative solution of the Heisenberg equation of 
motion. The method of mode-function expansion is also presented. 

We finish by presenting the linear response theory of a classical system, regarded as an 
approximation of a high-temperature quantum system. 



2.1 Schwinger-Dyson approach 



In order to understand the mechanism by which the system approaches equilibrium one needs 
to keep track of real time processes. The most general framework for dealing with field theory 
in a non-equilibrium real time settings is the Schwinger-Keldysh, or close time path (CTP) 
formalism. For a review se for example Ref. and references therein. 

Generally we need to evaluate the expectation value of time ordered products of Heisenberg 
fields for given initial density matrix p{to): 



(T$^(ti)$^(t2)...$//(tn)) := Tr {Tp(to)$//(ti)$//(t2)...$//(tn)} . 



(2.1) 



These n-point function can be obtained formally with the help of the formula 

1 



Z[0] i6j{xi) ■ ■ ■i6j{xr. 



j=0 



(2.2) 
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CHAPTER 2. GENERAL FORMALISM FOR NON-EQUILIBRIUM QFT 



tQ=-oo 



Figure 2.1: The Schwinger-Keldysh contour. 



from the generator 



(2.3) 



If we want to evaluate this quantity perturbatively then we keep to the usual way of doing 
this by splitting the Hamiltonian into a free and an interacting part and by switching to the 
interaction picture. The different pictures coincide at time to, which means pnito) = Piito)- In 
the interaction picture the density matrix is evolved from its value in the remote past by the 
unitary time evolution operator U{tb,ta) = exp —i[J^'' dtH'^^((^j{t))] according to the relation: 



p(to) = U{to, -oo)pi(t = -oo)f/+(to, -oo). 



(2.4) 



We denote the density matrix in the remote past by po that represents the density matrix of a 
free system if we switch on the interactions adiabatically. 

Using the relation between the fields in the interaction and Heisenberg picture, namely 
$H = {t,to)^i(t)U{t,tQ), Eq.( |2.4|) the cychc invariance of the trace and the properties of 
the time evolution operator we find: 

(T$H(ti)$H(t2)...<f/f(tn)) = Tr {poTU+itf, -oo)$,(ti)$/(t2)...$/(tn)f/(t/, -oo)} . (2.5) 

We arrived at this formula by inserting the unit operator 1 = f/(t„, tj)t/(t/, t„) after $/(t„). 
The time t/ is arbitrary but its value must be above the largest time argument of the n-point 
function to be evaluated. Eventually its value will be set to +oo. 

In order to deal with the U~^{tf, — oo) (for a nonequilibrium system no state of the past can 
be related with a state in the future) the Schwinger-Keldysh contour of Fig. 1 is introduced, 
the use of which is mandatory in an off-equilibrium situations. As one can see the contour goes 
first from — oo to t/ = +oo and then back from tf = +oo to — oo. The free fields (the fields in 
the interaction picture) appearing in the interaction Hamiltonian are living on the contour. For 
physical observables the time values are on the upper branch, yet in a self-consistent formalism 
both upper and lower branches will come into play at intermediate steps of the calculation. 
This fact explains the need for four Green functions. 

The time ordering on the contour is defined through the relation: 



(2.6) 



10 



2. 1 . SCHWINGER-DYSON APPROACH 



where 0c(a^o,Z/o) = 1, if Vo precedes xq on the contour, otherways it is zero. The relation above 
also means that on the upper branch we have chronological ordering while on the lower branch 
we have anti-chronological ordering. 

Summarising we have for the time-ordered product of fields : 



(T$H(ti)$//(t2)...$/f(tn)) = Tr <; poTcCxp 
where by definition 



z / dtHUHt)) 



$/(tl)$/(t2)...$/(tn) > , (2.7) 



n poo pC 

dt= dti- 

Jc J — oo J —c 



dU. 



{2.1 



Here 1 stands for the upper and 2 for the lower branch of the contour. 
This once again can be obtained from the generating functional: 



Z\j{x)] = exp 



-t / dtHL 



int 



6 



iSjit) 



expi J d'^xj{x)^i{x) 



(2.9) 



where j^d'^x = J^dt J d^x and j{x) is defined also on the contour and has different values on 
the upper (1) and lower (2) branch of it: 



j{x) = Qcitf,xo)ji{x) - Qc{xo,tf)j2{x). 



(2.10) 



Using the Wick theorem (see [^) which is an operator identity of the form: 



expz J d'^xj{x)^j{x) 



exp ( — / d'^x / d'^x' j{x)Gc{x — x')j{x') 

c J c 

4 



Tr <^ po : exp 



(2.11) 



i d xj{x)^i{x) 



one obtains the following expression for the generating functional: 

6 



Z[j{x)] = exp 



zo[j{x)mj{x)], 



[2.12) 



ZqIJ^x)] being the first factor and Nc[j{x)] the second factor on the r.h.s. of Eg. ( |2 . 1 1|) . 
In the preceding equations, Gc{x — x') is the vacuum or causal propagator 

G,{x-x') = (0|T,$,(a;)$,(x')|0) 



(2.13) 



and Njj{xy\ describes the "initial correlations". It's interesting to see that due to the fact that 
$/ satisfies the homogeneous equation (pq — — m^)$/(p) = 0, the correlations contribute to 
the propagator only on the mass-shell, and that each component of the propagator gets the 
same additional term. 
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At equilibrium and in the non-equilibrium case at least in the kinetic approximation (see 
p7| ) Nc[j{x)] can be expressed in terms of the corresponding thermal propagator Gt 



Nc[j{x)Y'' = exp J (fx J d*x'j{x)GT{x - x')j{x') ] . (2.14) 

So we can write: 

Zo[j{x)] = exp ( - / d'^xi I d^X2j{xi)A°^{xi - X2)j{x2) ) , (2.15) 



where A°(a;i — X2) is defined by 

A°(xi -X2) = Qcixo,yo){^iix)^iiy)) + Qciyo,xo){^i{y)^iix)) (2.16) 

and consists of two parts: a T-dependent an a T-independent part A°(,ti — X2) = Gc{xi — X2) + 
Gt{xi - X2). 

Using the expression Eq. ( |2.10| ) for the current and Eq. ( p.8| ) we obtain 

^ f ^ /"OO /»oo 

dt dt'j{x)A^^{Xi-X2)jix')=J2 dt dt']r{x)^^.Xxi-X2)]s{x') (2.17) 

; Jc 1 J —oo J —oo 



r,s=l ' 



with 

A[^{x,x') = {T<!>i{x)^i{x')) = Q{t-t')A>{x,x') + Q{t' ~t)A<{x,x') (2.18) 

Al,{x,x') = A<{x,x') (2.19) 

Al^{x,x') = A>(x,x) (2.20) 

A°2(a;,x') = {T*<^i{x)<!>i{x')) = Q{t' -t)A>{x,x') + e{t-t')A<{x,x') (2.21) 

where T is the time ordering, T* is the anti time-ordering, and the Wightman functions 
Aq{x,x'), Aq{x,x') stands for two distinct ordering of $(a;) and $(x') 

A<{x,x') = (<l>/(x')$/(x)), A>{x,x') = {^i{x)^i{x')). (2.22) 

The index of the propagator denotes the branch from which the first and the second time 
argument is taken. 

Extending the definition of the time ordering on the contour (|2.6|) to fields in Heisenberg 
picture we introduce the exact Green-function of fields that possess a non-vanishing expectation 
value to be treated as a classical or "mean" field, in the form: 

zA{x,y) = {TMxMy)) - mx))my)), (2.23) 

where the H subscript was omitted. Then, depending from which branch the arguments of the 
propagator are the same relation as in Eqs. ( |2.18| ) . . . ( [^.22| ) hold. 
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The two functions defined in ( p.22| ) are related to each-other through the commutator of 
the field: 

A>{x,y) - A<{x,y) = =: {D{x,y)), (2.24) 

where for the free case 

with Uk = {m? + h?)^/"^ and ko = Uk- 

The equation of motion for the 2-point function can be obtained in the form: 



[dl + m^]A{x,y) = -6^^\x,y)+ / d^x'U{x,x')A{x' ,y), (2.26) 

Jc 

[dl + m^]A{x,y) = -5^^\x,y)+ [ d^x'A{x,x')U{x',y), (2.27) 

Jc 

where for the case of scalar theory with a given interaction Lagrangian term Cint 

^^dVn(x,x')A(x',y) = {Tc^^Hy)) - (^)("^(?/))' (2-28) 

z [ d'x'A{x,x')U{x',y) = (T,^$(x)) - (^)($(x)), (2.29) 



as can be seen doing perturbative expansion in the r.h.s. of (p.28|) , (|2.29| ). 

The effect of the mean field can be separated from the self-energy by writing: 

n(x, y) = UMF{x)si'\x, y) + n>(x, y)e,(xo, z/o) + n<(x, 2/)e,(i/o, x^) (2.30) 

where 

{(5*^^^(x — y) for xq, yo from the upper branch, 
for xo, yo from different branches, (2-31) 

— 5^^-'(x — y) for xq, yo from the lower branch, 

and Y\.mf{x) is the mean field self-energy while YP''^^\x^y) is the collisional self-energy which 
provides the collision terms in the transport equations. In the case of a theory we will see at 
the end of Subsection ^.1.2| what exactly is liMpix) and the leading expression of 11^^^^ (x,?/) 



in a perturbative expansion. 

Due to the fact that Ao(x,|/), the free Green function satisfies the equation 

[dl + m^]Ao{x,y) = -5^^\x,y), (2.32) 

we can rewrite Eqs. (|2.26|) and ( p. 27 ) in the form of a Dyson-Schwinger equation: 

A = Ao - AqHA. (2.33) 
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Using an equation analogous to ( |2.16| ) corresponding to Heisenberg fields: 

A(x, y) = ee(xo, yo)A>{x, y) + Q^yo, xq) A<(x, y), (2.34) 

Eq. ( p.8|) , and the fact that xo,yo are on the upper branch of the contour but x'q could be as 
well on the upper as on the lower branch, we obtain: 

/yo 
rfVn>«)(x,x')[A<(a;',y) - A>{x',y)] + 
-oo 

/XQ 
d''x'[U>{x,x')-U<{x,x')]A>^<\x',y), (2.35) 
-oo 

/yo 
d'x'A>^<\x,x')[U<ix',y) - n>ix',y)] + 
-oo 

/XQ 
dV[A>(x,x') - A<(x,x')]n>«)(a;',y). (2.36) 
-oo 

Introducing the retarded and advanced quantities 

nnix,y) = (n>(x,y)-n<(x,y))e(xo-i/o), (2.37) 

An{x,y) = (A>(x,y)- A<(x,y))e(a;o-2/o), (2.38) 
nAix,y) = in<ix,y)-n>{x,y))eiyo-Xo), (2.39) 
AA{x,y) = {A<{x,y)-A>{x,y))eiyo-Xo) (2.40) 

we obtain the following equations of motions: 

[dl + m'' -nMFix)]A>^<\x,y) = J [n>«)(x, x') A^(x', y) 



+Un{x,x')A>^<\x',y)] (2.41) 
[d^y+m^ -UMF{x)]A>^<\x,y) = J d^x' [A>'^<\x,x')Ua{x' ,y) 

+AR{x,x')U>^<\x',y)] (2.42) 



where the time integration runs no more on the contour C but from — oo to +oo. 

Using the defining equation for the advanced and retarded two-point function and the 
commutation relation of the fields we can express the derivative of Ajk^a) function with A>(<\ 
for example 

dlAnix,y) = {dlA>ix,y) - dlA<ix,y)) 0(xo - yo) - 6^*\xo - yo). (2.43) 
Then using Eqs. (|2.41|) and (|2.42| ) an equation of motion for the advanced and retarded two- 



point function can be derived 

[dl + m^ - Umf{x)] An^A){x,y) = -d^^^^x - y) + ^ ciVn^(^)(x, a;')A^(^)(x', y), (2.44) 

[d^y+m^-UMFiy)]AR^A)ix,y) = -6^'^\x - y) + J d^x' A r^a) (x, x')Un^A) (x' ,y). (2.45) 
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Equations ( |2.41| ), ( p.42| ), ( |2.44| ), ( p.45| ) are exact and they are equivalent to the field equa- 



tions of motion and they are known as the Kadanoff-Baym equations first derived in the 
framework of non-relativistic many-body theory . 

In order to solve these equations one has to implement an approximation scheme. Usually 
this is done in a way to transform the general equations into much simpler kinetic equations for 
the two-point function. The approximation scheme involves gradient expansion, quasi-particle 
approximation and the perturbative expansion of the self-energy. We present it in Subsection 
following the treatment of Ref . |^9|] . We can consider the classical equation of motion for 



the two-point function as an approximation of the quantum Kadanoff-Baym equations. Using 
linear response theory this is presented in Section |2.2| . 

Another approximation is to perform a perturbative expansion of the self-energy. By doing 
this UO\x,y) can be expressed with the help of the two-point function A'**^*~) closing in this 
way the equation of motion for A-^'-*--* (note that A^(x, y) = A^{y,x)). 

2.1.1 Kinetic equations for the two-point functions 

In thermal equilibrium the two-point functions depend only on the relative coordinates = 
~ Vfi and are strongly peaked around = with a range of variation determined by the 
wavelength of a particle with a typical momentum k in the plasma. At high temperature k ^ T 
and the particle's thermal wavelength is Ay = ~ T. 

Out of equilibrium the two-point function depends on both coordinates, and if the deviations 
from equilibrium are slowly varying in space and time it pays out to introduce two new variables: 
the relative coordinate u^ = x^ — and the center-of-mass coordinate = ^h^Mr terms of 
which the kinetic equation can be derived. For slowly varying disturbances with A ^ At, one 
can expect that the u dependence of the two-point function is close to that of equilibrium and 
(9„ ~ A; ~ T and (9x ~ 1/A < T. 

For a function f{X,u) that varies slowly with X and is strongly peaked for u ^ one can 
approximate f{X + u, u) as 

/(X + u,u)^ f{X, u) + u^ ^^'^^^\ (2.46) 

i.e. only terms involving at most one derivative dx are kept. This approximation is referred to 
as the gradient expansion. 

Then one defines the Wigner transform 

A>(x,y) A>(X,p) := j rfwV^-A>(X + - ^) (2.47) 

A>(a:,2/) = j (|^e-^^^^-^)A>(^,p). (2.48) 
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The properties of the Wigner transformation, useful for our calculation, are the following: 



a{x)f{x,y) 



W.tr fl 
W.tr 



d'yA>ix,y)'^' [^dx+ip] A>iX,p), 



P) - i 



I dajX) df{X,p) 
2 dXi" dp^ 



(2.49) 
(2.50) 
, (2.51) 



aiy)n.,y) / |^a(.)/ + I) e-- . a(X)/(X,p) + 4^^, (2-52) 



d^zf{x,z)g{z,y) f{X,p)g{X,p) + ^ {f,g}pB. + 



(2.53) 



In the last formula (for the derivation of it we refer to Section 4.2. of Ref. [^) {f,g}pB 
denotes a Poisson bracket 



{/' 9}p.B. ■= 9pf ■ 9x9 - dxf ■ dpg, 



(2.54) 



and the dots mean that the gradient expansion was used. The first expressions on the r.h.s 
of Eqs. ( p.51| ), ( 2.52|) represent the exact result of the Wigner transform, while the second 
expressions the result from the gradient expansion. 

Of course, for other two-point functions and the self energies analogous relations hold. 

The kinetic equations for the two-point function A^*^^) is obtained by subtracting the Wigner 
transform of Eqs. ( 2.41 ) and ( 2.42| ) while for the two-point function Aj^^a) by adding the Wigner 
transform of Eqs. ( 2.44 ) and ( 2.45|) . 

Using the properties 



Ar{X,p) - Aa{X,p) = A>{X,p) - A<(X,p), 
Un{X,p) - Ua{X,p) = U>{X,p) - U<{X,p), 



(2.55) 
(2.56) 



that follows from the definition ( p.37|) . . . (|2.40| ) and 

An{X,p)+AAiX,p) = 2Re A^(X,p), n«(X,p) + n^(X,p) =2ReUR{X,p), (2.57) 



one obtains the following kinetic equation 
dReU\ A>«) dUMF 5A>«) 



2 ho^ + 



dp. 



dX. dp^' 



{ Re A^, n>«)}^^ = I (n<A> - A<n>) , (2.5^ 



where Ren(X,p) := Hmf + Ren^(X,p). 

This equation represents the quantum generalisation of the Boltzmann equation. The terms 



of this equation have the following physical interpretation (see Ref. |39]). The Wigner function 
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A plays the role of the phase-space distribution function. The drift term on the l.h.s. generalises 
the kinetic drift term by including two types of self energy corrections. One is due to the real 
part of the self-energy that acts as an effective potential whose space-time derivative provides 
the "force" (9^ Re 11) (5^ A*-). The second correction is due to the momentum dependence of 
the self-energy that modifies the velocity of the particles. The terms on the r.h.s. describe 
collisions, while the term containing the Poisson bracket accounts for the off-equilibrium shape 
of the spectral density as one can see in Eq. ( p.59|) . 

For the spectral density p{X,p) = iA^{X,p) — iA^{X,p) we obtain 

2 fp- + ^) ^ - ^ * = { Re A„, n> - n<}p, (2.59) 

One arrives at the quasi-particle approximation when the term on the r.h.s. is neglected. 
In this case the spectral density solution of Eq. ( p.59| ) read as 

p(X, p) = 27re(po)5 (p' - + Re U{X, p)) . (2.60) 

In a recent work by Aarts and Berges |^ non-equilibrium time evolution of the spectral 



function was studied numerically. There was observed that the spectral function develops a 
non- vanishing width. This points towards the need of relaxing the assumption of zero-width of 
the spectral density if a consistent quantum-Boltzmann equation is aimed. 

In the mean field approximation, that is when the collision term is neglected altogether we 
obtain: 



p.dx + \{d^X^MF)dl 



A<(X,p) = 0, (2.61) 



p(X, p) = 2neipo)5 {p^ - + Re Umf{X, p)) . (2.62) 

Neglecting the collision term corresponds in the Dyson-Schwinger formalism to the trun- 
cation of the hierarchy of the n-point function at the level of the 2-point function, all higher 
n-point functions are neglected. 

We will use this approximation in the thesis. 

2.1.2 Iterative solution of the Heisenberg equations 

We can arrive at the approximated Kadanoff-Baym equations with a different method too. 
We start with the equation of motion of the full quantum operator 

(92 + m2)|.(x,t) + ^l>=^(x,t) = 0, (2.63) 
6 

then we split the field into the sum of its average (f){x) and the quantum fluctuation around it 

<l(x, t) = 0(x, t) + y.(x, t), 0(x, t) = (<l>(x, t)). (2.64) 
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Upon averaging the Heisenberg equation one obtains the equation of the classical field 

and subtracting this from the original Heisenberg equation we obtain the equation from the 

fiuctuation ip{x). These two equations reads as: 



(^52 ^ ^2 ^ A^2(^) ^ ^{^\x)))<l>ix) = -^{ip'ix)), (2.65) 
d'+m' + + ^{^\x)))^ix) = [^\x)] - ^ [^'(x)] + ^{^'ix)Mx), (2.66) 



where [■] = ■ — (■)• We have added to both sides of the equation of (f{x) the term ^{Lp'^{x))ip{x). 
This will mean a resummation and provides the mean-field part of the self-energy, the IImf{x) 
when considering the equation of motion for {ip{x)ip{y)). 
Introducing the notations: 

j{x) := -i0(x) [^\x)] - \ [p\x)\ + \{^\x))^{x), (2.67) 

we follow Ref. and solve the equation 

{d^ + ^2 ^ A^2^^) ^ j^(^) ^ (2.68) 

in the form 

Lp{x) = ipo{x) - X dx'GR{x,x')j{x'), (2.69) 



with ipo the solution of the homogeneous equation 

+ rn' + ^<p\x) + ^{^\x)))m^) = 0, (2.70) 
and Or the retarded Green-function associated with this equation 

[d' + m' + ^02(x) + ^{^'{x)))GRix,y) = -5^'\x,y). 

Using the homogeneous equation and the commutator relation of the free field ^po{x) it's 
easy to see that: 

Gr{x, y) = {G>{x, y) - G<(x, y))Q{xo - yo), iG>{x, y) = (<^o(x)<^o(y)), (2.71) 

iG<{x,y) = {ipo{y)^,{x)). (2.72) 

For the equation of we calculate to hnear order in A the quantities {ip^i^x)) and {ip^i^x)) 
using equation ( p.69| ) and obtain: 

{^\x)) = {^l{x)) + 0{\% (2.73) 
{ip\x)) = {iflix)) + \ j d^x'(t){x')GR{x,x')[G>{x\xf 

+G>{x,x')G>{x\x) + G>{x,x'f] +0{\^). (2.74) 
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The equation for {(p{x)!f{y)) is obtained after multiplying Eq. ( |2.68| ) with (p{y) and taking 
the average. We need to evaluate 

xijixMy)) = -^mi^'i^My)) - ^i^HxMy)) + ^{^'ix)){^ixMy)), (2.75) 

making use of Eq. (|2.69| ). The 0{X) contribution is zero. 

Taking into account that {j[ipo{x')]ipo{x)) = up to (9(A^) we obtain: 

MjixMy)) = -X' I dy'Gn{y,y'){j[M^MMy')]) 
X^ f 

/ dx' Gr{x, x') {{j[(po{x')]ipo{x)(poiy)) + {Mx)j[Mx')]My))) 

X^ f f 

/ o?a;'G/j(x,x')(^(jbo(a;')]v^o(a;)^o(y)) + (<^o(a;)jbo(a;')]<^o(a;)v?o(2/)) 



+ (</'o(a;)jbo(a;')]</'o(y)) )• 



(2.76) 

The averages on the r.h.s. can be evaluated with the help of the Wick theorem leading to: 
-— / c/xV(x)0(x')[Gjj,(i/,x')G>(x,x')^ + G^j(x,x')G>(x',y)(^G>(x',x) 

-y / {GR{y.AG'^{x,x^f ^Gr{x,x')G'>{x^ ,y) (G'>(x',x)2 

+G>(x,a;')G>(a;',x) + G'>(x,a;')^)] . 

We observe the same structure in the last line of the above equation as the one appearing 
in the Eq. ( p.74| ). Using the definition ( p.71|) , the property G^{x,y) = G^{y,x) and some 
algebraic manipulation with G^'^^^ one arrives to express all terms with some power of G^^^\ 
One obtain: 



X{j{x)f{y)) 



X' 



X 



yo 



x) / d^x'G>{x,x'y G<{x',y)-G>{x',y) (j){x') 



— oo 

XO 



--(f){x) / rfV [G>{x,x'y - G<{x,x'f] G>{x',y)(f){x') 

2 J-oo 

/ d^x'G>{x,x')^ G<{x',y)-G>{x',y) 

r d^x' \G''{x,xy -G<{x,x'f]G>{x',y) + 0{X^). (2.77) 

Now we can replace the Green's function iG^{x,y) = {ipo{x)ipo{y)) with the exact one since 
the corrections will be 0{X^) implying (9(A'^) corrections in the equation above. 

Introducing for the one and two-loop self-energy components that emerges in a perturbation 
expansion in the broken phase the notation 



n 



>«), 

2a ' 



x,x') := -^0(x)G'>«)(a;,x')'0(a;'), U^i^\x,x') := -^G>^<\x, x')\ (2.78) 
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a) b) 

Figure 2.2: One-loop a) and two-loop b) self energies. 

up to order the self-energy takes the form: 



(2.79) 



The components of the self-energy labelled with a and b correspond to the graphs of Fig. |2]2 . 

With the help of the above relations we can write down to order the following set of 
equations: 

2 r^o 



A A \ A A"^ 

d' + m' + -02 + -{^^)yix) + -03(x) + -j dx' (n> (x, x') - n< (x, x')) 0(x') = 0,(2.80) 

d'x'U>{x,x') G<ix',y)-G>ix',y) 

-oo 

d^x'[U>{x,x') -U<{x,x')]G>{x',y). (2.81) 

-oo 



d' + m' + ^<l>'ix) + ^{^'ix)))G>ix,y) 



We see that Eq. ( g^SlD for the {(p{x)(p{y)) has the same structure as Eq. (p.35|) . Actually 
they coincide if A^'-^-' is identify with G-^^^^ and —]1mf{x) with ^^^(x) + ^{ip'^{x)) and if one 
expand the self-energy in Eq. (|2.35|) up to A^. 

If one express the self-energy 11^'^'^^ in term of G^ alone the equation for turns into an 
equation with the same structural form as the one obtained from the three-loop 2PI effective 
action by Berges and Cox |T^, although the relation between the two approximation schemes is 
not yet clear. With the help of this equation these authors has shown in a numerical simulation 
in 1 + 1 dimensions the thermalization of a quantum system. 



2.1.3 Mode-function expansion 

The method of mode-function expansion was used in the investigation of the evolution to- 
wards equilibrium of non-equilibrium states produced during preheating. This was done in the 
Hartee-Fock approximation in which the dynamics is described in terms of a mean-field and the 
fluctuations around it are characterised by the two-point function only. The two-point function 
can be constructed in terms of a complete set of mode functions that, if the mean field is homo- 
geneous, can be chosen as plane waves labeled by the wave vector k. During preheating, only 
mode functions in a narrow k-band will be excited by the time-dependent homogeneous mean 
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field, via the phenomenon called parametric resonance. In an early stage, the so called linear 
stage, when the occupation number of different modes, and consequently the fluctuation is not 
too large the dynamics of the system is well approximated by the equations of motion linearised 
with respect to the fluctuation. In the second stage, called back-reaction the fluctuation grows 
exponentially and one has to switch to a fully non-linear dynamics (see for example 0). 

The system eventually becomes stationary but due to the lack of scattering the method of 
the mode function expansion in a homogeneous background cannot lead to thermalization. This 
is reflected by the shape of the distribution function that does not approach the Bose-Einstein 
distribution in case of bosons showing instead resonant peaks . 

Recently it has been proved that considering inhomogeneous background for which the 
quantum modes that are represented by the mode function can scatter via their back-reaction 
on the background, the time evolution of particle distribution functions display for a flnite time 
interval approximate thermalization |^ . 

The mode-function expansion on an inhomogeneous background is needed when the density 
matrix describing the initial state of the system does not commute with the translation. In 
this case the quantum average will depend on the space-time position and the mean fleld is 
inhomogeneous . 

The method of mode-function expansion represents a way of solving the Eq. ( p. 70 ) in which 
the expression (y9Q(x, t)) need to be renormalised. We do so by subtracting its value at t = 0. 
Introducing 

m'=ml{T) + ^{0li^,O)) (2.82) 

we obtain 

d' + ml{T) + ^02(x,t) + ^MM) - (y.2(x,0))]]y.o(x,t) = 0. (2.83) 

The solution of this equation can be constructed by expanding ipoi^x, t) into a series with respect 
of an appropriately chosen complete set of time dependent functions, named mode-function with 
constant operatorial coefficients, labelled by "k": 

M^,t) = ^(afcV^fc(x,t) + a+V^^(x,t)), (2.84) 

k 

7r(x,t) = ^(afcV'fc(x,t) + a+V'fe(x,t)), (2.85) 

k 

where the creation and annihilation operators obey (by definition) the usual commutation rules: 

[afc,a^,] = 6k,k'- (2.86) 
The standard canonical commutator between y9(x, t) and 7r(x, t) requires the relation: 

J2 [M^, 0)V'fc(y, 0) - ^*(x, 0)My, 0)] = ^5(x - y), (2.87) 
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which is fulfilled by the choice: 



M^, 0) = ^,(x, 0) = ^e-^^\ (2.88) 

Here we choose 

u;l = k' + ml{T) + ^, = ^ / d'x<P{^, 0). (2.89) 

The T appearing in this relation is the initial temperature, fixed as part of the initial physical 
data of the system. 

Then an initial state corresponding to thermal equilibrium would be characterised by the 
number operator expectation value 

(4«fc)o = ^ = nM, (2.90) 

and all the higher moments expressed in terms of this. Using the thermal density matrix for 
the initial moment we find for the relevant field expectation values 

(y.2(x, t)) = \M^, t)\\2n{ujk) + 1), {^H^, t)) = 0. (2.91) 

k 

In the evolution of the order parameter, the contribution from the second moment at t = 
is absorbed partly into the renormalised squared mass, partly it contributes the temperature 
dependence of the mass. (This is the same renormalization we applied in the equation of (po). 
In this way we find the effective equation for the inhomogeneous mean field: 

d^ + mliT) + (|^fc(x,t)|2 - |^fc(x,0)|2) {2n{ujk) + l)]0(x,t) + ^(j)^{^,t) = h. (2.92) 

k 

The equations for the fc-modes are found from Eq. (|2.83|) when multiplying it by aj, and 
taking its expectation value: 

d^ + ml{T) + ^0'(x,t) + (|V^fc(x,t)|2 - |^fc(x,0)|2) {n{uk) + ^)]^/'p(x,t) = 0. (2.93) 

k 



This equation is the generalisation of the equations proposed in ^ for homogeneous 
classical background. The method of mode function expansion in an inhomogeneous background 



was applied in numerical simulations in Ref. [44 



2.2 Classical linear response theory 

We present now the classical linear response theory of Jakovac and Buchmiiller following closely 
their paper [4?]. Using this method, we can calculate in linear approximation the damping rate 
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of a field. When we compare in Subsection ^4.1.4] the analytical quantum damping rate with 
the classical one or the numerical value for the damping rate with the analytical one as we do 
in in Section we need the result for the classical damping rate in the broken phase. Here 
we present for simplicity the derivation of it for the symmetric phase. The extension of this 
method to a classical 0{N) scalar theory in the broken phase can be found in Appendix 0. 
In a classical theory at finite temperature T = (3~^ defined with the Hamiltonian 

H = j ^^^(Ivr^ + \{d<pf + \ml{k)<i? + ^Ki<P' + 3<p) , (2.94) 

that contains j{t,x) as an external source, the classical real time n-point functions are defined 
in the following way 

(0(ti, xi) . . . 0(t2, X2)}ci = ^ J P7rP0e-'^^('^'<^V(ti, xi)... <P{h, X2) , (2.95) 
with the partition function 

Z = j pTTp^e-^^^''^) . (2.96) 
In Eq. |2.95| (/)(t, x) is the solution of the equations of motion 

0(t,f) = 7r(t,f), (2.97) 
7T(t,x) = A(l){t,x)-ml{A)(l){t,x)-^(t)\t,x)-3{t,x), 

with the initial conditions 

0(to, x) = 0(x) , 7r(to, x) = 7r(x) . (2.98) 



In Eqs. ( ^.951 ) and ( ^.96] ) the integration is performed at t = to over the space of initial conditions 



= 0(a;),7r = 7r(x) (averaging over the initial conditions). 
As it is well know from the time prior to the discovery of the quantum theory, a classical 



theory is not well defined in the UV. As was shown in |T3, |2J] for a scalar classical theory of 
constant temperature T both for equal time ( or equivalently for static quantities) and unequal 
time correlation functions merely a mass renormalization suffices to render the theory finite. In 
Eq. ( p. 94 ) the mass term can be split into the sum of two terms: the renormalized classical mass 
and the divergent part of the mass, to be treated as a counterterm: m^;(A) = m\ + 5m? {A). 
Here A is an UV cutoff introduced to regularize the theory, and 5m? {h) is the sum of two terms, 
(5mi(A) and (5m2(A), corresponding to the one and two loop divergent diagrams, as we will see 
at the end of this Subsection. 

If the classical theory is considered to be the high-temperature approximation of the quan- 
tum theory, then the renormalized classical mass should be related to the tree-level mass of 
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the quantum theory. Such a relation between masses can be derived using the dimensional 
reduction by matching the result of some quantities calculated in both quantum and classical 
theories. For the exact form of the matching relations (the masses, the coupling constant and 
the expressions of 6mi{A) and 5m2(A)) I refer to Refs. ||22|, [17| . 

The j(t, x)-dependent solution of the equations of motion constructed with the retarded 
Green's function 

DR{t,x) = [ -^e'^^eit)-sm{uj^t) . (2.99) 
satisfies the integral equation {t,t' > to,x = (t,x)), 

0(x;j) =0o(x) + J d^x'Dnix - x')(^^(l)^ix';j) + jix')y (2.100) 

00 is the solution of the free (homogeneous) equations of motion which satisfies the boundary 
condition. 

As one can see from equation (|2.100| ) the classical solution depends on the external source 
in powers of which it can be expressed. The idea of linear response theory is to keep only the 
term linear in j{x). The linear connection between the solution of the equation of motion and 
the external source is given by the retarded response function defined through 

(2.101) 

j=0 



Hr{x-x') ""^"^'^^ 



5i{x' 

in the form 



j d^x']{x')HR{x-x'). (2.102) 



The retarded response function depends on the initial conditions 0(x) and 7r(x). Substitut- 
ing Eqs. ( p. loop into ( |2.101| ) one easily gets the integral equation which determines Hr, 



Hr{x - x') = Dr{x - x') + ^Aei j d^yDRix - y)(P\y; 0)HR{y - x') . (2.103) 
This equation can be solved iteratively, for example up to order A^; its solution reads: 
H^^\x - x') = Dr{x - x') + J d'yDRix - y)^\y- 0)DR{y - x') + (2.104) 
^Xlj d'y j d'zDRix-y)(P\y;0)DR{y-z)<p\z- 0)Dr{z - x') + 0{Xl), 



and is depicted in Fig. ^ 



In order to obtain the finite-temperature retarded response function denoted with D^, the 



ensemble average of Eq. (|2.100|) with respect to the initial conditions is taken. This yields 

{<P{x))ci = J d^x'j{x')Di{x - x') , (2.105) 
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Figure 2.3: The response function Hn expanded in powers of the classical field. Full lines 



denote retarded Green's functions Gr. The figure is from Ref. p7 
where 



- x') = {Hr{x - x'))a = V(t>VTie-f'"HR{x - x') . (2.106) 



can be evaluated by first expanding the solution ( p.l03| ) for in powers of Xd (cf. Fig. |2] 
and then performing the thermal average for each term. 
This method means that we have to evaluate 

Di{x - x') = {e'^^'^-^Hitix - x'))o,c, (2.107) 



using the iterative solution of Eq. (|2.103|) . In the equation above Hint is the interaction 



Hamiltonian, the subscript means thermal average taken with the free Hamiltonian, and c 
means that we have to take only connected graphs. 

If the classical theory of interest is only a tool to evaluate the leading contribution of a 
quantum theory then, the full result of the classical theory is irrelevant and one can make the 
approximation 

{e-P^^-^Hn{x - x'))o,c ~ {Hr{x - x'))o,c (2.108) 

because only to leading order in the coupling the results of the classical and quantum theory 
are related in the high-temperature limit. 

So, what is required is the evaluation of the thermal n-point functions of the type of (|2.95|) , 
namely 

(<^^(xi) . . . 4>\x2n))ci = 2"5^(a;i - X2) • • • S\x2n-i - X2n) + permutatlous . (2.109) 
where {t — t',x — x') is the free two-point function that reads p4|, 



S^t-t',x-a') = {<P{x)4>{x'))',, = T I ^y^(^-^-)l^^cos{Mt-t')). (2.110) 



From Fig. p.3| it can be easily seen that the finite temperature response function obtained 



after thermal averaging with Hq satisfies a Dyson-Schwinger equation 

Diix - x') = Dr{x + j j d^y'DRix - y)n,i{y - y')D<^J,{y' - x') . (2.111) 
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HJix) is the classical self energy and by looking at the Eq. ( |2.104| ) we see that up to 0{Xli) 
in perturbation theory it is given by: 



ni;^(y-y') = {<P'\y,0)) = '-^S{y-y')6{y-y') 



d 



(2.112) 



^\y, m\y\ 0)) = -fS{y - y')Dn{y - y')S{y - y'). (2.113) 



These two terms are shown also in Fig. ( |2.4| 



(a) 



(b) 



Figure 2.4: Perturbative expansion for the self-energy 11. Dashed lines denote thermal two-point 
functions S. The figure is from Ref. 



The contribution 0{\ci) (Fig. |2.4| a) is linearly divergent and the divergence can be removed 



by a mass renormalization yielding the counter term 



5ml = lAa5(0, 0) = h^iT ! . (2.114) 



2 " ' ' ' 2 " 7 (27r)3tu, 



The contribution 0(A^J (Fig. |2.4| b) is logarithmically divergent and can be rendered finite by 
a further mass renormalization 

5ml = -\>^liP [ d'xS'iO, x) = -\\lT^ [ d^O)^- , (2.115) 
b J S J UJ1UJ2UJ3 

where d^{p) is the integration measure, 

Contributions of higher order in Ac; are all finite. 

The damping rate, i.e., the imaginary part of the self-energy is finite. Using Eqs. ( p.99| ) and 
( |2.11CI| ) one easily obtains for the self-energy to order A^;, 



U{uj,p)^i = j dtd^x e'^^'e-'P^ U{t,x)ci (2.117) 



d^{p) 
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where = ±1, i = 1, . . . , 3. 

From Eq. ( |2.117| ) one obtains for the imaginary part of the self-energy {u > 0), 

IT f 1 

T{u,p)ci = -XliT^ / d^{p) [uJi6{u -UJ1-U2- UJ3) + u6{uj + ui - U2 - u^)] . (2.118) 

Z J UJ1UJ2UJ3 

When comparing the classical damping rate Td with the damping rate F of the full quantum 
theory it turns out that at very high temperatures, i.e., T 3> m, where the Bose-Einstein 
distribution function becomes 

f{u) ~ - > 1 (2.119) 

00 

to leading order in the coupling the damping rates are given by the classical theory. 
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Chapter 3 

One-component theory 



The purpose of this chapter is to develop and test the methods of investigation of the evolution 
of out-of-equilibrium systems. 

We start in the first Section with a new approach, a kinetic theory, aimed to substitute the 
quantum degrees of freedom, coupled to a classical field theory. This approach can be regarded 
to be at an intermediate level between a fully quantum treatment and its high-temperature 
classical approximation. This method of treating the short wavelength quantum modes has 
been used in gauge theories in numerical works aiming to measure the hot sphaleron rate, when 
it was realised that quantum degrees of freedom play an important role in the determination 
of this rate. 

We present the connection of this way of treating the quantum modes with the Schwinger- 
Dyson approach in linear response approximation. 

In the second Section we present a numerical investigation of the non-linear dynamics in 
a classical $^ theory far from equilibrium. The system starts from a metastable state and its 
evolution is followed as it undergoes a first order phase-transition to the true vacuum. 



3.1 Classical kinetic theory of Landau Damping 

Relevance of combined classical field theory and kinetic theory to the high temperature linear 
response theory of quantum gauge fields has been demonstrated by several authors ||9, 5C . 
These authors have shown that for high enough temperature the leading non-equilibrium trans- 
port effects of large wave number |k| ^ T fluctuations can be reproduced by a kinetic theory 
of the corresponding response functions. This approach has been applied to the study of the 
QCD plasma |^ and proved successful in reproducing the contribution of hard loops to 
the Green's functions of low |k| ^ T modes in Abelian and non-Abelian gauge theories. 

For gauge fields the central object which allows explicit construction of the kinetic Boltz- 
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mann equation is the Wong-equation |Q , since it determines the force exerted on non-Abehan 
charged particles by external fields. 

For selfinteracting scalar fields the notion of the external field and of the force felt by the 
modes with low wave number is not so obvious. The physical idea based on the renormalisation 
group is to represent the high wave number fiuctuations in form of a gas of massive particles, 
and study the effect of the gas on the dynamics of the slow modes. The derivation of the 
Boltzmann equation for the gas in the low wave number background field was presented for 



the scalar fields in |5J] using a non-relativistic potential picture. The authors have found the 
effective action which accounts for the loop contribution of high-/c fluctuations to the Green's 
functions of the low-k modes. 

The authors of Ref . |^ have assumed positive curvature for the potential felt by the effective 
particles at the origin. The above mentioned gauge investigations have also assumed no scalar 
background flelds and the restoration of all symmetries. However, certain scalar models with 
speciflc internal symmetries are known where symmetry is not restored up to high temperatures 
[56[| and the question of the derivation of a kinetic equation is still of interest. If the 
scalar theory is thought to be part of a non-Abelian gauge+Higgs system, the non-perturbative 



treatment of the coupled one-point and two-point Schwinger-Dyson equations leads to a 
high temperature solution with non-vanishing vacuum expectation value for the scalar fleld. 
Therefore it is of interest to explore the physical consequences of the presence of a constant 
non-vanishing background from the point of view of the kinetic behaviour of the effective gas. 

Motivated by the success of the classical kinetic theory in describing the effect of the hard 
quantum modes on the dynamics of the soft (low-k) ones in the context of a gauge theory in 
the symmetric phase, we present in this Section a fully relativistic Lagrangian formalism for 
the kinetic theory of scalar flelds for the case of non-zero spatial fleld average. 

The construction starts by proposing a Lagrangian for the effective particles, representing 
the high-k modes on the low-k fluctuations. An advantage of this proposition is that it leads 
to the induced source density of the low-k fluctuations directly, without any reference to the 
quantum theory. Evidence for the correctness of the effective Lagrangian can be presented 
by comparing its consequences with the results of the corresponding quantum calculations. 
In Subsection |3.1.1| a fully relativistic rederivation of the results of Ref. ||5^ is presented for 
the kinetic theory of the flelds with large spatial momenta. The emerging correction to the 
Bose-Einstein phase-space distribution function of high-k modes is used in Subsection |3.1.2| 
to compute corrections to the classical fleld equations of the low momentum modes. Finally 
the Landau damping coefficient for the off-shell scalar fluctuations is computed. This effect is 
present only in the broken symmetry phase of the scalar theory. 
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3.1.1 Kinetic theory of particles with background field dependent 
mass 

The effective gas of high-frequency fluctuations is kicked out of thermal equilibrium if an in- 
homogeneous low frequency background fluctuation is present. This state of the gas induces 
a source term into the wave equation of the low-k modes. In a scalar $^-theory a unified 
description of the two effects can be attempted if in addition to the action Sd describing the 
classical low-frequency dynamics an appropriate action can be introduced for a gas particle 
(that represents the high frequency quantum modes) coupled to the low-frequency field along 
its trajectory C,^{t). 

We split the $-field into two terms: 

^x) = / j^,me-"'-+ -— ^$(fc)e-'^- = 0(x,A) + y.(x,A), (3.1) 

Jo [^T^r Ja i^'^) 

and assuming the spontaneous generation of a non-zero constant average background field (p 
for the slow modes, we shift the right hand side of the above expression by 

$(x) = + </)(x) + (^(x). (3.2) 

The wave equation for (j){x) will be modified by an induced source density, determined as a 
statistical average over the phase space distribution of the gas particles representing the fast 
modes. 

In this Subsection we give the expression of the force felt by the effective particles repre- 
senting the high-/c modes in a background. For this we write down the approximate expression 
of the Lagrangian quadratic in the fast mode (p 

We stress that the is of negative sign. Leaving out fro Eq. (|3.3| ) the selfinteraction of the cp 
is at the basis of replacing the parts containing if fields in the Lagrangian with a gas of massive 
particle interacting only with the background. 

We propose an effective action describing the system consisting of a massive gas of particles 
interacting with a scalar background: 

Seff = SM Mi^S, mir))]dT, (3.4) 

where Sci[(f)] is simply the action of a 0^-theory understood with a cut-off A. The second term 
is the action for a gas of relativistic scalar particles superimposed on the background, (p is the 
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constant scalar condensate, and is the long wavelength fluctuation amplitude on the top of 
the average field. The local mass is determined by the field values along the path C,i{T) of the 
i-th particle: 

ML[^,mr))]=m' + ^i^ + <j){Ur))r, (3.5) 

where is the squared mass parameter of the theory, which can be negative. 

Variation of the second piece of the action in Eq. ( p.4|) with respect to the field variable (f){x) 
should yield the induced source term to the wave equation of the low frequency modes. This 
induced source can be expressed with the distribution of the gas particles. The distribution 
can be derived from the solution of a Boltzmann equation describing the gas in the background 
field 0. 

Variation of Eq. (|3.4| ) with respect to the particle trajectory provides the EOM of the 
particles together with the expression of the force exerted on the particles by the external field 
(p. This information is used in deriving the coUisionless Boltzmann-equation. 

Concentrating only on one particle, with the usual variational procedure |Q one can derive 



the canonical momentum from the translational invariance of the action, and the variable 
rest mass M from the invariance of the action with respect to the variation of the proper time 
r: 



M 



^£ -L 



(3.6) 
(3.7) 



The kinetic momentum is defined through the relation = M^^ and the force is found from 



the relation F., 



dr 



The equation of motion of the particle is: 



d_ 
dr 



dL_ 

Wu 



where L is the sum of the free particle Lagrangian and the interaction Lagrangian. 

In our case the kinetic momentum coincides with the canonical one and the EOM can be 
written in the form: 

(3.9) 



dp, 1 OMUO 



The coUisionless Boltzmann-equation arises by the application of the chain rule to the proper- 
time derivative of the one-particle phase space distribution: 



dn{x,p) 
dr 



p-d+ MioS,(P{x)] 



dp, d 
dr dp,- 



n{x,p) = 0. 



(3.10) 
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{p ■ d)n{x,p) + -(0 + (^{x))d^4>{x) = 0. (3.11) 



In the next Subsection we will show that a statistical average introduce a distribution function 
for relativistic particles. Appendix ^ shows how the distribution function is introduced in the 
context of relativistic kinetic theory. 

Exploiting Eq. (|3.9|) one writes the coUisionless Boltzmann-equation for the gas of these 
particles in the form: 

Here one follows the usual perturbative procedure to find corrections in the weak coupling 
limit A ^ 1 to the equilibrium Bose-Einstein distribution by writing the distribution function 
as a series in powers of A: 

n{x,p)=no{p) + Xni{(f){x),p) + ..., ^o(p) = ^^^^^ _ ^ Q(po - A). (3.12) 

We have slightly modified the equilibrium Bose-Einstein distribution function that depends 
only on the pq, in order to account for the fact, that the particles represents the degrees of 
freedom above the A-scale. 

The formal solution for ni{x,p) := ni{(j){x),p) which is linear in (j){x) is easily found to be 

n.ix,p) = 4T^dM-)'-¥^- (3.13) 



2{p-d) ' dp 



In the next subsection we shall calculate corrections to the linearised field equation for the 
long wavelength fluctuations. These corrections arise from the statistically averaged source 
term emerging through the ^-dependence of the particle action ( |3.4| ). This completes the 
selfconsistent construction of the wave equation for fields interacting with a relativistic gas of 
particles. 

3.1.2 The wave equation of the slow modes 

The variation with respect to of Eq. ( p.4| ) provides the wave equation of the \ow-k modes 

9^ + + -0 ) (f)(x) + m^(j) 
6 / 



It contains the source term arising from the gas particles through the 0-dependence of their 
local mass. We arrived at this expression keeping only terms linear in (p. 

In order to express the source term with the help of the one-particle distribution function 
n{x,p), we perform a statistical average over the full momentum space and in a small volume 
around the point x. This procedure introduces n{x,p) through the relation: 

d^p 



^eff 



I (2^''^'^'^^ = ^ / ^^E^^'Hx-e.(r))). (3.15) 



33 



CHAPTER 3. ONE-COMPONENT $4 THEORY 



Here we introduced "^e// = ""^^ + f "average" or "effective" mass square of a particle. 

With the help of it we approximate the dispersion relation of a particle as Pq ~ "^e// + P^? that 
is the effect of the long wavelength fluctuation is taken into account only in the force felt by 
the particle. Also, we keep only the 0{X^) term in the expression arising in Eq. (|3.14|) 



when the averaging procedure mentioned before is performed. These two approximations has 
no effect on the non-local part of the induced current. We get the following "macroscopic" 
source density: 

Jav{x) = ^J (2^)3pp ('^ + (3.16) 

The induced linear response to 0(x) is found using the perturbative solution of the Boltz- 
mann-equation presented in the previous subsection by retaining terms which produce linear 
functional dependence of jav on (p. One finds for the leading terms of the induced source the 
following expression 

Ja.{x) = ^0 / j/^MPo)+J^J{x)+ji]^{x), (3.17) 
2 J (27r)^po 



with 



i!^Kx) = h{x) I j^f^MPo), and j^{x) = ^ / j/^nMx),p). (3.18) 
2 J (27r)^po 2 J (27r)^po 



The first term in Eq. (|3.17| ) which is proportional to the average field enters in the 



equation that determines the expectation value of the (p 

+ + 7^Mpo) = 0. (3.19) 
6 2 J {2TTypo 

This contribution to the thermal mass in the high-T limit gives the correct limiting value 
for the one-component scalar theory, but it can not account for the T = renormalisation. 
The damping effect arises from the piece of current proportional with ni{x,p) 



iS'W = p^^*^W^. (3.20) 

The imaginary part of the linear response function E defined through j2u{x) = T,(j){x) is eval- 
uated in Appendix ^ by the use of the principal value theorem, when the e-prescription of 
Landau is applied to the {p ■ d)~^ operator on the right hand side of (|3.20|) . For an explicit 
expression one may assume that represents an off-mass-shell fluctuation characterised by the 
4- vector: (fco, k), that is 0(x) = (f){k ■ x). 
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3.1.3 Connection with the formaUsm of Danielewicz and Mrow- 
czynski 

Before explicitly evaluating the linear response function, let us make connection to the classical 
kinetic theory for self-interacting scalar fields derived first by Danielewicz and Mrowczyhski 
|60| , and presented in the introduction and to present an alternative but equivalent method 
for the derivation of the effect of the quantum (high-k) modes on the EOM of classical (low-k) 
modes. 

We sketch the steps followed in the derivation of the effective equations for the propagation 
of a non-thermal signal on a thermalized background: 

1. Decomposition of the fields (pi into the sum of high-frequency, thermalized {(fi) and 
low-frequency, non-thermal {(pi) components; 

2. Derivation of the equations of motion for (pi in the background of the thermal components; 

3. Averaging the equations over the thermal background, retaining only contributions from 
the two-point functions of the thermalized fields; 

4. Approximating the two-point functions resulting from step 3. by expressions with at most 
linear functional dependence on the non-thermal fields, what is sufficient for the calculation of 
the linear response function of the theory: 

S{ip{x)(p{y)) 



{ip{x)ip{y)) ^ {vix)ip{y))\^^Q+ I dz 



' </>=o 

= + MxMy))^'\ (3.21) 

We start from the Heisenberg equation of motion for the quantum field 

(a^ + m2)$ + -<l>3(x) = 0, (3.22) 
6 

then we split $(x) into the sum of terms with low {po < A) and high {po > A) frequency 
Fourier- components, that is ^(x) = (p{x) + (p{x), {(p{x)) = 0. The classical equation of motion 
for the low frequency component (p{x) is the following: 

(a^ + m^)4){x) + -[4>\x) + 30(x)(^2(x) + 3^'^{x)ip{x) + ip%x)] = 0. (3.23) 
6 

The effective equation of motion arises upon averaging over the (quantum) fluctuations of the 
high frequency field (p{x). Dropping the three-point function one obtains: 

((9^ + m^U(x) + ^ \^^(x) + 34)(x)(^'^(x))] = 0. (3.24) 
6 L J 

The last term on the left hand side represents the source induced by the action of the high 
frequency modes. It is a functional of (p{x). The action of the effective field theory may be 
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reconstructed from this equation. This approach is closely related to the Thermal Renormali- 
sation Group equation of D'Attanasio and Pietroni . 

In the broken phase the non-zero average value spontaneously generated below Tc is sepa- 
rated from the low-frequency part, (f){x) = (f) + (pix). The expectation value (p is determined by 
the effective equation 

m' + ^{v\x))^'^ + ^^' = 0. (3.25) 
2 

Our present goal is to determine the effective linear dynamics of the 0-field, therefore it is 
sufficient to study the linearised effective equation for (f){x): 

id' + ^^'mx) = -^^{^'ix))('\ (3.26) 

In this equation cj) is the solution of ( p. 25] ), for this reason the mass term is m^+^(^^+((y9^(x))*^°^ = 
=: M^. The linear response theory of is contained in the induced current, determined by 
(y,2(x))«. 

For the computation of the leading effect of the high-frequency modes in the low frequency 
projection of the equation of motion it is sufficient to study the two-point function A^(x, y) : = 
{ip{x)^{y)). 

As in Eqs. ( p.26|) , (|2.27|) we can derive 

{dl + M\x)) A>(a;, y) = 0, (9j + M\y)) A>(x, y) = 0, (3.27) 

where M^(x) = rn? + A0^(x)/2. 

After adding and subtracting the Wigner transform Q of the two equations appearing in 
(|3.27|) one arrives at the following equations for A(X, p) 

+ M\X))a{X,p) = 0, (3.28) 
/ d IdM^(X) 9 \ , , , , , 

('"•a^^2^aJ^)^(^-''>-°- f^-^^) 

The quantity appearing in the induced source is related to the Wigner transform of the two- 
point function through the relation 

{V\x))('^ = I ^A«(x = X,p). (3.30) 

The important limitation on the range of validity of the effective dynamics is expressed by 
the assumption that the second derivative with respect to X is negligible relative to and 
on the left hand side of Eq. (|3.28|). Then this equation is transformed simply into a local 



^For the introduction of the Wigner transform see Subsection |2.1.1| . 



36 



3. 1 . CLASSICAL KINETIC THEORY OF LANDA U DAMPING 



mass-shell condition, while Eq. ( |3.29| ) can be interpreted as a Boltzmann-equation for the 
phase-space "distribution function" A{X,p). Comparing Eq. ( p. 291 ) to Eq. ( |3.10| ) suggests the 
relation 

1 

-( 
2 

which supports the form of the particle Lagrangian 



M{X)F, = -d,M\X), (3.31) 



Lparticle = -M[0(x)], (3.32) 

defined in Subsection 

The background-independent solution A'^°^(p) is given by the well-known free correlator, 
slightly modified to account for the lower frequency cut appearing in the Fourier series expansion 
of <l>(x): 

AW(p) = (e(po) + n(bo|))2vr5(p2 _ M^), n(po) = ^^J_^ e(bo| - A). (3.33) 

We notice here, that in view of Eq. ( p.60|) or equivalently because is the solution of Eq. 
( p.25| ) (or Eq. (|3.19|) ) we have to use in place of the tree level mass square ?^e// one-loop 
resumed (Hartree) mass square = 10^. 

An iterated solution of Eq.( p.29D starting from (|3.33|) yields 



A.'.(.Y.p)^-kp-ax)-'^^. (3.34) 

By taking the Fourier-transform of the induced source jind{x) = — ■^0(v5^(x))^^'', and using the 
explicit form for A^^\p) and the relation d/dp^ = 2p^d/dp'^ one easily finds 

^Uk) = j (^^2vr5(p - M )^-- + j (3.35) 

The first term is the non-local contribution. One easily recognises after performing the po 
integration with help of the 5 function that it is exactly the Fourier transform of the r.h.s. of 
Eq. 

The second term calculated explicitly in the Appendix ^ is a local contribution, and it 
accounts for the renormalization of the self coupling A. 

3.1.4 The induced source 

The imaginary part of Eq. ( p.2(J| ) or ( p.35| ) determines the rate of the Landau-damping. Simple 
integration steps lead to (see Eq. [A. 121 ) 

lmjiko,k) = lmJ:m = -^j^mmM''kl)J^ dt-, 

= ^^mLi Q(^^-^)' io = l/Vl-(V|k|)2. (3.36) 
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This integral is zero if A > Mto, but for A < Mto it gives 



Imj(A;o,k) 



IGvr 



m- 



1 



,/3M/^l-(fco/|k|)2 



^e(|k|2-A;2)0fM 



A 



A' 2 



(3.37) 



independent of the value of A. 

The result has a very transparent interpretation. In the HTL-limit, when only the modes 
with much higher frequencies than any mass scale in the theory are taken into account, no 
Landau-damping arises. The effective theory is local! This also means that in a theory with 
only one massless degree of freedom, i.e. M = 0, there is no Landau damping. We will see, 
however, in Chapter [4.1| that in the 0(N) theory the massless Goldstone modes "feel" the 
presence of Landau damping due to the massive Higgs modes they interact with. In the case 
of 0(N) model the modes can transform into each other, this is the reason why we give up the 
attempt to describe the model with a kinetic theory. We will use instead a generalisation of 
the way presented in Subsection |3.L3| . 

Going beyond HTL, {ko ^ A ^ M) the 1-loop exact self-energy contribution, and so the 



1-loop exact Landau damping rate is calculated in the next Chapter in Subsection [4.1.3 
Using the formula 



n{s)ds = T ln(l - e"^*) - T ln(l - e"^'') 
one can rewrite the first term from the first line of Eq. (f4.32|) in the form 



(3.38) 



ImS 




4M2 



ki - |kP 



(3.39) 



This coincides with the result previously derived by Boyanovsky et al. in Ref. [^. The result 
of Eq. ( |3.37| ) holds for small values of and k as shown in Subsection [4.1. 6| . 



3.1.5 Discussion 

We have presented two equivalent methods of treating the effective dynamics of the low- 
frequency fluctuations of self-interacting scalar fields in the broken phase of the theory. The 
dynamics is proved to be non-local if the effect of fluctuation modes below the mass scale 
M is also taken into account. A fully local representation was proposed by superimposing a 
relativistic gas with specially chosen field dependent mass on the original field theory. 

The range of validity of the fully local version of the effective model can be ascertained 
only from its comparison with the result of lowering the separation scale A in the detailed 
integration over the fluctuations with different frequencies. From the comparison one learns 
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that the combined kinetic plus field theory is equivalent to the effective theory for the modes 
with ko <^ A < M, that is its validity goes beyond the Hard Thermal Loop approximation. 

There is yet another way of constructing a local theory, namely by introducing an auxiliary 
field W{x,\) defined through the equation 

V ■ kW{k,^^) = kocpik), (3.40) 

where v = (l,p/po)- Then by looking at Eq. (|3.13| ) one sees that fi{k,p) the deviation from 



equilibrium of the one particle distribution function is given by the equation 

Mk.p)^-lw(k..)i!^. (3.41) 
The non-local part of A^^\k,p) from Eq. ( p.34| )is given by the equation 

A^'\k,p) = -\^W{k, V) 2^5{p' - Mo^)^^. (3.42) 



We will see in Section [4.1| that the auxiliary field can be introduced even in a more general 
term, when we don't use the gradient expansion and a Boltzmannian kinetic evolution for these 
fields doesn't exist. 

We have also shown that in the broken phase of scalar theories and in the low-fco region the 
Landau-type damping phenomenon occurs according to an effective classical kinetic theory. In 
this context, the clue to its existence is the presence of the ~ 00 term in the fluctuating part 
of the local mass, which is responsible for the emergence of a linear source-amplitude relation 
with non-zero imaginary part. 

The extension of our discussion to the iV-component scalar flelds in the broken phase will 



be presented in Section 4.1. 
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3.2 Effective order parameter dynamics and the decay 
of a metastable vacuum state far from equilibrium 

The reaching of equihbrium from a metastable state involves a large number of interesting 
effects from instabilities observed in the mixed phase of first order phase transitions to the 



inflation in the early Universe [Q. The final state is reached in an irreversible process. 



The decay of metastable states is usually discussed in the framework of the nucleation 
scenario ||6^. It has been implemented in the form of saddle point expansions in classical 
p5| , and quantum systems [Q. This large amplitude instability is the first of the possible 
instabilities, suggested by mean field analysis. It consists of the creation of a bubble of the true 
vacuum larger than the critical size, embedded into the false one. 

Another possibility, the instability against fluctuations with infinitesimal amplitude leads 
to the spinodal phase separation. A recent observation made it clear that soft fluctuations of 
these inhomogeneous unstable modes generate in equilibrium the Maxwell construction by their 
tree-level contribution to the renormalization group flow The fluctuation induced flatness 



of the effective potential suggests the dominance of spinodal phase separation in equilibrium. 

Since the type of instability one observes, might depend essentially on the time scale of 
the observation, a detailed investigation of the time evolution can separate the effects of the 
two kinds of instabilities. This is made possible by large scale computer simulations of the 
thermalisation process in closed systems. 

Whether the relevant mechanism for a first order phase transition is the formation of bub- 
bles of the new phase, as described by thermal nucleation theory, or the gradual change of a 
large region of the sample, due to small amplitude spinodal instabilities described by spinodal 
decomposition is also an intriguing question in heavy ion physics where the actual expansion 
rate of the plasma may favour one or the other scenario [^] . 

Another question, left open by the Maxwell construction in equilibrium [^], concerns the 



structure of the vacuum with spontaneous symmetry breaking. In fact, the effective potential 
is related to the probability distribution of the order parameter and the Maxwell-cut applied 
to the former suggests that the latter is also degenerate in the mixed phase. Either we accept 
that the vacua with spontaneously broken symmetry correspond overwhelmingly to the mixed 
phase or a dynamical mechanism is seeked to eliminate the mixed phase from among the final 
states of the time evolution. 

In cosmology different slow-roll scenarios of inflation are being considered. Recent studies 
of the dynamics of inflaton fields with large number of components (large N limit) displayed 



for a first time a dynamical version of the Maxwell construction |70i . The Hartree type solution 



of the quantum dynamical equations leads to the conclusion that the order parameter might 
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get rest with finite probability at any value smaller than the position of the stable minimum of 
the tree level effective potential, corresponding to the stabilisation of a mixed state. 

Detailed investigations of the thermalisation phenomenon were performed also in noisy- 
damped systems, coupled to an external heat bath The relaxation to thermal 
equilibrium of the space averaged scalar field (the order parameter) starting from metastable 
initial values has been thoroughly investigated. In these simulations the damping coefficient is 
treated as an external control parameter. Using the numerical solution of the corresponding 
Langevin-equations the validity range of the analytical results for the homogenous nucleation 
mechanism has been explored. 

We focus on an alternative description of the decay process of the metastable vacuum state. 
The process is described exclusively in terms of the homogenous order parameter (OP) mode. 
The evolution of the OP is studied in interaction with the rest of the system as described 
by the reversible dynamical equations of motion of the full system. Careful analysis of its 
dynamics allows us to explore the effects of both kinds of the above mentioned instabilities. 
The transition of the order parameter from the metastable to the stable vacuum is induced by 
a homogenous external "magnetic" field, whose strength is systematically reduced. No random 
noise is introduced to represent any external heat bath, the friction coefficient of the effective 
order parameter dynamics is determined internally. 

Our model, a spacelike lattice regulated classical scalar $^ field theory in its broken sym- 



metry phase is introduced in Subsection |3.2.1| . In Subsection p.2.2| we describe the evolution of 
the system starting from order parameter values near a metastable point which relaxes first to 
this state under the combined effect of parametric resonances and spinodal instabilities. The 
second stage of the transition to the stable ground state is the actual focus of our discussion. 

Characteristic intervals of the observed order parameter trajectory are reinterpreted as being 
the solutions of some effective point-particle equation of motion, which displays dissipation 
effects explicitly. Our approach can be understood also as the real-time version of the lowest 
mode approximation used for the estimation of finite size dependences in Euclidean field theory 



75, 76 1 . In this sense our approach can be considered also as the numerical implementation of 



a real time renormalisation group strategy. 

One of our principal goals is to reconstruct the thermodynamics of the classical "OP- 
ensemble" on the (meta)stable branches of the OP-trajectory (Subsection |3.2.3|) . Its dissipative 
dynamical equations near equilibria will be established. On the transition trajectory we shall 
elaborate on the presence of the Maxwell construction in the effective dynamics describing the 
motion after nucleation (Subsection p.2.4| ). The statistical aspects of the approach to the equi- 
librium are established for reference and comparison in Subsection p.2.5| . The conclusions of 



this investigation are summarised in Subsection 3.2.6 
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3.2.1 Classical cut-ofF $^-theory on lattice 



The energy functional of a classical system in a two-dimensional box of size Ld coupled to an 
external magnetic field of strength hd is of the following form: 



d<^d 
did 



1 



1 



24 



(3.43) 



The index d is introduced to distinguish the dimensionfuU quantities from the dimensionless 
ones, defined by the relations (for rn? < 0): 



t = td\m\ 



X = Xd\m\ 
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6 m 



(3.44) 



For the spatial discretisation one introduces a lattice of size Ld = Nad = ^(^j^^ 
The energy functional of the lattice system can be written as 
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6\m\ 
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(3.45) 



(Here we have introduced the dimensionless time-step at, which should be chosen much smaller 
than a, and n denotes the lattice site vectors.) The equation of motion to be solved numerically 
is the following: 

+ at) + $n(t - at) - 2$„(t) - 4 E(*n+i(^) + ^^n-iit) ' 2$n(t)) 



+a2(-$„ + $3 -/,) = 0. 

The initial conditions for Eq. ( |3.46D were chosen as 

$n(t = 0) = 0, ^^(t = 0) = % + ^^i. 



(3.46) 



(3.47) 



The random variable ^ is distributed evenly on the interval (—1/2, 1/2). Therefore the starting 
OP-value is $o- The energy density E/Na^ is controlled through the magnitude of $i. In this 
study we have chosen $o = 0.815 and $i = A/\/Q. The latter corresponds to a temperature 
value Tj = 0.57 in the metastable equilibrium (from Eq.( 3.52D ). This value is much below the 
critical temperature of the system (Tc ~ 1.5Tj). It has been checked that at this energy density 
all other choices of $o > 0, for fixed /i, find a unique metastable equilibrium. 

Eq.( 3.46 ) was solved with a = 1 and with typical at values in the range (0.01-0.09). It has 
been checked that the statistical characteristics of the time evolution is not sensitive to the 
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variation of at, though the "release" time (the moment of the transition from metastabihty to 
the true ground state) in any single run with given initial conditions might change considerably 
under the variation of at /a. Three lattice sizes were systematically explored: = 64, 128, 256. 
Several single runs were performed also for = 512 and = 1000 with the aim to analyze in 
more detail some self-averaging physical quantities on different portions of the trajectory under 
the assumption of the ergodicity of the system. The magnetic field h inducing the transition 
was varied in the range h G —(0, 0.08)/v/6- For the reconstruction of the effective potential felt 
by the OP also positive values were chosen up to /i = 0.5/^/6. The smaller the value of was 
fixed, the longer the "release" times have grown on the average. For this reason also the runs 
were prolongated with decreasing \h\, and for the smallest \h\ the length of a run reached up 
to (10^ — 10'')|m|~^ until the transition took place. 

For a careful comparison of our transition statistics with the generally used statistical nu- 
cleation theory, and also for understanding the systematics of its change when h has been 
diminished a large number of [h, N) pairs were used in this analysis. Altogether 24 422 tran- 
sition events have been recorded (for A^ = 64: 16 908, A^ = 128: 2 903, A^ = 256: 4 411). For 
the largest systems at the smallest h the event rate was 1-2/day/ 400 MHz-processor. 



3.2.2 Time-history of the order parameter 



A typical OP-history is displayed in Fig. |3.1| . In the same figure we show also the history of 
the OP mean square (MS)-fluctuation (($2) _ ($)2) and of its third moment ((($ - ($))^)). 
The evolution of the non-zero k modes is demonstrated in Fig. 3^, where the averaged kinetic 
energy content of the |k| < 2.5 and |k| > 2.5 regions is followed. Although the separation value 
is somewhat arbitrary, namely it divides into two nearly equal groups the spatial frequencies 
available in the lattice system, the figure demonstrates the most important features of the 
evolution of the power in the low-|k| and high-|k| modes. 

In general, five qualitatively distinct parts of the trajectory can be distinguished, although 
some of the first three might be missing for some initial configurations and/or magnetic field 
strengths. 

The OP-motion usually starts with large amplitude damped oscillations. The "white noise" 
initial condition of Eq.( p.47|) corresponds to a k-independent Fourier amphtude distribution, 
therefore the initial distribution of the kinetic energy is ~ cj^dkl). During this period, in the 
power spectrum of the kinetic energy, first a single sharp peak shows up at a resonating |k|-value 
(|k| ~ 1.5), which breaks up into several peaks (|k| < 1) at later times due to the non-linear 
interaction of the modes, see Fig. (|3.2| ). At the end of the first period the whole |k| < 1 range 
gets increased power, the |k| > 1.5 part of the power spectrum does not seem to change. 

Next a slow, almost linear (modulated) decrease of the OP follows. At the same time its 
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Figure 3.1: The time evolution of the order parameter, its MS fluctuation and the third moment. 
The example is selected from runs on a = 512 lattice with h = — 0.04/-\/6 external source 
strength. 
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Figure 3.2: Early time power spectrum. 
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Figure 3.3: The time evolution of the kinetic energy content of the |k| > 2.5 and |k| < 2.5 
regions averaged over the corresponding |k|-intervals. The example is the same as in Fig. 1. 
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MS fluctuation increases linearly. In Fig. p.3| an energy flow towards the low-|k| part can be 
observed, which proceeds through the excitation of single modes in this part. As a result 
the minimum of the effective potential is continously shifted to smaller <l>-values, as if the 
temperature would gradually increase (see the left part of Figj3.1|). The full (microscopic) 
kinetic energy density shows in this period less than 5% variation. In view of the picture based 
on the Maxwell construction it is rather surprising that independently of the initial conditions 
the OP converges towards a well-defined absolute value, depending only on the total energy 
density. 

On the third portion the average value of the OP and its moments stay constant (see the 
right picture in Fig. p.l|) . The average energy content of the low and high-|k| part of the spectra 
is nearly the same. This suggests the establishment of a sort of thermal (meta)equilibrium. 

In terms of the terminology introduced for the inflation, the first period leading to this 
(quasi) -stationary state can be called preheating, and the second reheating. Directly before 
the moment of the transition to the true vacuum a peak appears in the power spectrum in 
the narrow neighbourhood of k = with varying position in time. On the snapshots of the 
real space configurations a set of randomly distributed small bubbles of the true ground state 
appear with a radius increasing in time, until one of the bubbles exceeds the critical size. 

The fourth portion of the motion is the transition itself. The value of the OP MS-fluctuation 
increases by about a factor of three and the temporal width of this transient increase measures 
very well the transition time. The transition time decreases on larger lattices, the height of the 
jump in the OP-fluctuation is not sensitive to the lattice size. Also the third reduced moment 
shows a characteristic variation. An increase of the temperature proceeds smoothly during the 
transition of the OP to its stable value (Fig. |3.3| ). The slight separation of the two curves in 



Fig. |3.3| gives a feeling on the degree of uniformity of the temperature variation of the different 
modes. 

The last portion of the trajectory represents stable (thermalised) oscillations around the 
true ground state. Here a complete equilibration of the power spectrum can be observed (see 



Fig. 3.3) corresponding to a somewhat increased temperature (T ~ 0.6). 



3.2.3 Motion near the (meta)stable point 

Our analysis of the motion around the (meta) stable value of the order parameter explores the 
consequences of assuming the ergodicity hypothesis for a sufficiently long finite time interval, 
after the system has already reached the equilibrium. The equilibrium is characterised by a 
limiting probability density in the configuration space. The averaging with this density should 
provide the same value as the one yielded by a single long time evolution when subsequent 
configurations are used in constructing the statistics. Our "statistical system" is now a single 
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degree of freedom, the order parameter of the lattice system, interacting with all other {k ^ 0) 
modes. 

The effective equation of motion can be thought to result from the application of a "molec- 
ular dynamical renormalization group" , to our microscopical equations. The blocking in space 
is performed by projecting the field configuration $(x, t) on the OP It represents the 

infrared (IR) end point of such a blocking whose effective theory is now reconstructed from the 
actual time dependence found numerically. Combining ergodicity of the full system with the 
renormalization group (RG) concept we arrive to the conclusion that ensemble averages of any 
OP-function coincide with time averages of the same function. 

In view of the RG concept we look for an effective equation of motion where the value of 
OP is determined exclusively by its values preceding in time. Assuming the absence of long 
memory effects a "gradient" expansion in time can be envisaged, leading to a local differential 
equation. We introduce into this phenomenological "Newton-type" equation a term violating 
time-reversal invariance. We are not able to derive it from the original system, we only wish 
to test its presence. The sustained motion of OP, however, requires, in this case the presence 
of a random "force" term, too. The deterministic part of the "force" is expected to be related 
to the equilibrium effective potential, since this object determines the stationary probability 
distribution for the OP near equilibrium. 

The above considerations lead us to write down the equation of motion, which represents a 
linear relation between the acceleration and the velocity of the order parameter: 

$d + Vdi^A -hd + ^^^^^ = Cd, (3.48) 

where (d is a noise term. In the corresponding dimensionless equation of motion the following 
new rescaled quantities will appear: 

Vd = Imlr], Cd = C\mf^. (3.49) 

The fitting procedure for the coefficients on the left hand side of Eq. ( p.48| ) was the following. 
For a given interval of time It the region of the order parameter space visited by the system 
{$(t)|t G If} was divided into small bins. Having defined the time set 

T$j^ = {t G /t|$fe < $(t) < + A$} corresponding to a given bin, the linear relation ^(tf,) = 
—rj{^b)^{tb) — /(${,) was fitted with the method of least squares using the ^(tft), $(^6) data 
measured at time moments tf, belonging to T<j,j_. We have obtained in this way the coefficient 
functions ?7($fe), f{^b)- Once the functions //($) and /($) = —h + V^g($) are determined, we 
evaluate for each time t the expression $(t) + ?7($(t))$(t) + /($(t)). Its actual value determines 
the random noise function ((t), whose statistical features (autocorrelation) should be extracted 
from the data. 
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The effective force /($) calculated from the time-average of the oscillatory motion around 
the equilibrium, is expected to agree with the force coming from the theoretically determined 
finite temperature effective potential calculated perturbatively in the cut-off two-dimensional 



field theory for some appropriately chosen value of the temperature With one-loop accu- 
racy the expected equahty reads: 



f (^^) measured h -\- 



F($) + ^F"(<I>) fl + log- 



+ 0(T2). (3.50) 

theory 



The expressions on the right hand side of this equation are connected to the dimensionfuU 
quantities of the original one-loop computation in the following way: 

^ ,T,, V{<^) = -U' + U' = -^V,, A = -. (3.51) 



6|mP 2 4 6|mp a 

The dimensionless temperature is defined by the time-average of the kinetic energy based on 
the assumption that in the effective theory of the OP it has the usual expression in terms of 

hmi fdt'lm\x)r^l¥' = '^. (3.52) 

One should note that this definition in the dimensionfuU version implies a "Boltzmann-cons- 
tant", |mp multiplying T^. 

By measuring the order parameter average for different values of the external field h we ob- 
tain the magnetisation curve of the system. In the numerical work we restarted the computation 
at different values of the external magnetic field but one might as well change h adiabatically 
and measure the force law at each (quasi)equilibrium point. The results should agree in the 
stable regime, up to the phenomenon of the hysteresis. The resulting curve can be viewed as 
the numerical Legendre transformation by identifying the external source h with the derivative 
of the effective potential, 

h,measured 

) = 0. (3.53) 

It was found that f{^)measured IS always vanishing at $ = {^) h,measured, thus the force acting 
on the order parameter at the equilibrium position {^)h induced by the external source h is 
indeed always —V^jrj{{^)h)- 

The two sides of the relation ( |3.50 ) are shown in Fig. |3.4| . To the values of the external 



field used in preparing Fig.|3 = -0.5, -0.02, 0, 0.5, 1, 2) single runs were selected by the 
requirement that the system stayed in the metastable vacuum up to the time 10^. The force 
shown in Fig.|3.4| demonstrates that not only the equilibrium positions but also the fluctuations 



of the OP are governed by the static effective potential. Similar measurements, performed on 
100 X 100 lattices showed no finite size effects in the stable regime, h > 0. The error of /($) 
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Figure 3.4: The force /($) as the function of the order parameter on 40 x 40 lattices with 
T = 1/30 and 1/3. The force is measured by shifting the center of motion to different $ values 
by applying appropriate h fields to the system. The finite size effects are negligible, see text. 
For comparison we display the force arising from the derivative of the one-loop and for T = 1/3 
also of the two-loop expressions of the effective potential. 



is not shown in the Figure, since the typical values (~ 0.004 for T = 1/3 and ~ 0.002 for 
T = 1/30) are too small to be displayed. 



Another piece in the effective equation of motion ( p.48|) is the friction term whose presence 



indicates the dynamical breakdown of the time reversal symmetry. The friction coefficient 
ri{^) proves clearly no n- vanishing and shows only weak dependence on the actual value of 
the OP around its equilibrium position. The breakdown of the time-reflection symmetry in a 
closed system must arise only in presence of infinitely many degrees of freedom. Till then only 
statements on the Poincare-time can be made. Thus our non-zero results for r] require further 
clarifications. 

The point is that there are two types of infinities, controlled by the temporal UV and the IR 
cutoffs, respectively. The spontaneous symmetry breaking is driven by the IR modes, and the 
non-trivial minima of the potential energy arise from the presence of infinitely many degrees 
of freedom in the IR (thermodynamical limit). On the contrary, the dynamical symmetry 



breaking is the result of the effects of the derivative terms in the action and the infinitely 
many UV modes (continuum limit). The breakdown of the time reversal, being related to a 
linear time-derivative term in the effective equation of motion, should come from the UV, the 
short time behaviour of the system. 

In fact, one expects no friction when the UV cutoff, at is so small, that not enough energy 
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Figure 3.5: The friction ri{{^)h=o) as the function of the timehke lattice spacing bt = nat- The 
time scale of the fastest mode is r ^ 2 and 0.22 for a = 1 and 0.1, respectively, according to 
the perturbation expansion. Typical error bars appear on the curves both for the IR and UV 
regimes. 

can be dissipated during such a short time. In quantitative terms one should have 

at>r = —^^ + 0{T), (3.54) 



where r is the time scale of the fastest mode in the system. The right hand side relation of 
Eg. ( p. 541) gives an estimate of the maximal frequency from the free dispersion relation Pq = 
4(sin^ PxO./^ + sin^ pya/2) /a^ + M'^(T) of the lattice Hamiltonian system for fixed a in the limit 
at — > 0. The fast modes absorb energy from the OP in a single time step for at > r, and the 
friction term should appear in the effective equation of motion. 
With help of the temporal blocking 



<l>(t)^$„(t) = - y ^t + kat) (3.55) 

k=-n/2 

one can construct the trajectories $n(^) corresponding to larger values of the time cut-off, 
ht = nat- Such blocking was performed in time up to n = 2000 and a discrete time equation of 



motion of the form Eg. ( 3.48 ) was reconstructed for the blocked trajectory $„(t). The non-trivial 



dependence of r] on the temporal cutoff at is shown in Fig. pT5 . 

One can distinguish two regimes separated by a crossover apparently independent of the 
lattice size, located at bt ~ 0.2r. A scaling behavior is observed on the UV side, where ri{bt) 
tends to zero with a critical exponent close to 1 and the dissipative force being constant. This 
result should be independent of the actual blocking details. In the other regime, on the IR side 
T] goes over first into a ferindependent regime, corresponding to the saturation of the energy 
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transfer from the OP and giving a stable, microscopic definition of the friction coefficient. 
Finally, in the far IR part a qualitatively different oscillatory behaviour sets in. 

The location of the crossover from the UV scaling regime to the plateau can be understood 
by writing the fluctuation dissipation theorem in the corresponding discretised form: {('^) = 
2riT/bt. The linearly increasing regime of T]{bt) implies constant second moment for the noise. 
When ri{bt) reaches the plateau the second moment of the noise decreases like 1/bt. The 
crossover therefore is located at the autocorrelation time scale of the noise. 

A qualitative interpretation of the oscillatory IR regime can be based on the observed 
small amplitude beating phenomenon in the OP trajectory. This can be recognized by closer 



inspection of the left side of Fig. |3.1| , which persists further also on the right side. It is reflected 
in the OP autocorrelation function, too, since in the course of the blocking an interference effect 
occurs on the right hand side of Eq.( |3.55| ) due to this regularity. This feature is relevant to the 
value of f]{bt), responsible for the decay of all kinds of fluctuations. 

The appareance of peaks in r]{bt) at both the maximal destructive and constructive inter- 
ferences can be modeled semi-quantitatively by identifying the beating part of the OP-motion 
Re(5<I> = Re(5$o exp(za;t) with the stationary solution of a single weakly damped driven har- 
monic oscillator, 5$ + 7^5$ + uJq6^ = foexp{iujt). The blocking acts on the trajectory 6^{t) 
as 5^{t) u6^{t), where u = {exp{iujbt) — l)/iujbt. It does not change the relative phase of 
the driving force and of the blocked oscillation amplitude, leading to a relation between the 
parameters of the original and the blocked equation of motion: 

— uj + irfuj + uJq = —uj u +ifiuJu + ujQ. (3.56) 

The friction coefficient for the blocked trajectory turns out to be = [r] + uj\m.{v?)) /'Re{u) . 
It is easy to see that Re(M) is vanishing at maximal constructive and destructive interferences. 
This provides singularities in f]. Whenever Re(u) = we have Im(u^) = and numerator 
changes sign in the vicinity of the singularity. Thus fj > Q apart for a short time interval 
around the singularities where the non-harmonic features should stabilise the fluctuations and 
keep > 0, as observed in our simulation. 

The a dependence appearing in Fig.^]5| arises from the following two effects. One is that 
for larger a the maximal oscillation frequency is smaller and the time resolution of the system 
becomes cruder. Another is that larger a represent bigger physical volume, many more soft 
modes and less harmonic system, which tends to invalidate the simple picture based on a single 
harmonic mode. As a result the effects of oscillatory nature will be smeared, as one clearly 
recognizes in the flgure. 

The quality of any proposed deterministic equation of motion (e.g. equations similar to 
Eq.( p.48| ) with zero on the right hand side), can be judged by the amplitude and the auto- 
correlation of its error term, the noise term (. The amplitude of the noise was found at least 
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two-three order of magnitude below the average level of the force as fitted to Eq. ( |3.48| ). The 
autocorrelation function of the noise of Eq. ( p.48| ) appeared to be local, approximately of 
the form {Cd{t)Cd{t')) ~ ^"{t — t')- The status of the fluctuation-dissipation theorem will be 
investigated for more complex field theoretical systems in future investigations. 

The selfconsistency of the definition of the temperature in Eq. (p.52|) and also the establish- 
ment of the thermal equilibrium can be tested further by plotting histograms for the following 
quantities: 



Eh 



Ep 
E, 



1 /ci$Y 

2 \ dt) 

2 ^ 

Ek + Ep, 



(3.57) 



where /x^ is the slope of the force as the function of the order parameter, determined numerically. 



Typical results for the energy-histograms are shown in Fig. |3.6| . It shows perfect agreement of all 
slopes and good agreement with the expectations based on equipartition of the energy between 
the kinetic and the potential parts. We have checked for several temperatures, that the tem- 
perature determined by these histograms and from the |k|-spectra of the kinetic energy agree. 
This piece of information is parallel to the recent detailed investigations of the thermalisation 
in (1 + l)-dimensional $'^-theories ||79| , |80 . 



1000 




Figure 3.6: The histogram of the energies Et, Ek and Ep of Eq.( (|3.57| )) on a 256 x 256 lattice 
with h = 2/VQ- 



It is worthwhile noting that the relaxation from a given initial condition to the thermally 
distributed state like in this figure takes at least one order of magnitude longer time for systems 
in the symmetric phase compared to the spontaneously broken case. At higher energy density 
(temperature) one would expect larger collision frequency, therefore shorter thermalisation time. 
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The opposite result hints to the importance of slow, soft modes, whose presence is due to the 
symmetry breaking mechanism. We shall argue in the next Subsection that these modes are 
responsible for the realization of the Maxwell-cut in the potential term of the effective equation 
of motion for the OP. What we find remarkable is that these modes are present not only in the 
mixed phase but also near (meta) stable equilibria, among the dynamical fluctuations around the 
ordered vacuum in contrast to the massive perturbative excitation spectrum in the equilibrium. 

We complete the analysis of this subsection with the remark, that after the transition to the 
stable vacuum shortly a thermalised distribution is recovered for the OP at somewhat increased 
temperature which agrees with the value given by the equipartition. 



3.2.4 Jump from the metastable to the stable vacuum 

As it has been emphasized in the Introduction, there are (at least) two different descriptions of 
the transition from the metastable state to the stable one. In the first approach an expansion 
around the lower pass, the critical bubble (bounce) configuration yields a detailed space-time 
picture and the transition rate by means of the analytic continuation of the potential experi- 
enced by the OP near the (meta)stable position |66| . 

Another possibility is to provide for the OP a probabilistic description, obtained by the 
elimination of all other degrees of freedom in some kind of blocking procedure. This description 
is based usually on some underlying Master equation for the probability distribution of the OP 
and the resulting Fokker-Planck equation [^. The transition to the stable state appears in 
this approach formally as a tunneling solution of the Fokker-Planck equation. The probabilistic 
feature of the dynamics of the OP is supposed to arise from near the (meta) stable equilibrium 
assuming a statistical ensemble of the initial conditions. 

In our simulation we find results analogous with the predictions of the probabilistic descrip- 
tion by analysing the OP-motion starting from a single, well defined initial condition. One 
might wonder at this point if it is possible to understand the probabilistic tunneling of the 
OP by following the system from a unique initial condition. The self averaging in time can 
not be used for this argument since such a transition occurs only once during the evolution in 
a (quasi) irreversible manner. We have to develop a third approach to the metastable— ^^stable 
transition. 

In general, the effective equation of motion for $(t) reflects the typical landscape of the 
microscopic potential energy functional around the actual point $(x, t) in the configuration 
space. The classical origin of what appears as a tunneling on the Fokker-Planck level must 
be the arrival of <l>(x, t) to the vicinity of some narrow valley opening up towards the stable 
vacuum. In traversing this valley the landscape changes and the typical fluctuations will be 
different from those felt in the metastable regime. The constants parametrising the effective 
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Figure 3.7: The force /($) as the function of the order parameter on 100 x 100 lattice with 
h = -0.04/^6. Each curve corresponds to time intervals of length 500|m| ^. The interval 
corresponding to the continous curve ends 1000|m|~^ before the transition, the long dashed 
curve refers to the next time interval and the short dashed curve is the force measured during 
the last time interval preceding the transition. One clearly observes the bending down left from 
the central OP-value in the metastable regime. Typical error bars are shown on the leftmost 
points of the three curves. 
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equation of motion must reflect this change. 

Our goal in this Subsection is to construct an effective description of the transition to the 
stable state by carefully tracing the time evolution of the OP. This will be achieved by project- 
ing the microscopic equation of motion onto the homogeneous mode and phenomenologically 
parametrising it similarly to Eq.( |3.48D . 

As long as the system is far from the narrow valley of the instability the force is time inde- 
pendent and agrees with the force derived from the perturbative effective potential according 
to the part of Fig. |3.4| corresponding to $ > 0.9. When the system arrives close to the entrance 
of the unstable valley, ($ ~ 0.8), the soft modes start to be important. This is reflected in 
the slight glitch in the leftmost piece of the measured force law in Fig. 3^. A sequence of 



glitches results in a situation depicted in Fig. |3.7| . It shows the force acting on the OP in three 
successive time intervals preceding the event of tunneling for h = —0.0A/\/Q. The /($) curves 
determined using Eq. (|3.48| ) in the disjoint time subintervals coincide within error bars for all, 
but the last one. In this last interval preceding directly the transition towards the direction 
of negative $ values, the fitted force bends down and its average becomes (a small positive) 
constant. This is characteristic feature of the instants when the system finds the entrance into 
the unstable potential valley. The vanishing of the force is an indication for the dynamical 
realization of the Maxwell-cut. The OP moves fast through the valley and the method of fitting 
the trajectory to Eq. (|3.48|) for finding the force fails due to the insufficient statistics. 
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The fluctuation moments depicted in Fig. tell a bit more about this region. The 
increased values of the moments, the renormalized coupling constants in Wilsonian sense at 
vanishing momentum, indicate the enhanced importance of the soft interactions as the OP 
tunnels through the mean field potential barrier. This softening makes the OP fluctuating with 
larger amplitude. The valley of instability is in a surprising manner flatter than the typical 
landscape around equilibrium. The flatness along the motion of the OP (the mode k = 
in momentum space) comes from the Maxwell-cut. The average curvature of the potential 
in the transverse directions, i.e. for modes with |k| ^ 0, can be estimated by the second 
functional derivative of the two dimensional effective action with cutoff |k|, which can be taken 
as k^ + V^'($). The increased MS fluctuation of the OP corresponds to a decrease in V^^q{^) 
in the valley. It pushes down V^{^) also for small non- vanishing |k|, in the low momentum 
regime which is expected to be the most influenced by the changing landscape. 

The lesson of Fig. is that the potential itself should be considered as a fluctuating 
quantity. This helps to translate into more quantitative terms the above qualitative points. We 
discuss the projection of the microscopic equation of motion ( 3.46| ) onto the zero momentum 
sector, 

= $- $ + $3 + 3^^$ + - h = ^ + V(^^^{^), (3.58) 

where the symbol (/^"^ means the space average of {Tp^ = 0, $(t, x) = $(t) + (p{-K,t)). The 
instant potential introduced in the second line contains a deterministic piece, which is the sum 
of the tree-level potential and of the slowly varying part the second and third moments (eg. 
(v'"^)cZei)- This last feature clearly appears graphycally in Fig.^TT] and a simple model based 
on a two-phase model of the transition period will be constructed below to account for it. The 
remaining oscillating pieces of the moments provide the probabilistic fluctuating contribution 
to the instant potential: 

^inst('^) = - + + (^''U^ + li^'^U^' + Co<f + ICi^'- (3.59) 
The additive (Co) and the multiplicative (Ci) noises are given by the differences 

Co(t) = ^''{t)-{^''u, 

Ci(t) = ^''{t)-{^''u. (3.60) 

Note that no friction terms can be introduced in a natural way into the system (|3.58| - |3l60|) . 
Therefore we face the intriguing question, how irreversibility is realised in such a system. The 
time-correlation matrix Ci{t)Cj{t + r) appears to us to be the key object for its investigation, 
to which we plan to return in the future. 

In order to gain more insight how this works we build into Eq. ( |3.58| ) the consequences of the 
mixed two-phase picture of the phase transformation. The microscopical basis for this picture 
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is provided by the thermal nucleation whose quantitative discussion is given in the Subsection 

More specifically, we assume that the space can be splitted into sharp domains (neglecting 
the thickness of the walls in between), where the field is the sum of the constant background 
values $o± and the fiuctuations if>± around it, 

$±(x,t) = <l>o± + ^±(x,t). (3.61) 

We assume local equilibrium in both phases, based on the smooth evolution of the temperature 
as displayed in Fig .[3 .31. 

The actual value of the order parameter is determined by the surface ratio p(t) occupied by 
the stable phase: 

$(t) = p{t)%. + (1 - p{t))%+. (3.62) 

Simple calculation then yields 

(^VW = ~ y (Wi^'' - + - $^(t) + (3.63) 



^0+ ~ ^0- ^ 



+ ^li^,tf - 3$(t)(^2^)rfet(t) - $'(t), (3.64) 

where the volume averages should be read off the corresponding equilibria on the two sides of 

the transition. If one takes the values of ^o±, from the respective equilibria determined in 

the same simulation, a quite accurate description of the shape of the two fiuctuation moments 

arises in the whole transition region and its close neighbourhood using $(t) to parametrize their 

V 

t-dependence. (Note that ^\{x.,^) does not depend on time.) The deterministic part of the 
moments are therefore well-defined functions of $. In this way, a simple explicit construction can 
be given for the effective noisy equation of the OP-motion, if the above pieces are supplemented 
by the correlation characteristics of the two kinds of noises. 

The RMS (root mean square) fluctuations found numerically on the transition part of the 
trajectory and its close neighbourhood are the following: 



Qit) = 0.0063(10), VC?W = 0.0054(10), foriV = 64 (3.65) 

(its magnitude increases with the size of the system). The magnitude of their equilibrium cross 
correlation was found ~ 10~^. 

As a corollary of this construction one can demonstrate the absence of the deterministic 
part of the acceleration of the order parameter in the transition period. 
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Substituting Eqs.( p.63D , ( |3.64| ) into Eq.( |3.58| ), one finds for tlie deterministic part of the 
force, 



' <l>o+-$o- / $o+-$o- 



The average of the equations of motion in the respective equihbria, 

(<l>i(x,t))-<l>o±-/i = (3.67) 



imphes the vanishing of the deterministic force in Eq. ( |3.66| ), when exploiting the equahty of 



the volume and the ensemble average in this case. Eq. ( |3.58| ) transforms into 

$(t)+Co(t)+3Ci(t)$(t) = 0. (3.68) 

This is the dynamical realization of the Maxwell-construction holding when the mixed phase 
model with local equilibrium is valid. 



If the force is calculated from the full instant potential, ^('^')inst (|3-59|) , it depends 

parametrically on the moments (f'^ (t) and (f^ (t). The approximate trajectory '^'inst(^)) de- 
fined by minimising V^('^')inst with respect to $, where the moments are taken from the nu- 
merically determined time evolution, reproduced accurately the observed OP, $(t) with the 
following RMS/unit time: 



$(t) - $insJt))2 = 0.0014(5), forAr = 64. (3.69) 

This construction interpretes the OP-trajectory as a continous deformation of the instant po- 
tential with the OP "sitting" permanently in the actual minimum. Notice that such motion is 
possible only if the OP continously undergoes some sort of dissipation. 

The vanishing of the acceleration was checked by comparing the computer generated OP 
trajectory in the transition region with a ballistic motion in a viscous medium. In particular 
it was tested that the ratio h/^ is nearly time and h independent for small enough h. The 
dynamical friction measured by the above ratio tends to a constant with decreasing h at fixed 
lattice size. Although the error for its value in each individual transition is rather small, the 
central values obtained in different runs fluctuate quite strongly, which leads eventually to a 
large error of the mean calculated as an average of runs with different initial configurations 
having the same energy density. 

It is worth to note that the order of magnitude of these "renormalized" , i.e. IR determined 
values of the friction coefficient agree with the peak value in Fig. ^]5| for = 64. The "running" 
friction coefficient determined by the blocking in time, however, drops as bt is increased and 
takes extremely small values at bt ~ 20, the average time length of the ballistic fits for the 
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transition. This serves an example for the dependence of the renormahzation group flow on the 
details of the blocking. It remains to be understood whether the agreement between the peak 
value and the ballistic fit is an accident or follows from the internal dynamics. 



3.2.5 Connection with the nucleation picture 



In this subsection we analyse the transition using the more conventional statistical approach of 
thermal nucleation theory. 

From the study (on lattices up to = 512) of the detailed microscopical field configurations 
it turned out that the phase transformation starts by the nucleation of a single bubble of the 
stable ($ < 0) phase. Late-coming further large bubbles are aggregated to it as well as the 
small size (consisting of ribubue < 50 — 60 joint sites) bubbles. Following the nucleation the 
expansion rate of the large bubble governs the rate of change of the order parameter. To a very 
good approximation each individual transition could have been characterised by a constant 
value of $ in this interval. Two mechanisms are known to lead to this behavior. The first is 
scattering of hard waves ("particles") off the bubble wall, while in the second the expansion 
velocity is limited by the diffusive aggregation of smaller bubbles 
to be the dominant in the kinetic Ising model . 



. The latter process seem 
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Figure 3.8: The normalised "release" time histogram counting the frequency of a certain time 
moment-bin when the order parameter first takes a negative value. The figure represents the 
statistics based on 750 events at h = — 0.08/-\/6 on an iV = 64 lattice. 



The statistics of the "release" time of the supercritical bubble from the metastable state 
shows at first sight rather peculiar characteristics. The binned histogram for larger values 
of h shows an asymmetric peaked structure which apparently deviates from the exponential 
distribution characterising the thermal nucleation scenario (Fig. ^.SI ). The very early transitions 
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(t < tmax) seem to be suppressed for small values of h. We expect they correspond to transitions, 
which happen before the system reaches the metastable state. 

The fall-off for t > tmax starts nearly exponentially, but the histogram develops a very 
long tail when h ^ (). The bigger sample is used for estimating the probability distribution, 
the more suppressed is the weight of these events in the normalised distribution. In practice 
this means a longer time interval where a good linear fit can be obtained to the log-linear 
histogram. Eventually, the separation of a clean exponential signal is possible, and the slope 
can be compared with predictions of the nucleation theory. 

Knowing the frequency histogram of the time t that elapsed before the system has escaped 
from the metastable region by nucleating a growing bubble of the stable state, we can construct 
the nucleation rate per unit area: T{h). Indeed, for a lattice size L we can read off F by fitting 
to the expected form Phit) ~ exp{—tV{h)L'^). 

The nucleation rate is proportional to the volume of the system L. In Fig.|3.9| we see evidence 



for the size independence of the transition rate per unit surface. The logarithm of this quantity 



can be estimated following the standard nucleation theory |6^, |66|, |8J] where the nucleation 
rate per unit area is given by: 



r ~ e 



-S2/T 



with ^2 = - S{<^f{h)), 



(3.70) 



$B is the bounce solution and S the Euclidean action S = J d^x \_\{d^(^Y -\- with 
defined below in Eq. (|3.74D. 
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Figure 3.9: Nucleation rate for unit lattice surface versus the reciprocal of the magnetic field 
for different lattice sizes. The error bars are smaller than the size of the symbols representing 
the data points. 



In the so called thin wall approximation, when the energy density difference beween the 
true and the false vacuums is tiny the coefficient ^2 can be calculated analiticaly and in the 
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relation 

\nT = K-^. (3.71) 

5*2 is determined approximately by the action of the bounce solution, connecting the two vacua. 
In the investigation below we have studied at fixed temperature the /i-dependence of S2/T, 
assuming the /i-independence of K, although this might be not the case. 

The bounce is a solution of the Euclidean Euler-Lagrange equation of a theory defined by 
the action S'($) = / c/V (|(V$) + |m2<|)2 + A<|.4 _^ g^^h that (see fo exemple |5§) 

1. $B approaches the false vacuum $/(/?-) at the infinity, 

2. $B is not a constant, 

3. S{^b) is the minimal amoung the solutions obeying 1. and 2. . 

Since the nucleation proceeds through approximately spherical bubbles we need an 0(2)- 
invariant solution of the equation of motion i.e. the solution of 

1^$ A 



+ -— = m'^ + -^' + h, (3.72) 
dr"^ r ar b 

with the boundary conditions ensuring that the field be regular at the origin ($'^(0) = 0) and 
in the false vacuum at spatial infinity ($5(00) = $j (/;,)). 

As it is well-known, in the thin wall approximation the differential equation ( p. 721 ) reduces 

where t/+($) is the a symmetric function of $, with minima at = ±y^|m|, obtained by 
replacing /i by in the tree level potential of Lagrangian and dropping the constant term: 



24 24 V A / 2 A 



Multiplying equation ( |3.73D by $' one can see that it admits a first integral [|($')^ — = 
whose value is determined by the condition that $(00) = and ^ 
finiteness of the energy. So, we have 



because of the 



2($')'-f/+($) = -f/+($/) = 0. (3.75) 

This equation determines a monotonic $ which goes from $t to $/. The solution depends on 
a single integration constant f, a point at which the $ takes the average of its two extreme 
values and $/: 

= = r — r. (3.76) 
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Performing the integral we obtain the well-known kink solution 



^(r) 



— tanh 



|m| 

71 



r — r] 



(3.77) 



Expressed analytically the thin wall approximation of the bounce looks like: 

$j, < r < f - Ar 



(3.78) 



$tua/i(r), r — Ar < r < r + Ar 
$/, r > f + Ar 

which describes a bubble of radius f of the stable phase $t embedded in the metastable phase 



tanh 



V2 



r — r 



describes the bubble wall separating the two phases. 
The unknown parameter f, the radius of the ball, can be determined by computing the 
action in the above approximation and demanding that it be stationary under variations of f. 
We have to calculate the action in three distinct regions according to the solution ( p.78 ). 
Outside the wall, $ = hence, S2outside = 0. 



Inside the wall, $ = $t hence, S2inside = Stt^^ {^t ^ ^f) 
Within the wall, in the thin-wall approximation. 



27rr^/i<l>/ 



m\ 



2wall 



27rr / dr 



^<i>'L, + f/+($^„,) 



- f/+($/) 
2'KfSi. 



(3.79) 



So, we can conclude that the action asscociated with the bounce in the thin-wall approximation 
is: 

(3.80) 



{S2)d = -2'Kr^h\/ ^\m\ + 2'KrSi, 



where the Si is the energy of the kink 



^1 



dr 



1 2 



\m\ 



dr- 



cosh^ f ^ 



4^2 



A 



(3.81) 



The first term on the right hand side of Eq. 



corresponds to the volume energy of a 



bubble of radius r, the second one to its surface energy. 

From the stationarity of 5*2 with respect of the variation of f (^ = 0) we can calculate f: 

|2 

(3.82 



^1 



2h\m\ 



\m\ 



6 hV3 Vx' 

For the critical bubble size a very simple expression is obtained for the exponent of the rate in 
terms of dimensionless quantities: 

52(thin wall) IQtt \m\^ 1 Att 1 



T 



^X^I^hdTd 9 hT- 



(3.83) 
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The linear dependence on 1/h is fulfilled in our numerical calculations very well, but the 
predicted action is more than one order of magnitude larger than what can be derived from 
the slope of Fig. |3.9| : 5*2 /T (measured) ~ 0.1/{hT). If one relaxes the thin wall approximation 
and solves numerically the two-dimensional bounce equation directly for several h, one finds 
52(bounce)/T = 0.78/{hT). 

We mention here that for finding the bounce solution we have to choose an initial value for 
$ > ^min{h) evolve it and see what's happening at a later x. If $ > ^max{h) then the solution 
overshoots and if 4^ < then the solution undershoots. We have to tune the initial value of $ 

ax 

and the steps of the fourth order Runge-Kutta algorithm we used, in order to achieve that the 
solution stays very close to ^max{h), as long as possible. . Given the bounce solution we can 
simply evaluate its action by performing a numerical integration . 

Further improvement can be obtained by applying the temperature corrected effective poten- 
tial in the bounce equation. We have determined the parameters of the T-dependent potential 
directly from our numerical calculation in the following way. The restoring force was measured 
for a certain h around both the stable and the metastable minima as described in Subsection 
3.2.3| . Next an interpolating fit has been constructed of the form: /igg+m^g^ + As^^ + Agg^^/G. 
It turned out that in the best fit A3 ~ is fulfilled always, while the values of the other co- 
efficients are not too far from their tree level values. Then a bounce solution can be built 
on the corresponding (real!) potential which includes the temperature corrections. The final 
result is quite close to the measured value of the rate logarithm: 5*2 (T, bounce) /T = 0.29/(/iT). 
The same real interpolation can be built on the neighbourhood of the minima of the 2-loop T- 
dependent effective potential, leading to S'2(r, bounce) /T = 0.2/(/iT). By common experience 
in the surface tension simulations a factor of 2-3 difference in S2 is expected to arise relative to 
the mean field theory. 

We conclude, that the finite temperature corrections are important for quantitative treat- 
ment of the nucleation rate within the thermal nucleation theory. 

Next we turn to the discussion of the possible nucleation threshold at small h. The average 
"release time" increases for fixed lattice size with decreasing h (Fig. |3.10| ). On a lattice of fixed 
size we found small h values for which we could not detect any transition, what makes very 
probable the existence of a threshold value of the external field hth{N). We did not attempt to 
locate this value beyond the simple hyperbolic fits to the few largest {Ueiease) values. The value 
of hth{N) remains stable when the smallest h is excluded from the fit, therefore we conclude 



that a threshold magnetic field exists similar to the case of the kinetic Ising model [£3|. The 
value of hth decreases when the lattice size is increased. Intuitively we expect h{oo) = 0, but 
with the three lattices studied by us this conjecture cannot be demonstrated. 
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Figure 3.10: The average "release" time as a function of h measured on = 64,128,256 
lattices. 

3.2.6 Conclusions 

In this Section we have investigated in detail the decay of the false vacuum in a classical lattice 
field theory. This serves as a demonstrative example of the nonlinear relaxation of a far out 
of equilibrium classical cut-off field theory. Our investigation was based dominantly on the 
effective theory of the order parameter. Two versions of the effective theory were reconstructed 
from its trajectory derived from the microscopical equations of the theory. The first refers to 
the (meta)stable branch of the motion. The second one which takes into account the existence 
of a mixed phase during the transition period describes very well the transition together with 
its neighbourhood. The first equation has the form of conventional mechanical motion taking 
place in a dissipative noisy environment. The dissipation, the dynamical breakdown of the time 
reversal symmetry was found only for times longer than the minimal microscopic time scale of 
the system, the autocorrelation time of the noise. 

Our results offer a"dualistic" resolution of the competition between the nucleation and 
the spinodal phase separation mechanisms in establishing the true equilibrium. On the one 
hand, we find that the statistical features of the decay of the false vacuum agree with the 
results obtained by expanding around the critical bubble. The microscopic mapping of the field 
configurations during the relaxation supports the bubble creation scenario. Alternatively, the 
effective OP-thcory displays the presence of soft modes and produces dynamically a Maxwell- 
cut when the time dependence of the transition trajectory is described in the effective OP 
theory. We find that the larger is the system the smaller is the external field which is able to 
produce the instability. For infinite systems an infinitesimal field pushes the system through 
the Maxwell-cut, where no force is experienced by the OP. Therefore it will not stop before 
reaching the true homogenous ground state passing by the mixed states with constant velocity. 



N=64 
N=128 
N=256 
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Chapter 4 

The damping of the Goldstone modes 



This Chapter is devoted to the dynamics of the Goldstone modes. First we derive a set of 
effective equations for the long- wavelength modes of an 0{N) symmetrical scalar field theory. 
Next, we analytically determine the damping rate of the Goldstone modes as a function of the 
wave number k. Next, in an 0(2) symmetrical scalar field theory we determine the same rate 
in the presence of a small explicit symmetry breaking parameter h. We study the behaviour of 
the damping rate as /i — > 0. 

By performing a 2 + 1 dimensional numerical simulation in an 0(2) symmetric classical field 
theory, we check the large time asymptotic power law of the relaxation of an arbitrary field 
configuration. The manifestation of Mermin- Wagner theorem is also discussed. 



4.1 Effective theory for the soft modes in 0(N) model 

In this section we extend our discussion to the broken phase dynamics of the N-component 
scalar field theory with self-interaction. This model is relevant to the dynamics of the tt — a- 
system (A^ = 4) [^] and is actively investigated in connection with the phenomenon of the 
disoriented chiral condensate PBl BTl R3, R^. The appearance, the evolution and the damping 



of the low frequency fluctuations and instabilities of the chiral order parameter are carefully 
studied in these papers. Less attention is paid to the damping of the Goldstone-modes. Though 
in the realistic case the explicit breaking of the 0{N) symmetry leads to massive "Goldstone" - 
bosons, it is of interest to see what is the intrinsic dynamics of these excitations in the "ideal" 
spontaneously broken case. It turns out that the damping of the Goldstone-fluctuations with 
frequency ko is essentially different in the respective domains 

fco < — < M < T, and — < A;o < M < T. (4.1) 
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Our conclusion is that both the on-shell dissipation time and the large time asymptotics of the 
Goldstone fields is qualitatively different in the first region in comparison to the fiuctuations of 
the order parameter. 

The physical picture is very appealing. Radial fiuctuations along the order parameter relax 
fairly quickly, what freezes fast the length of the vacuum expectation value of the $ field. The 
second stage of the relaxation consists of the slow rotation of the order parameter, which is 
described by the relaxation of collectively excited low momentum Goldstone fiuctuation modes. 

The high-temperature one-loop quantum dynamics of scalar models has been investigated 
very actively recently partly as a kind of theoretical laboratory for developing powerful calcu- 
lational methods, partly with the aim to provide answers to important questions of infiatory 
cosmology and the physics of heavy ion collisions. Our work was substantially infiuenced by 
References 0, |91|, |2|, |2[ . 

We first derive in Subsection the effective equations for the low frequency modes. The 
coefficients of the terms correcting the classical equations are determined by various n-point 
functions of the fast modes. For the effective dynamics we work out the modification of the 
linear part of the equations (the self-energy operator), which is fully determined by the two-point 
function of the high-frequency fields. A detailed study of this quantity appears in Subsection 



4.1.2| . Stated in more technical terms, we compute one-loop contributions to the two-point 
functions with full T = propagators and restrict the temperature dependent contribution 
to the po > A modes. This approach follows the finite temperature renormalization group 



transformation scheme of D'Attanasio and Pietroni In Subsection [4.1.3| we shall discuss 



in particular, whether the linear response of the high frequency modes to long wavelength 
fiuctuations allows a classical kinetic theory interpretation. Using the results obtained for 



the two-point function of the theory, in Subsection [4.1.6| we present the explicit effective field 
equations, with help of auxiliary fields introduced to handle the nonlocal nature of the effective 
dynamics. Our results are summarised in Subsection [4.1.7| . In Appendix ^ some results of the 
main text are rederived with help of the conventional perturbation theory in a form explicitly 
showing its equivalence to the Dyson- Schwinger treatment. Appendix |C| presents the iterative 
solution of the dynamics of the classical 0(A^)-model. Some relevant integrals are explicitly 
evaluated in Appendix 0. 

4.1.1 The effective equations of motion for the slow modes 

The Lagrangian of the system is given by 

^ = ^(W^-^-W-^(($.)^)^- (4.2) 
64 



4.1. EFFECTIVE THEORY FOR THE SOFT MODES IN 0(N) $^ MODEL 



The aim of this investigation is to integrate out the effect oi po > A high frequency fluctuations. 
From the point of view of the final result it turns out to be important whether A < M or A > M, 
where M is the mass scale spontaneously generated in the broken symmetry phase. In the spirit 
of the renormalization group the value of A will be lowered gradually and the importance of 
passing the scale M will become evident in this process. Our actual interest will concentrate 
on the temporal variation of the lowest frequency fluctuations ko <^ A. 

We separate in the starting Lagrangian the high-frequency modes {ipa{x) with frequency 
u; > A) and the low-frequency modes {(f)a{x) with frequency u < A) 



(4.3) 



Averaging over the high frequency fluctuations gives {(fa{x)) = which yields also (paix) 
($a(x)). This separation leads to 



i(90J^ + i(5^J^-^ 



X 
24 



24 



(4.4) 



From Eq. ( [4.4| ) the equations for the slow modes can be derived 

A 



24 



fa{4>bfb) + 40a(06)^ + 4:lfa{fb)'^ 

+4:^a{4>by + Ha{<PbVb) + ^a{^bf 



0. 



(4.5) 



The effect of the high-frequency fluctuations on the slow ones is obtained by averaging the 
equations with respect to their statistics. At one-loop level accuracy {(pa^bfc) = and one 
arrives at 

id' + m')^a + ^UkY + ^U^a^b) + ^MiVbf) = 0. (4.6) 
6 3 6 

We introduce at this point the conventional notation: 



(4.7) 



Below no summation will be understood when repeated indices appear without the explicit sum- 
mation symbol. We shall also use specific pieces extracted from the above two-point functions 
defined as 

6Aab{x,y) 2 



^abi^^y) = Aab(a;,i/)U=o' 



Ail' 



{x, y)= dz 



\<f>=0 



In the broken symmetry phase one has to separate the nonzero average value from the slowly 
varying field, and write the equation only for the fluctuating part: 



■>al 



+ 



(4.9) 
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where we have chosen the direction of the average to point along the a = 1 direction. We 
shall analyse the resulting equations for the a = 1 and the remaining a = i ^ 1 components 
separately. Since the main effect of the high frequency modes we are interested in is the 
modification of the mass term (self-energy contribution), therefore we shall restrict our study 
to the linearised equation of (f)a{x)- One has to take into account that the two-point functions 
{(Pa'^b) depend on the background (f)a{x). For the linearised equations it is sufficient to compute 
them only up to terms linear in the background. 
The two equations are 

id' + ^0^)01 (x) + ^0?(x) + + ^0i(x) 5^0^(x) + Mx) = 0, 

d%{x) + ^0,(a:)0i(a;) + ^0,(x) ^^^^x) + ^x) =0, (z ^ 1) (4.10) 



a 



with the induced currents 



Mx) = ^0Aff(x,a;) + ^(iV-l)0A«(x,a;), 
A / . I'n'i . ^ . (n^ , ^ \ , A - 



Mx) = -(Ai°)(x,x)-AS?(x,x))0. + -0A;})(x,x). (4.11) 

We have simplified the equations using the fact Aj has no 0-independent piece, since with 
zero background the correlators are diagonal. The equations ([4.1 1|) express the linear response 



of the fast modes to the presence of a low frequency background producing an effective source 
term to their classical dynamical equations. Also one has not to forget that (p is now the 
solution of the constant part of the equations: 

m20+ ^ (3AS^(x,x) + (AT - l)Ai°)(x,x)) + ^0=^ = 0. (4.12) 
4.1.2 The two-point function of the fast modes 



Varying Eg. (|4.4|) with respect to (fix) one arrives at the equations of motion of the fast modes 
in the background of 0(x): 

(9^ + m^)ipa + ^ (80a(0fe<^fe) + 4:ipa{ipbY + 40a(06)^ 



+ 8{4>m)^a + 4:M^bf + 4<^a(0fe)') = 0. (4.13) 



The equations of motion linearised in the high frequency fields can be written in the form 



A 



6 



At 



{d' + m')Sab + 2<f)a(j)b + {(pcYSab } ^b = ~ T^M^I^bY ■ (4.14) 



6 
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We introduce the notation 



L 



(4.15) 



and apply Eqs. ([4.14|) and ( |4.15|) to the fields appearing in the definition of the two-point function 
Aacix,y): 



dlAac{x,y) = -ml^{x)Abc{x,y), dlAadx^y) = -ml^{y)Aab{x,y). 



(4.16) 



In the derivation of these homogeneous equations we have exploited that {(f)a{x){(f)c{x)y(pb{y)) 
= 0. After the Wigner-transformation 



A(X,p) = J d%e'P-Aix,y), 



u = X — y 



(4.17) 



one arrives at the exact linear equations for the Wigner transforms. (The use of the Latin 
letters k,p, q is reserved for the 4-momenta and p ■ k, kq, etc. denotes their Minkovskian scalar 
products.) For the two-point functions diagonal in 0{N) indices one finds in the broken phase 
to linear order in the background (after using Eq. ([4.9|) ) 



-9i +ip-dx-p^ + M^] Aaa{X, p) + K 

where for the Goldstone bosons one has 
while for the heavy ("Higgs") mode 



^^^>i{q)Aaa{X,p- -)e 
^2^^i'^^i(l)^UX,P+^)e 



d^q 



2' 

^\-igX 



3' 



2 -^12 



Ai = A. 



0, (4.18) 



(4.19) 



(4.20) 



Here the tree level mass relations were used, as they correspond to the actual order of the 
perturbation theory. This procedure can be improved using Eq. ( ^.12| ) for 0, what is equivalent 
to the resummation of the perturbative series. 

After performing the Fourier-transformation also with respect to the center-of-mass coordi- 
nate and subtracting the resulting two equations one arrives at 



2p-kAnik,p) 
2p ■ kAii{k,p) 



-^0 I ^^'^i(^) -<1,P+ ^) - An{k-q,p- -) 



X 
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For the mixed two-point function Aji in the same steps a similar equation is derived: 



(2vr) 



A,,(fc - + |) - Au{k -q,p- |)] ^k). (4.22) 



In the kinematical region, where the center-of-mass variation is very slow one can assume 
q ^ k. (If there would be no X-dependence one would find A{k — q,P) ~ 6{k — q)). If in addition 
one restricts the variation of p in the high frequency fields by the relation <^ A ~ 



Ml' 



one can approximate the integrals in these equations by a factorized form. If the functions 
^aaik,p — |) are replaced by the first two terms of their power series with respect to k/2, one 
recovers the drift equations proposed originally by Mrowczyhski and Danielewicz |]BD[. In case 
of the diagonal two-point functions these equations are identical to the coUisionless Boltzmann 
equation for a gas of scalar particles, whose masses (M^ + X(f)(j)i{x) and + A00i(x)/3) are 
determined by the 0i field: 



p-dxAaa{X,p) + 



hdpAaa{X,p)=0. 



(4.23) 



Below we shall discuss an alternative to this expansion, which does not exploit the assumption 
of slow X-variation. At the end we will be in position to assess the validity of the assumption 



which led to the Boltzmannian kinetic equations ( [4.23| ). 

At weak coupling one can attempt the recursive solution of Eqs. ( |4.21| ) and ( |4.22| ). The 
starting point of the iteration is the assumption that the high frequency modes are close to 
thermal equilibrium, therefore one has for the starting 2-point functions 



A^{p) = 27r6{p') {&{po)+n{\po\)), 
Af^ip) = 2n6{p' - M',) (0(po) + n{\Po\)) , 
h{x) = n{x)Q{x - A) = ^ Q(x - A). 



1 



(4.24) 



Since the starting distributions are independent of X, the integrals in Eq. ( [4.18 ) factorize 
exactly and one has for the first corrections the following explicit expressions: 



A?(A;,p) 



.^A.00i(fc) [Ai°)(p + ^) - Ai°)(p - 



1 



aS?(p 



k.l 



(4.25) 



The retarded nature (ie. forward time evolution) will be taken into account with the Landau 
prescription fco + (see Appendix ^ for the equivalent formulas in the conventional 

perturbation theory). 
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4.1.3 Non-equilibrium linear dynamics of the slow modes 

The quantum improved, linearised equations of motion are (see eq. (|4.10|) ) 



(d^ + Ml)ct)a{x) + Ja{x) = 0. 

The (retarded) self-energy function is introduced through the relation 

J nix) = I d^Uaix - y) (paiy), 



(4.26) 



(4.27) 



where, after Fourier-transformation one obtains, (see eq. ( |4.25| )) 



ni(A;) 



n,(A;) 



-^(AS°'(p)-A<»'(rt 



(4.28) 



These integrals can be calculated (or at least reduced to ID integrals). Some details of the 
calculations can be found in Appendix O, here we just state the results: 



ni(A;) 



2 2 

^\ R,{k,M,)+(^] {N-l)R,{k,0) 



^ 2 



(4.29) 



The explicit form of their real parts is the following: 

oo 

ReRi{k,M)= I duj [l + 2fi{uj)) 



., 2|p||k| , ., 2|p||k 



M 
oo 



Rei?i(fc,M)= j duo 

M 



[l + n{uj))A{ 



■2ujko + k'^' 

2|p||k| 

2ujko + k^ + M^' 



-2ujko + A;2 
+ n{uj)A{ 



2|p||k| 



2ujko + k^ + M^' 



+ j d\p\ 



:i+^(ipi))^( 



2 


Pl 


|k 




-2 


P 


ko 


+ A;2 


-M2 



) + ^(|plM( 



2 


Pl 


|k 




2 


P 




+ P 


-M2 



(4.30) 



where |pp = a;^ — and 



Aix) 



1 

47r2|k 



■arth(x) 



^TT^Ikl 



In 



1 + x 

1 — X 



(4.31) 
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The expressions of the imaginary parts look as follows: 

Pc 

lmRi{k,M) ~^ 



Pc 



lmRi{k,M) 



47r|k| 



167r|k| 



e(-A;2) jdPh{P) + e(P-4M2) J dP{l + niko-P) + n{P)) 

Pc~ko ko/2 

Q++ko \Q-\ 
eiM^-k"^) j dPn{P)+ j dPfi{P) 
Q+ \Q--ko\ 



+e{P-M^) j dP {1 + fi{ko - P) + n{P)) 
Q+ 



(4.32) 



where 



|k| / 4M2 fco 



k^-M^ 



2(A;o + |k|) 



Q- 



k^-M^ 
2(A;o-|k|)' 



(4.33) 



The real part is divergent at zero temperature which cancels against the coupling constant 
counterterm contribution. Care has to be taken, however, when implementing the regulariza- 
tion, to maintain the Lorentz invariance for the pieces containing T = parts of the propagators 
A'^'^\ which is manifest in the original form in Eq. ( [4.28| ). 

It is remarkable that the usual domain of Landau damping (|kp > |A;oP) is apparently 
extended up to Mf + |kp > |fcoP- An analogous situation has been noted and interpreted 



recently for the propagation of a light fermion in heavy scalar plasma |^3[. Below we find for 
the damping of soft on-shell Goldstone-modes the same interpretation. 

The linearised equations of motion can be analysed from several points of view. One is 
the determination of the dispersion relations. This describes the quasiparticles, and physically 
corresponds to the "dressing" of an (external) particle passing through the thermal medium. 
The other point of view is the field evolution ||90|, when one follows the solution of the field 
equations developing from given initial conditions (history). 



4.1.4 Dispersion relations for on-shell waves 

The position of the poles is determined by the equation 



(4.34) 



We split ko into real and imaginary parts: k^ = u — iT, and assume F <^ a;. Then the 
perturbative solution can be written as 



lUq + Ren(u;o, k). 



Imn(L(Jo, k) 
2^0-0 



(4.35) 
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where uJq = + k^. The corresponding time dependence for fixed wave vector k and given 
initial amphtudes {dt(f){t = 0, k) = P(k),0(t = 0, k) = -F(k)) is found using Eq. ( [4.50| ) and 
taking into account only the poles 

k) = \p(k) + F(k) cos ujt] e-^K (4.36) 

L u J 

Goldstone modes. The tree level dispersion relation has no mass gap. The second equation 
of Eq. (|]2D shows that Ui{k = 0) = 0, that is no mass gap is created radiatively neither at 
finite temperature; this is the manifestation of the Goldstone theorem. The on-shell imaginary 
part, as Eq. ( 4.29| ) shows, comes from the continuation of the Landau-damping extending up 



to k"^ < M^. Going back to Appendix (Eq.f p.Sj) ), one finds that the imaginary part receives 
contribution from the collision of a hard thermal Goldstone particle with distribution n(po) 
with the soft external Goldstone wave of momentum k producing a hard Higgs-particle minus 
the inverse reaction, when a hard thermal Higgs with momentum distribution n(po + ko) decays 
into a soft and a hard Goldstone. The two contributions can be combined into a single integral 
leading to the following expression for the damping rate: 

_ 2 Af2/4|k| + |k| 

r,(k)=f^^ -i- / dph{p). (4.37) 



^ 3 J 327rk2 

M2/4|k| 

If |k| ^ Ml we can write with a good approximation 

This contribution survives the IR cut, if A < M^/4|k|. 

The result ( [4.38|) is interesting from several points of view. First, in HTL approximation. 



where we neglect all the masses, we would found Fj ~ ri(0) = 0. On the other hand the 
classical approximation corresponds to the substitution n{x) —^T/x'm Eq.( [4.38| ), which results 
(see Appendix |CD in 

AT 

pclass _ _ /4_3g^ 

247r ^ ' 

However, for very small momenta (|k| <^ Mf/4T) the correct result is exponentially small, 
deviating considerably from the classical result. The purely classical simulations therefore 
cannot reproduce the correct Goldstone dynamics in the long wavelength region, they would 
overestimate the damping rate. 

The most important consequence of the form of the Goldstone damping is the generation 
of a new dynamical length scale Mg"^ = AT/Mf. The components of the Goldstone condensate 
with longer wave-length can not (exponentially) decay, they survive for a longer time. 



71 



CHAPTER 4. THE DAMPING OF THE GOLDSTONE MODES 



To get more insight in the dynamics of the Goldstone modes we have calculated the damping 
rate r(|k|) for classical on-shell Goldstone modes of the 0(2) symmetric scalar fields propagating 
in a thermal medium of the broken symmetry phase taking into account the effect of the explicit 
symmetry breaking. We have specialised the discussion of the linear response function to the 
0(2) case just to simplify some formulae. The conclusions appear to be of intrinsic validity for 
the dynamical effect of the Goldstone phenomenon. 

The study of the Goldstone-propagation through thermal medium is of particular interest 
for the interpretation of the pion-sigma dynamics in heavy ion collisions. Several recent field 
theoretical studies have dealt with this problem ^ mainly concentrating on the evo- 

lution of the heavy (sigma or Higgs) component. Though it has been observed in that 
homogeneous pion condensates do not decay, the origins of this statement have not been fully 
explored. The effect of explicit symmetry breaking was taken into account by some authors by 
treating the pions as effectively massive degrees of freedom, without exploiting the consequences 
of the approximate symmetry [^, ^ . 

The Lagrangian of the 0(2) symmetric scalar field theory with explicit symmetry breaking 
differs from the one presented in Eq. (|4.2| ) only by an additional h(f)i term. We will follow what 
happens with the damping rate r(|k|) as the strength of the explicit symmetry breaking goes 
to zero. 



The difference in the linearised equations of motion for 02 Eq- ( 4.12 ) is that due to a 
modified equation for the equilibrium (static) value of the vacuum expectation, one obtains 

+ - + ^0 (v.?) W + ^^{^ir = 0. (4.40) 
5 2 

Also it induces the following mass parameters for the Higgs and Goldstone modes: 

M,^ = ^0^ + i, M2^ = i. (4.41) 



The effective linear wave equation for the field $2 looks like: 
{-k^+'^^Uk) = ^02(^)/^(AS?(p)-Ag(^) 



+ (^6\\ik^ f Ag)(p + A:/2)-AS?(p-fc/2) 

For k = the right hand side of this equation vanishes, what reflects the effect of the Goldstone- 
theorem on the self-energy function of 02. 

The decay rate of the 02-field is determined by the imaginary part of the self-energy function 
i.e. the coefficient of 02 appearing in the right hand side of ( 4.42| ). 

The kinematical analysis of the range of variables contributing to this integral is much more 
involved than in the case when the symmetry breaking happens exclusively spontaneously but 
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it goes along the same lines with the presentation in Appendix 0. After a straightforward but 
tedious procedure we find for the imaginary part of the self-energy function: 



Im n.Jk,,, |k|) 



3 / 167r|k| 



X 



e{-P) ds {n{s) - n{s - ko)) + e{-k^) ds {n{s + ko) - n{s)) 



+e{k^)e{M^ - M^)e ((Mi - Ma)^ - fc^) / ds {n{s) - n{s - ko)) 



+Q{k^)Q{M^ - Ml)Q ((Ma - Mi)^ - A;^) j ds {n{s + ko) - n{s)) 



with 



- G(fc^)e (fc2 - (Ml + Ma)^) j ds {1+ n{s) + n{ko - s)) 

cc- 

^ (ko {e - Ml + Ml) ± |k|^(P- M2 + M|)2-4PM| 



(4.43) 



(4.44) 



In case of on-shell propagation [k^ — |kp = M|) and M2 < Mi, only the fourth term of Eq. 
( [4.4ij ) contributes. With its help for the decay rate one finds 

■a++fco a.-+ko 

is nis) - j is nis) . (4.45) 

a.— 

For M-2 = and fco 7^ it reproduces the result of Eq. ( |4.38|) for the damping rate of a 
Goldstone wave without explicit symmetry breaking. 

For M2 7^ 0, the damping rate is a continuous function of |kp and can be expanded around 
|k| = 0: 



(A0/3)2e(Mi - 2M2 
327r|k|v/|k|2 + M| 



r(|k| = o) 



Ml 



-0 (Ml - 2M2) Ml - AMI X 

2887r ^ ^ M| V 1 2 

/?Mf\ / (3{Ml + Ml) 



exp 



exp 



— exp 



2M2 

P{MI + Ml) 



2M2 

-1 



exp 



f3Ml 



1 



-1 



(4.46) 



2M2 J J \ \ 2Af2 
In particular, it can be evaluated for T ^ M2,Mi, with a result which exactly coincides 



with the classical expression for the decay rate, arising from Eq.( 4.45 ) with the replacement 
n{s) nci{s) = T/s. 
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However, for AI2 — > the expression ( ^.46| ) does not follow the classical theory, but vanishes 



non-analytically as exp(— /5Mf /2M2). This result demonstrates that not only the self-energy 
function, but also the decay rate vanishes for |k| = 0, when the strength of the explicit symmetry 
breaking goes to zero. 

The vanishing with h of the decay rate of a homogeneous Goldstone modulation of the 
vacuum is in conformity with the expectation according to which the equivalence of the vac- 
uum states in a system with spontaneous symmetry breaking should be manifest also in the 
dynamical evolution of fluctuations and external signals: when an external signal superimposed 
on the actual vacuum realizes the transition over to an equivalent vacuum state no relaxation 
to the initial vacuum should occur (absence of the restoring force). 

The quantity T{h, |k|) arises at this order of the perturbation theory as the difference of 
the rates of two processes: the transformation of the Goldstone-wave into a Higgs-wave with 
the absorption of a (high-frequency) thermal Goldstone-particle and its inverse. From the 
macroscopic point of view our result presents the rate of transformation of a Goldstone-signal 
into Higgs-signal when propagating through the thermalised medium. In the approximation 
when one assumes that a single act of transformation is not followed by any further interaction 
with the thermal bath, it gives also the inverse life-time of the Goldstone-wave. However, if 
one anticipates further multiple transformations between the Goldstone and the Higgs forms of 
propagation, then the damping rate of the original Goldstone signal will considerably increase 
for small, but finite h due to the more efficient damping of the ko ^ T Higgs-waves. However, 
for h = the full rate still vanishes, since the Goldstone-to-Higgs transformation rate itself 



becomes zero. An analysis of A. Jakovac presents a systematic way for taking into account 
the effect of the larger Higgs-width in the Goldstone-propagator. 



Higgs modes. From Eq. ( [4. 281) one can easily read off the on-shell [k^ = k^ + M^) value of 



the imaginary part of the self-energy. This is entirely due to Higgs scattering into a soft + hard 
Goldstone-pair, and leads to 

|k|/2 

TmHi XMf .^^ ^, /• , ^.ko ^, , ^.ko 



Ti = = ^j^TTHl^^^ - ^) dPil + ni^ -P)+ n(^ + P)). (4.47) 

2ko 9D7rA;o|k| J 2 2 







In the k -C Mi limit it simplifies to 



AMi , , /I Ml, , 
Fi = — ^(AT - 1) - + n(-i) . (4.48) 



967r ' ' V 2 '2 

This survives the cut if 2 A < Mi. 

Since Mi < T we can perform a high temperature expansion. This means at the same time 
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that the classical approximation is applicable. We find 



AT 

r. = — (iV-l), (4.49) 



which is (iV - l)/2r: 



class 



4.1.5 Solution of the initial value problem of the slow fields 

Equation( ^4.26|) describes an integro-differential equation. For its solution we have to know the 
complete past history of the fields. In general we can introduce Za{x) = J^^^Qd'^yIla{x—y)(f)a{y), 
then for Xq > 

Since all singularities are in the lower half plane, (paix) defined by this relation vanishes for 
Xq = t < 0. z{x) can be used for setting the initial conditions |Q, as it can describe jumps in 



the field as well as in its time derivative. For example (f)a{x) = for xq < and = 0) = Fa, 
dt(pa(t = 0) = Pa corresponds to 

Zaix) = -Pa(x)5(xo) - Fa(x)(5'(Xo), Zaik) = tkoFaik) - Pa(k). (4.51) 

Direct substitution of Eq. (^4.501) into Eq. (|4.26| ) shows that it satisfies the effective homogeneous 
wave equation. By evaluating the ko integral for t = (for details, see below) one also can 
demonstrate that it fulfils the initial conditions set above. 

In general, we are faced with the computation of Fourier transforms of functions analytic 
on the upper complex plane. The procedure of extracting their large time asymptotics was 
investigated already in Refs. pT| , pO| . For t > the fco integration contour can be closed with 



an infinite semi-circle in the lower half-plane and we pick up the contribution of the cuts and 
poles inside the closed integration path: 



oo 



fit) = 1^ fiko) e-^'^"* = HZ{uj))e-^-* + E /'(^o) (^.52) 

— oo UJ epoles /ecuts / 

where Z^u) is the residuum of the physical poles of the /-function p = i Disc/ is the discon- 
tinuity along the cut, / = [ui, UJ2] is the support of the cut. The second term can be computed 
again by completing the original integration interval to a closed contour (see Fig. [4.1| ). The 
discontinuity itself may have poles (but no cuts!) which contribute in the same way as the 
"normal" poles. After these poles are "encircled", there are two straight contours - parallel 
to the imaginary axis - left, starting at the two ends of the cut and running in the interval 
^The poles below a cut do not contribute to this formula, they lie on the unphysical Riemann sheet. 
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+ 



= poles + 



Figure 4.1: This figure shows how to pick up the contribution of a cut (thick hne). 

[uj,ll! — zoo], where u is the end (or starting) point (threshold). After this analysis one arrives 
at the generic form 



fit) 



oo 

E iT)e-'-' I ^ piuj - zy) e-y\ (4.53) 



E(-.Z(^))e-* 
tjepoles wethresh. o 

where the — sign is to be applied for the starting, the + sign for the end point of the cut, 
and we have to take into account all poles, also the ones on the unphysical Riemann sheet. 
Expanding p around u into power series of y, the y integration can be performed. The term 
~ y^ of the expansion contributes to the time dependence t^"^^. The large time behaviour is 
dominated by the lowest power term of the expansion. If p cannot be power expanded then the 
damping for large times is faster than a power-law. 

In our case the position of the poles is determined by the dispersion relations (see previous 
subsection), and for given initial values we find for the pole contributions 



'jpoieit, k) = Z(k) P{kf-^ + F{k) COS cut 

UJ 



-rt 



(4.54) 



This expression coincides with eq. ( [4.36| ) apart from the wave function renormalization Z(k). 
The time dependence of the quasi-particle pole contribution was analysed before. 

For the cuts the factorized form can be used for the spectral function: pa = ZaP^, where p^ 
is the spectral function of the propagator: 

-2Imna 



(P - M2 - Renj2 + (Imnj2 • ^^'^^^ 

The thresholds and the leading large-time behaviour are determined by Imll, Eq. (^4.32|) . The 
imaginary part of the Goldstone self energy is non-analytic only at fc^ = 0, but also there the 
non-analytic piece vanishes as ~ exp(— M^/(2T|A;o — |k||)). This finally leads to a damping 
faster than any power which can be seen after a saddle point analysis 



■'cut A 



I3MI 



1/4 



P(k) 



sin 


k\t 




k 





+ F(k) cos|k|t 



(4.56) 



The cut contribution for the Higgs damping, on the other hand, is similar to the zero temper- 
ature case 



■'cut,! 



(t,k) 



AT 



24Mi (7rMit)3/2 



P(k) 



sin(2Mit - 7r/4) 
2Mi 



+ F(k) cos(2Mit - 7r/4) 
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Alk 

4 



967rMi vrMit 



sin( k 


t + 7r/2) 




k 





P(k) ° ^' 7 ^ ^ + F(k) COs(|k|t + 77/2) 



(4.57) 



The first term comes from tlie tliresliold of the Higgs pair production, the second from the 
threshold of the Goldstone pair production or their Landau damping. The Landau damping of 
the Higgs particles does not contribute to the power law decay. 

The terms decaying as some powers of time will dominate the time evolution after the period 
of the exponential decay for the Higgs bosons (|k|^^ <^ ^^i^)- Similar is the case for small 
Goldstone domains (|k|~^ ^ ^2"^); there the exponential decay is followed by a ~ exp(— ^/t) 
time evolution. In case of large Goldstone domains (|k|~^ ^ ^2~^)? however, because of the 
exponentially small damping rate, the situation is reversed: the ~ exp(— -\/t) behaviour will be 
dominant for intermediate times, while the amplitude is reduced by a factor of Z^^. Only for 
very long times will become the exponential damping term, arising from the Goldstone-pole, 
relevant. Its action will erase completely the large size domains. 



4.1.6 Non-local dynamics of the slow modes 



The calculation of the effective non-local evolution of the low-frequency modes can be based 
on Eq.( [4.1(j| ) using the explicit expressions for the Fourier transforms of the induced currents 
defined in Eg. ( 4. 11 ): 



+ {N-l) 



{2'KYp-k 

A202 I- 



Mk) Afhp-k/2) 
1 



36 



Mk) 



d'^p 



{2'KYp-k 
2p-k 



(27r)4 2p-k + Ml 



A^SiP + k/2) 
Mk) [Af(p-A:/2)-Ai°)(p + fc/2) 
Mk) fAS?(j,-fc/2)-Af (p + A:/2) 



(4.58) 



Implicitly in all Wigner-transforms A'^^-' the Bose-Einstein factor is understood with a low 
frequency cutoff, well separating the modes treated classically from the almost thermalised high 
frequency part of the fluctuation spectra. 

For the purpose of numerical investigations the above non-local form of the induced currents 
is difficult to use. In this subsection we discuss the introduction of auxiliary fields making the 
numerical solution of the non-local dynamics easier to implement. We shall not take into 
account the nonlinear piece of the source induced at one-loop level, since in weak coupling the 
leading nonlinear effect comes from the tree-level cubic term. The consistent inclusion of the 
higher power induced sources will not lead to qualitatively new, leading effects unlike the linear 
source, which is responsible for damping. We leave this extension for future investigations. 

We concentrate first on the induced current Ji. By appropriate shifts of the integration 
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variable p it can be rewritten as 



Ji{k) 



2 



(27r)4 



1 



2p • + fc2 2p-k-k^ 



i2n) 



1 



2p-k + k^ 2p-k-k^ 



Af (4.59) 



After explicitly performing the po integration one arrives at the following expression which has 
a clear interpretation as a specific statistical average: 



■Mk) 



dk) 



d^p 1 
(2^2^ 



;i + 2n(po)) 



1 



+ {N-1 



d^p 1 



18 7 (27r)3 2po 



;i + 2n(po)) 



2p-k + k^ 2p-k-k^ 
1 1 



2p-k + k^ 2p-k-k'^ 



P0 = (p2+Af2)l/2 

1 (4.60) 

Po=|p| 



The temperature independent piece corresponds to the T = renormalization of A, fixed at the 
scale fc^. It is absorbed into the mass term of therefore we retain in the integrand only the 
terms proportional to the cutoff Bose-Einstein factors. 

We introduce two complex auxiliary "on-shell" fields W°'{x, p), a = 1, They correspond to 
the two different mass-shell conditions appearing in the above equation, and fulfil the equations: 



(2p- A;-fc2)iy«(fc,p) =<Pi{k). 



(4.61) 



Clearly, Ji(x) can be expressed as a well defined combination of the thermal averages of these 
fields: 



Ji(x) = AV {W\x, p)) + {W\x, p)*) + 



A^-1 



9 



{{W\x,v)) + {W\x,Y>Y))). (4.62) 



Thermal averages are defined by the usual formula 

d^p 1 ^ 



{W\x,p)) 



-h{po)W%x,p) 



P0 = {p2+M2)l/2 



a = 1A. 



(4.63) 



(27r)3 2po 

For the relevant combination one can use in place of Eq. ( 4.61|) the equations arising after the 
combination W"'{x^ p) — W°'{x, p)* is eliminated: 



[{2p ■ kf - [k'^f] (iy"(fc, p) + W\-k, p)*) = 2P(j)i{k). 



(4.64) 



Equations ([4.10 ), ( 4.62 ) and ( 4.64|) represent that form of the non-local dynamics which is best 
adapted for numerical solution. 

For very small values of k^ a simplified form can be derived which coincides with the result 
of the conventional kinetic treatment of the Higgs-modes presented in Subsection |3.1.3[ One 
arrives at this approximate expression of Ji{k) if one performs an appropriate three- momentum 
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shift p— s>p=i=k/2on the p variable in Eq. ( |4.6CI| ). One finds the following expressions for the 
denominators, accurate to linear order in ko: 



2p-k±k^ 
Po 



2A;oV(pTk/2)2 + M2 - 2p ■ k± k^ ^ 2p ■ k{l ± ko/2po) + C(A;|j, 
pk 



Po =F 



2po 



(4.65) 



Performing the expansion of the denominators and of the arguments of n(po)/Po to linear order 
in k, ko one finds the following approximate expression for Ji : 



A202 



Mk) 



d^p 



(2n] 



ldn{po)kopo , 1/1 dn{po 

o + ^{Po) -PO 



jog dpo p - k Po \2 



iV- 1 1 



dn{po) kopo 
dpo p-k ' Po \2 



o + ^(Po) 



Po- 



dn{po 
dpo 




P0 = (p2+A/2)l/2 

(t)i{k). (4.66) 



The IR-cutoff A ^ |k| is important to prevent the singularity of the contribution from the 
Goldstone-modes. The contributions from the second terms in each line are independent of 
fc, therefore they are simply understood as a finite shift in the squared mass of (pi. With these 
steps one arrives at the final expression for the nonlocal part of the induced "Higgs" current: 



Ji{k) 



A202 
4 



d^p 
(2^ 



1 dn{po 



.Po dpo 



P0 = {p2+A^2)l/2 



+ 



- 1 1 dn{po 



9 p^ dpo 



Po=|p| 



V ■ k 



koMk)- (4.67) 



Here = (l,p/po)- In this case it is sufficient to introduce only two real auxiliary fields 
W^{x,\) and W^{x,\) in order to make the dynamical equations local. They are defined 
through the equations: 

V ■ kWik, v) = koMk), a = l,i. (4.68) 

The weighting factors in the integrals over the auxiliary field can be interpreted as the deviations 
from equilibrium of the Higgs and Goldstone particle distributions, (5ni(p) and 6ni{p), due the 
scattering of the high frequency particles off the (pi condensate. Then one writes 



Ji{x) 



^S{W\x, v)) + - — -S{W%x, v)) 



(4.69) 



4 V ^ ^ ' " 9 

where means a phase space integral weighted with the non-equilibrium part of the relevant 
distributions. 

For this case it is easy to construct the energy-momentum vector corresponding to the 
nonlocal piece of the induced current. One can follow the procedure proposed by Blaizot and 



^ The proper solution of the effective dynamics of the heavy field surrounded by massless Goldstone-quanta 
will require the application of methods analogous to the Bloch-Nordsieck resummation at T = |^ . We thank 
M. Simionato for an enlightening discussion on this point. 
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lancu ||5^ and investigate the divergence of the (/x, 0) component of the energy-momentum 
tensor: 

d,T,^,uce, = -Jid'K (4.70) 

The task is to transform the right hand side into the form of a divergence. Using Eqs. ([4.67|) 
and ( [4.68|) one finds the following expression: 



4 J (27r)3 
A^- 1 



1 dn{po 



{vd,)W\x,^^) 



P0 = (p2+M2)l/2 



+ 



9 



Po "Po 



Po=|p| 



(4.71) 



From here one can read off the induced energy-momentum function of the approximate Higgs 
dynamics: 



rpflO 

induced 



(27r)3 
- 1 1 dn{po 



1 dn{po) 
_pI dpo 



Po=|p| 



P0 = {p2+M2)l/2 



(4.72) 



9 pI dpo 

For the exact (essentially non-local) dynamics we did not attempt the construction of an 
energy functional, but its existence for the above limiting case hints for the Hamiltonian nature 
also of the full one-loop dynamics as expressed in terms of the auxiliary variables. 

The analysis of the induced Goldstone source, appearing in Eq.( 4.58| ) goes in fully analogous 
steps. The temperature dependent part to which a statistical interpretation can be linked is 
rewritten with help of two complex auxiliary fields V°'{x,p),a = l,i as 

d'^p 



(2vr)= 



-Hpo = (p' + M2)i/2)Reyi(a;,p) - ln{po = |p|)Re V^(x, p)! , (4.73) 
Po Po J 



where the auxiliary fields fulfil the equations 

{2p ■ k + k'^ + M^)V\k, p) 
{2p-k-k^ + M^)V\k, p) 



{2p-k + e)<Pi{k), 
{2p-k-e)(t)i{k). 



(4.74) 



In this case even for very small values of ko we were not able to derive any limiting case in 
which one could treat the time evolution of the low frequency Goldstone modes as a truly 
Boltzmannian kinetic evolution. 

The solution of the system of equations ( [4.10|) , (|4.62| ), (f4.64|) , (|4.73| ), ( [4.74|) requires the 
specification of initial conditions also for the auxiliary fields W"" + W°'*,V°'. If one uses the 
"past history" condition for the physical fields 0i(x) = const, (pi{x) = const for t < 0, then the 
linear integral equation form of the equations for the auxiliary fields and their retarded nature 
imply vanishing W"", V"' for t = 0. 
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Before concluding this Section we show that the method of introducing auxihary fields is 
equivalent with the method of mode-function expansion. We choose = 1 for simplicity. 

We start from the linearised equation for the quantum fluctuation that is also linear in 
the background 0: 

d^ + m^ + -^A ip + -#(^ = 0. (4.75) 



If the amplitude of the background fluctuations is small, we can try to expand with respect to 
it. As a first step we can stop after the linear approximation. The different orders read 

0th order: (9^ + m|)(/?o = 0, (4.76) 

1st. order: {d^ + m'\j)ipi + Xcfxfxpo = 0, 

where mjj = m? + 

The 0th order solution is a harmonic function 

V^o(^) = E ^ (^ke-^'^^ + aie^'^) , (4.77) 



where fco = = ^/k^ + rn^j. At 1st order we perform a Fourier transformation obtaining in 
k-space 

q ^ 

When transforming back in x-space we introduce the auxiliary field W{x, q) defined by 

A /" ^ V ^~'^'''^(^) = \ /" ^ V ^~''''^'^(^) (A 70] 



2n ^ {k + qY - m\ J 2n ^ + 2kq 
and we obtain 

M^) = -Yl Ke"''"W^(a;, q) + a+e'«"iy*(x, q)) . (4.80) 



q 

from this form it is evident that W satisfies 



{dl~2iqd,)W{x,q) = X<j){x), (4.81) 

that corresponds to Eq. ( ^.61| ). 

The complete solution up to linear order in (j) reads cf. Eq. ( [4. 77] ) and Eq. ( [4.80| ) 

^ = J2— (ake~*^^'(l - W{x, k)) + a+e'^'\l - W{x, k))*) . (4.82) 



k ^ 



This expression makes the connection with the mode-expansion since by looking at Eq. ( |2.84] ) 
we can identify 

Mx) = e'^^'^Xl - W{x, k)). (4.83) 
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4.1.7 Conclusions 

In this Section we have performed a complete 1-loop analysis of the time evolution of low 
frequency {ko <^ Mi) field configurations in the broken phase of the 0{N) symmetric scalar 
field theory. 

We have analysed exphcitly the case when T > Mi, which occurs, for instance, in the 
vicinity of the second order phase transition of the model. Here we could compare our results 
with the results arising from the dynamical equations of the classical 0{N) field theory with 
appropriately chosen parameters. For the "Higgs" particle we have found complete agreement 
for the |k| <C Mi < T modes. Related to this is the fact that we were able to show that 
the exact one-loop equations derived for the auxiliary fields VF" have a simple Boltzmannian 
coUisionless kinetic form for small |k|. 

However, the analysis gave different conclusions for the linear response of the Goldstone- 

modes. The agreement with the classical theory is restricted to the interval Ml /AT <C |k| <S 

M2 

Ml. Below the new characteristic scale M2 = the damping of the Goldstone modes 
becomes exponentially small and it vanishes non-analytically for |k| — > 0. This is a very 
natural manifestation of the Goldstone-theorem in a dynamical situation: no homogeneous 
ground state will relax to any "rotated" nearby configuration. In the light of this suggestive 
picture it is not surprising that we could not find a classical kinetic interpretation of the exact 
one-loop dynamics of the Goldstone modes even for very small values of |k| . Introducing a small 
explicit symmetry breaking term ~ /i we have checked that the damping of the Goldstone modes 
with wave number |k| — > vanishes non-analytically as /i — > 0. This is a convincing argument 
demonstrating that the behaviour outlined above is a manifestation of the 0(N) symmetry. 

Besides the exponential damping analysed above, {1 — Z) fraction of the initial configuration 
follows a different time evolution determined by the particle production and Landau-damping 
cuts. In case of the lowest wave number Higgs fiuctuations this leads to the result that the 
exponential regime will be followed by a power-decay for times t > Mf ^. In case of the 
Goldstone modes the cut contribution is ~ exp (—A^/M^), and for modes above the scale M2 
this will dominate for large times. Below the scale M2 the cut contribution will be observable 
for intermediate times, while the pole dominated damping, because of the exponentially small 
damping rate, becomes relevant only for very long times. 

It will be interesting to see through numerical investigations, if the onset of the nonlinear 
regime will infiuence the damping scenario of the Goldstone fiuctuations. The effect of interac- 
tion among the high-A; modes (in addition to their scattering off the low-A; background) on the 
effective theory merits also further study. 
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4.2 Nonlinear relaxation of classical 0(2) symmetric 
fields in 2 + 1 dimensions 

In this Section we wish to gain insight into the vahdity of the classical linear response theory, by 
comparing its results to exact solutions. A similar attempt was already done in |80, 0]. There 
the early, far from equilibrium time evolution was confronted with the linear response results, 
and a relaxation slower than expected has been observed. We will show later (see Fig. 4^) that 



the decay in the linear regime realizes the fastest relaxation rate of this quantity during the 
whole time evolution. Exploration of the range of validity of the linear regime is possible only 
by following the evolution very long, and averaging over initial conditions in order to raise the 
signal over the noise level. 

4.2.1 The model and its numerical solution 

We concentrate on the 0(2) symmetric classical scalar theory with a small explicit breaking 
term in its Lagrangian: 

C = ^(9,$i)^ + lid,<!>,r - - - ^ + ^')' + (^-S^) 

where h controls the explicit symmetry breaking, and < 0. We have studied the time 
evolution of a system, discretized on lattice in two space dimensions at so low energy density 
(see below) that according to the mean field analysis this would correspond to the broken 
symmetry phase. 

The dynamics eventually drives the system towards thermal equilibrium. For h ^ 0, the 
equilibrium state has large magnetisation, nearly corresponding to the minimum of the classical 
potential. Fluctuations around this state are naturally divided into Goldstone (light) and Higgs 
(heavy) excitations. These excitations experience an effective potential, which agrees very well 
with the result of finite volume perturbation theory. The relaxation into this state is the main 
subject of this investigation. We shall present a detailed comparison of the exact time evolution 
with the linear response theory. 

A second relaxation process can be initiated from the equilibrium with /i 7^ 0, if the magnetic 
field is switched off. The final h = equilibrium, however, obeys in the thermodynamical 
limit the Mermin-Wagner-Hohenberg theorem ^S\, which states the absence of spontaneous 



magnetisation in two dimensional systems with continuous symmetry. For equilibrium two 
dimensional finite volume systems even for vanishing external source, at very low temperature 



there is a non-zero magnetisation with a well defined direction as shown in |99]. The finite 



volume magnetisation selects the direction with respect to which we can define a parallel and 
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a perpendicular mode. This natural choice of base suggests the use of the Goldstone boson 
terminology in the h = case too. The way the finite magnetisation disappears as the volume 
goes to infinity will be touched upon shortly at the end of our discussion. 

The exact time evolution of the lattice system was studied by introducing dimensionless 
field variables $i y/a^i and rescaling further these fields like <|)^6|caled _ $-1^ 2a/A/6 and 
^rescaled _ jiyfxjQ^ This sets the bare coupling to A'^*'^'^^^*^'^ = 6. We have chosen m?a? = — 1, 
therefore all dimensionful quantities are actually expressed in units of \m\. Squared lattices of 
size ranging from 64 x 64 to 512 x 512 were used. 

The explicit symmetry breaking parameter h is chosen in the range ... 0.0025/ Ve. Ini- 
tially the system is at rest near the origin of the $-space. A small amplitude white noise 
determines the field in every lattice point ($i(x, t = 0) = 771 (x), $2(x, t = 0) = ri2{'x.), $i(x, t = 
0) = $2(x, t = 0) = 0, {rii(x.)rij{y)) ~ ^x.yi^jj)- This means, that the space average of $i, the 
two-component order parameter (OP) ($1^, $2^) is very close initially to zero, which is (ap- 
proximately) a local maximum of the effective potential. The roll-down towards the minimum 
is oriented by the external field h. It is the potential energy difference of the initial and the final 
states, which is redistributed between the modes. The magnitude of the noise was so small, 
that the final energy density of the system is exclusively determined by this difference in the 
potential energy. 

Numerical diagonalization of the measured 2x2 matrix of equal time fluctuations (Cj j(t) = 

V V V 

$ j (x, t ) $j (x, t) — $ j (x, t) $ j (x, t ) ) gave as uncorrelated degrees of freedom the radial (Higgs) 
and angular (Goldstone) components of ($i,$2): 

$ //(x, t) = ^<l>2(x,t) + <l>i(x,t), $g(x, t) = $^,0 arctan 1^1^, (4.85) 
^H,o = \m\>yQ/X + h/2\m\'^ + 0{h^) being the minimum of the bare potential valley. This state- 

V V V 

ment means that ^H(t)^G{t) ~ ^H(t) $g(^) for all t. The boundedness of $2 fluctuations 
follows from the IR-regulating effect of the explicit symmetry breaking term. The position of 
the tree level minimum scaled by \m\ is ^h,o = ^ + h/2 + 0{h?), and the scaled squared masses 
read as = 2 + h/2 + 0{h'^) for the radial mode and itlq = h + 0{h?) for the angular one. 



4.2.2 Relaxation to the h equilibrium 

During and immediately after the roll-down the main mechanism of the fast excitation of 
spatial fluctuations is the parametric resonance, as observed and studied by many authors, e.g. 



100, 101|. Also the dephasing of different oscillation modes can be observed as reported in 



151, H6|. The time scales of this stage are characterised by the cutoff and by the masses of the 



radial and angular modes. The measurement of these latter is discussed later. 
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Figure 4.2: Time evolution of |k| power spectra of Higgs (left) and Goldstone (right) field 
components (512 x 512 lattice, h = 0.0025/ ) 



The early time evolution can be characterised by the variation of the kinetic power spectra 
of the Higgs and Goldstone fields (see Fig. [4.21 ). After dephasing the potential and kinetic 
energy of each mode is balanced in itself, therefore the energy content of different modes may 
be characterised by the kinetic power spectrum 

Ph,g{\^\) =^''\ (4.86) 



where the averaging should be taken over the polar angle in momentum space. In the final 
equilibrium state this quantity does not depend neither on |k|, nor on which field it refers to. 
Its value corresponds to the temperature, which was measured in our case to be Ph,g{\^\) = 
rpkm _ o.l25±0.0005. This kinetic temperature can obviously be defined also out-of-equilibrium 

as T'^^" = ^jj^^ ^ 

The initial parametric peak in the power spectrum of the Higgs field disappears rapidly. 
It proceeds mainly through decays into pairs of the lighter Goldstone modes. This damping 
is reflected by an energy transfer from the Higgs to the Goldstone modes. In addition, the 
excitation of spinodally unstable Higgs modes also contributes to the relaxation, as it was 
shown in Subsection pj.2.2| . 

The parametric low |k| Goldstone peak, however, survives as long as t = 10^ . . . 10^. During 
its long relaxation the energy is slowly transferred back to the Higgs modes. This process 
qualitatively corresponds to off-shell emission of Higgs waves. The final equilibrium is reached 
in a slow, non-exponential relaxation of T^*" to T^™. 

The mechanism sketched above is suggested by the structure of the analytic formulae de- 
rived in the perturbative relaxation analysis, which should be relevant at least to the large 
time asymptotics of the evolution. The forward time evolution of an initial configuration is 
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determined by the classical self-energy function n(x, t): 



dkn 



2n k"^ -mjjQ- IlH,G{k) 



-ikfjt 



(4.87) 



where z{k) = ikoF{k.) — P(k) is determined by the corresponding initial configuration -F(k) = 

$(t= o,k), p{k) = dMt = o,k). 

For the calculation of the Fourier transform of n(x, t) the procedure described in |24, ^ 
was used. (For the quantum treatments, which lead for large occupation numbers to the same 
results, see [|10|, ||). 

First one makes the shift $1 — $1 + $ in order to describe the constant equilibrium 
background, where $ is the classical ensemble average of $/f(x). Its direction is actually selected 
by the explicit symmetry breaking. Then in the linear approximation one finds $h(x) ~ '^'i(x) 
and $g(x) ~ $2(x)$i/,o/^- We refer to Appendix for the detailed derivation of the self 
energies, which read now for the Higgs and Goldstone fields, as follows: 
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sin UJG,qt cos iUH,pt 
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(2vr)2 
sin ^JH,^t coscug.p^ 



UJ. 



G,p 



(4.89) 



where ^^/^ ^ 



m 



"IjIq + |kp. The relevant mass values were calculated numerically from the 
relaxed field configurations with the method developed in Subsection 3.2.3| . At the parameter 
values specified in our simulations the typical masses are amu 1.3 and aniG ^ 0.033. 

The imaginary part of the self energy, which accounts for the damping phenomena, may 
be extracted using the principal value theorem. There is a technically important difference 
between our present formulae and the ones derived in Subsection 4.1.4 . It follows from the 
fact that the system investigated in this section is defined in two space dimensions. This 
circumstance modifies the volume element in the momentum space and an integrand of more 
complicated analytical structure appears when calculating the imaginary part of the classical 
self energy function. The structure of the cuts in the Higgs and Goldstone self-energy functions 



is presented in Fig. 4.3 



The explicit expressions of the imaginary parts are the following: 
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Figure 4.3: The analytical structure of the classical Higgs and Goldstone self-energy functions. 
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(4.92) 



fc4 + 4|k|2m|^C' 



It is the off-shell damping coming from the contributions of the cuts to ([4.87|) which has a 
direct impact on the late time evolution of the OP we are analysing. We followed the method, 
developed in |102|| , which determines the leading power law tail of the relaxation. 

Contrary to the 3-1-1 dimensional case in 2 -|- 1 dimensions the integrand of Eq. ( |4.87| ) 
has extra cuts. These can, however, be directed in a way that either they do not contribute 
(Goldstone case) or their contribution is suppressed by a factor of e"™^* (Higgs case). The 
large time asymptotics of the two field components is the following: 

1 fP{k) cos(tfiH,k) 
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P(k) cos(tf]G,k) 
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F(k) sin(if^H,k) 
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(4.93) 



(4.94) 



m 



H/G 



and 



where the threshold frequency values, appearing above are given by ^"jj/G k ~ 
^H±G,k = + {niH ± mG^. 

In order to obtain numerical evidence for the presence of the predicted power law tail, we 
measured the quadratic spatial fluctuation moments of both fields as a function of time. They 

Y / y\ 2 

are defined as Fluct{t) = $(x, t)^ — (<l>(x, t) j . Substituting ( |4.94| ) into this definition we 
get ~ like relaxation for late times. This behaviour is indeed observed as the fitted power 



shows in Fig. 4.4 for the example of the Goldstone mode. 




7 8 9 10 11 12 

log(t) 



Figure 4.4: Relaxation of Goldstone spatial fluctuations with time on (e-based) log-log scale 
plot. (256 X 256 lattice, h = 0.0001/^6, tug = (62 ± 5) ■ 10-^ = 1.329 ± 0.001 ■ 10-^, the 
average over 69 runs is shown.) 



The formulation of the theory upon which the above formulae were obtained actually uses 
ensemble averages over the initial data of the classical evolution. We found, that although 
every single run seems to relax even quantitatively in the same way, this late time evolution 
can be extracted from the noisy data only, if this averaging is performed indeed (116 runs 
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were involved). We could recognise a ~ decay for Higgs modes too, but because of the 
low signal/noise rate we could not do quantitative analysis. In particular the analysis is made 
difficult by the fact, that Higgs modes decay rapidly and vanish in the noise before the power 
law tail would be reached. 

Relaxation behaviour of OP ($i:/(k=0)) was also investigated, and — in accordance with the 
expectation above — an oscillation damped by ~ ^-i±o.o2 ^^^g fQ^^d. Moreover, this oscillation, 
is observed around a value, monotonically approaching its equilibrium value as ~ This 
latter behaviour is explained by the fact, that the exact equation of motion for OP contains 
Fluct{t) as a time dependent parameter ||100| , |103|| . Its ~ like behaviour is inherited by the 



slowly varying part of OP. (The errors of the exponents come from averaging over 92 runs on 
a 256 X 256 lattice.) 

4.2.3 The h equilibrium 

When one follows the evolution of the system for long times {t = 10^), there seems no further 
relaxation to take place, i.e. the initial signal has been lost in the thermal noise. In order to 
convince ourselves of having arrived to thermal equilibrium we compare the measured values of 
and the masses mc and rriH to the analytical estimates of these quantities coming from 
the one-loop lattice effective potential evaluated for the same size lattice as used in the simu- 
lation. In the analytic expression we used the measured kinetic temperature. The equilibrium 
masses rriG and rriH were determined in the simulation both from the corresponding correlation 
functions and by fitting the oscillatory motion around the equilibrium of the corresponding 
OP-components as described in Subsection p.2.3| and [|103[ . 



We have checked that the system with the present initial conditions is deeply in the Coulomb 
phase i.e. the temperature was about four times smaller than the Kosterlitz-Thouless critical 
temperature Tkt- In the absence of explicit symmetry breaking on the critical line between 
T = and Tkt the Goldstone correlation length is expected to diverge with the lattice size. The 
explicit symmetry breaking parameter h acts as an infrared regulator. Indeed, the measured 
Goldstone correlation length is found to be proportional to L, but for the largest sizes (L = 256), 
where the IR cutoff h begins to dominate. On the other hand, the inverse correlation length in 
the Higgs direction is finite, its value being 3% smaller than the two-loop mass and 6% smaller 
than the one- loop value. 

A very good agreement of with $ is found where 5> is given by the minimum of the 
effective potential in the radial direction, the relative deviation being (9(10~^). The discrepancy 
between the measured (with the method described in Subsection p.2.3| or ||103|j ) and calculated 



masses was less than 1% and 5% for the Higgs and Goldstone modes, respectively, the pertur- 
bative values being systematically smaller. Our measurements of the Goldstone mass became 
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Figure 4.5: Numerical and analytical values for Goldstone on-shell decay rate as a function of 
|k| (128 X 128 lattice, h = 0.0025/^6, mo = 0.0319, ttih = 1-33, $ = 0.96, T = 0.125) 



very noisy for small values of h < 0.00001/^6. 

In the equilibrium system one can proceed to "experiments", which check the correctness of 
the on-shell decay rates computed in the linear response theory. These are the simple zeros of 
the denominator of Eg. ( [4. 87]) which determine the exponential damping of on-shell excitations. 
The damping rates are obtained by substituting |k| = into Eqs. ( [4.9U| ) and ( |4.91| ). To leading 
order in (A$/6)^T the rates read (assuming rriH > ^mc) as 
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(4.95) 
(4.96) 



On equilibrium configurations single k-modes have been superimposed with an amplitude 
between 0.01 . . .0.2. Time evolution of these excited modes showed perfect exponential decay. 
The exponents did, however, depend on the amplitude. For small amplitudes the fit was 
unstable, the mode was quickly lost in the noise. For bigger values nonlinear effects did show 
up. Therefore we were looking for a plateau in the amplitude dependence of the decay rates 
between these two extremes. Then its value has been compared to the analytical estimates 



(see Fig. [4.5|) . The error bars show how well-defined plateau could be found. This numerical 
experiment was performed both for Higgs and Goldstone waves, but the errors for Higgs decay 
were too large again. The value of Th computed from the linear response theory was reproduced 
within an accuracy of ±50%, but we can not say anything about its predicted k-dependence. 
The reason for this seems to be, that the evolution becomes very soon nonlinear as we try to 
increase the excitation amplitude higher than the noise level. 
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As was indicated in ||104|| finite size effects may be highly important in systems with 
Goldstone-hke excitations. This circumstance made necessary the use of finite volume per- 
turbation theory in the analysis above. Neither the lattice summed perturbative effective 
potential, nor the numerically found mass values showed any L dependence if L > 64. We have 
checked for the L independence in all cases where the exact evolution and linear analysis were 
confronted (e.g. the values of damping rates). 



4.2.4 Diffusion to the h = equilibrium 

The most care must be taken when the macroscopic magnetisation M{L) is considered in a 
system with h = 0. For finite volume, the ensemble average computed from the low temperature 



spin- wave approximation is non-zero [99| even in the absence of explicit symmetry breaking: 

M^{L) = (^'')' + {^"^y ~ (2L)-^/(2-*'). (4.97) 

For each realization of the canonical ensemble observed in Monte Carlo simulations the mag- 
netisation vector has a well-defined direction 9 = arctan (^^2^ /^i^^ 

It is an interesting question, by what mechanism the MWH-theorem is realized as L ^ 00. 
It was shown in Monte Carlo simulations that the erasure of the order is realized by the diffusion 



like displacement of the direction of ^ p9 



In our real time study, for the investigation of the onset of the finite volume version of MWH- 
theorem a number of runs were continued from the equilibrium state reached with h ^ 0, after 
the magnetic field was switched off. In each individual system OP begins to circle around the 
origin with an average radius M{L). This motion is probably due to a random nonzero angular 
momentum of the initial configuration. Subtracting the angular momentum of OP, the relevant 
diffusion-like motion remains. MC studies employing first order Monte Carlo "time" -evolution 
of non-equilibrium one dimensional systems have shown that the sign for SSB is the exponential 



decrease of the diffusion constant with L ||105|| . In agreement with the expectations based on the 
MWH-theorem we find in the present case that the diffusion constant decreases with increasing 
L following a power law with the exponent — 1.16 ± 0.1. 



4.2.5 Conclusions 

In conclusion we can state that the numerical study of the Hamiltonian dynamics of the classical 
0(2) symmetric scalar model in 2+1 dimensions provides a non-trivial check of our understand- 
ing of real time relaxation phenomena. Numerical results both for the late-time OP-asymptotics 
and for the decay rate of on-shell waves (with well-defined wave vector k) were found to be 
in agreement with the linear response theory. In the h ^ equilibrium the agreement of 
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the masses coming from the perturbatively calculated effective potential with the numerically 
established excitation masses was also verified. Finally, we have demonstrated the real time 
manifestation of the finite volume MWH-theorem by measuring the large L asymptotics of the 
angular diffusion rate of the macroscopic order parameter. 
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Chapter 5 

Real time dynamics of the Higgs effect 



The real-time transformation of light gauge fields into massive intermediate vector bosons is a 
subject of increasing cosmo logical interest. 

Recently, Garcia-Bellido et al. ||12[ (see also [|1^) assumed that the coupled inflaton+Higgs 
system was far out of equilibrium when the energy density of the expanding Universe has passed 
the point corresponding to the thermal Higgs transition. With the supplementary condition 
that the reheating temperature after the exit from the inflationary period did not exceed the 
electroweak critical temperature, they demonstrated in a (l+l)-dimensional toy model, that 
large enough matter-antimatter asymmetry could have survived till today. The existence of 
non-equilibrium Higgs transitions in 3+1 dimensions has been demonstrated by Rajantie et al. 
106|| . In these investigations the classical mean field equations were used with renormalised 



and temperature-dependent couplings. 

Traditionally the evolution of the baryon asymmetry through the electroweak phase tran- 
sition, accompanying the onset of the Higgs effect has been thoroughly discussed, under the 
assumption that the system is in thermal equilibrium, with specific emphasis on the high- 
temperature sphaleron rate ||107| , |108| , |109|| . Near equilibrium, the gauge-invariant HTL action 
dominates the influence of the thermalized quantum fluctuations on the motion of the mean 
fields with k ~ 0{gT) [|110| , |111| , |112| . This term was taken into account in real-time simulations 
of the sphaleron rate [|108|| (see also ||113|| ). 



The framework for the theoretical study of the baryon asymmetry has been somewhat 
modified with the realization that, within the Standard Model, the Higgs effect sets in via 
smooth phase transformation, characterised by an analytic variation of the order parameter 
P, |114] , |115| , |116|| . Under such circumstances the expectation value of the Higgs field is different 



from zero also above the electroweak energy scale. In this context it is remarkable that the 



leading HTL correction is insensitive to the presence of the scalar condensate ||117| , |108|| . 

In this Chapter we go one step beyond the HTL correction. Generalising the discussion of 
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the damping rates of scalar fields in the previous chapters, we complete the kinetic theory of 
the pure gauge fields by appropriate terms reflecting the presence of scalar fields. Our aim is 
to find the leading effect due to the presence of an arbitrary constant background field: <!'. In 
equilibrium, such a background is sufficient for the calculation of the effective potential (see 
for example fT^), but not the full effective action. Similarly here, we restrict somewhat the 
generality of our correction terms to the equations of motion, namely the corrections to the 
$^-y4^ vertex will be determined neglecting the non-local effects in the scalar field. 

The high-temperature limit of the resulting expression of the induced current shows that the 
next-to-leading "mass and vertex" corrections can be uniquely decomposed into gauge-fixing 
invariant and gauge-fixing dependent modes. In this way we can propose an expression of wider 
applicability for the Landau damping effect also including the effects of a scalar condensate. We 
propose a physical gauge for the numerical solution of these equations where the gauge-fixing 
dependent mode does not propagate. 

Baacke and his collaborators | |118| , |119| , |12U|| have applied, in a series of papers, a comple- 



mentary approach to the Higgs effective action without assuming the presence of a general 
gauge background. 

A non-trivial mean gauge field configuration reacts back on the scalar condensate. Our final 
goal is to derive a coupled set of equations for the scalar condensate and the mean gauge field, 
which is appropriate for studying the real-time onset of the Higgs effect and the variation of 
the sphaleron generation and decay rate. 

5.1 Mean Field Equations in the Abelian Higgs model 

The Lagrangian of the model in the 0(2) notation is the following: 

C = -\f,,F^'' + \{d,^f - Kn"^^ + e [^2d,^i - ^18,^2) > + - ^A^^^f^ (5-1) 

where $ = ($1, $2). We split the fields into a mean field (A^, $) and a fluctuation contribution 

A^ = A^ + a^, $ = <!> + (a^) = ((p) = (5.2) 

at any time (averaging is understood with respect to the initial density matrix). We assume 
that the scalar background is constant and that it points to the $1 direction, its value being $. 
We use in the sequel the notations my/ = e$, = + A$^/2, = + A$^/6, although 
these are not the vacuum values. 

We fix a 't Hooft gauge by changing the Lagrangian as 

C - ^{d^A'' + imw^2f - c((9^ + + e^mw^i)c, (5.3) 
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where c is the ghost field. In this way the field (p2 receives an additional mass contribution 



w 



The operatorial equation of motion (EOM) separately induces EOMs for the mean fields 
and the fiuctuations. The average of the operatorial EOM for the gauge field reads as 



A^{x)+f;\x) = Q 



(5.4) 



with the induced current 



where j^(a;) = ip2{,x)d^j_ipi{x) — ({)i{x)d^({>2{x). We will perform the calculations at the one-loop 
level, when the last term does not contribute. The subtraction of ( |5.4| ) from the full equation 
yields the EOM for the fiuctuations. At one loop it is sufficient to consider only the equations 
linearised in the fiuctuations. On the other hand, since we want to calculate f^'^ in the linear 
response approximation, we can neglect all terms non-linear in A. These assumptions make the 
equations very simple: 



-(9^ + m^)^^ 



a:'{x) - 2e^^Lpx{x)A^,{x) = 0, 



(9^ + ttih) v9i(x) + 2eA{x) ■ d(p2{x) + eip2{x)dA{x) - 2e^^A{x) ■ a{x) = 0, 
(9^ + ttiq + ^m^) (p2ix) — 2eA{x) ■ d(pi{x) — e(pi{x)dA{x) = 0. 



(5.6) 



The ghost fields follow a free EOM, so that they do not infiuence the present calculation. These 
equations are solved to linear order in the A background: 



a,{x) = af{x)-2e'^jd'zG^^^{x-z)A^{z)vf\z), 



ipi{x) 



ipf\x) + e / (fzG'^ix - z) 2A>'{z)d^^'^'{z) + {dA){z)^'^> {z) - 2e^A^'{z)a'-;^\z) 



(0), 



,(0)/ 



Mx) = vf\x)-e Id'zG^ix-z) 2A'^{z)d,v^;'\z) + {dA){z)iff\z) 



(5.7) 

where the superscript zero denotes the solutions of the free EOM, and G^'s are the free retarded 
Green functions^. 

For the induced current ( |5.5| ) we need the expectation values of certain local products of 
the fiuctuating fields. For example, it directly follows from the above equations (with {AB)q = 
(^(0)5(0))), that 



{Mx)df^Mx)) 



eld'zl dlG^{x-z) 2A^{z){^2{x)dl^2{z)), + {dA){z){^2{x)v2{z)), 



G^{x-z) 



2A 



^{z){d,v^{x)d,V^{z)), + {^A){z){^,^^{x)^^{z)),^^ }. (5.8) 



^Here the definition KG^ = —5 is used, where K denotes the free kernel. 
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We define the local products in a symmetric way (i.e. ^{ip2{x)d^ipi{x) + d^(pi{x)ip2{x))) and 
introduce 

-z) = -{<^{x)<^{z) + <^{z)<^{x))o. (5.9) 
Then the above expression can be written in Fourier space as 

{^2d,^i){Q) = -eA'^iQ) I ^ p,{2p - Q\ [Gf (p) A2(g -p) + G^{Q - p)A,{p)] . (5.10) 

The evaluation of other expectation values goes along the same line, finally giving 

e{j,){Q) = -e'A^Q) |^ (2p - Q),{2p - Q), [Gf (p)A2(g - p) + G^{Q - p)A,{p)\ , 

2e'^{a,^^){Q) = -\e'm%,A\Q) J ^ [G^^MQ - P) + G^{Q - p)A,^{p)] , 

e'A,{^') = e'A,{Q) J ^ [A,{p) + A^ip)] . (5.11) 
Assuming that the free fiuctuations are in thermal equilibrium, the propagators can be related 



to the corresponding spectral functions ||121|| , which are the discontinuities of the free kernels 
of eq. (|5.6|). Introducing 

A^.{p) = 2n(l + n{\po\))6{p'-m'), and G^P) = (5-12) 

\ 2 / p'' — + lepo 

where n is the Bose-Einstein distribution, we have 

^i = ^m%^ ^2 = A„2 +^„2 , A^^ = -5(^j,A„2 + [ A^2^ - A^„2^ j , (5.13) 

and analogously for the corresponding G^'s. 

5.2 High-temperature expansion beyond the HTL ap- 
proximation 

The leading HTL term of j^"*^ in the high-temperature expansion comes from e(j^) by neglecting 



all the masses ||117|| . Our aim is to calculate the first subleading term proportional to instead 
of in the high-temperature expansion. 

We emphasise that j^^'^iQ) = J*"'^(— Q)*, therefore all partial contributions that are odd 
under hermitian Q-refiection can be freely omitted. 



A useful formula, which simplifies the calculations, is the following ||122|| : 

/^(Q,m^M2) = l0-F{Q,p){GUp)AM2{Q-p) + Gfj.{Q~p)A^.{p)) 

d^P r.^^ Q „,A„2(p- f) - AM2(p + f) 



FiQ, ^ - P) \' '..2 (5.14) 



(27r)4 '^'2 2pQ + m2 -M2 
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The tensorial structures appearing in when cast into the form ( p.l4| ) imply the appearance 

in F of the p-dependent terms p^Pu,PfiQu + PuQfj. and also of terms independent of p. 

5.2.1 Even terms 

We start the discussion with the terms even under p-reflection and introduce the notation 
f{p) = {F{Q, Q/2 — p) + F{Q, Q/2 + p))/2. Here the — term in the denominator can 
be neglected since in the mass expansion 

/(p)(M^-..')^°'^-!'-e°'''^''-0 (5.15) 



(2ir)4 " " ' (2pQ)2 

because of the odd p — > ~p behaviour of the integrand. Then we only need to expand the 
difference of A's of the numerator. The numerator can be expanded with respect to Q, when 
low-momentum mean fields are considered. The leading term of the gradient expansion gives 
zero, again because of the p-odd integrand. 
The first non-zero contribution is therefore 

Now, we treat separately the two actual cases: / = Pf^Pu and / = constant. The / = p^p^ 
case contributes to the full expression of 

2e'A%Q)Qf [ ^ A(Ai + A2 + - A.„2 ). (5.17) 

We introduce the field-strength tensor by A^Qp = —iF^^ + A'^Q'^ . The local term proportional 
to A reads as 

After partial integration the first two terms cancel with e^A^l^ip^) in the induced current. What 
remains is a contribution vanishing at zero m^. In the mass expansion they give 

= _e^(l - i),-n^AM (5.19) 

For each Am2 the field strength contribution of (|5.17| ) is rewritten with the help of the relation 



s' 



= -2ps— + s„^2.S(p' - m') (5.20) 



in the form 
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The first term drops out because of tlie antisymmetry of F. In the second term we perform the 
mass expansion. After adding the different contributions we find 



(1) = _4,e2F' 



(27r)4 2pQ dpo 



U^.A". (5.22) 



The first term of this expression is the usual HTL contribution, the second is a mass correction 
to it. In the mass correction, originally, there was mf + + — C.m'^, but because of 
ml = rriQ + C,m^ the ^ dependence drops out. 

The / = constant contribution to the induced current appears as 

-4e^»f.A.W)«'/|fl^^. (5.23) 

Since it is already proportional to we have dropped the mass dependence of the integral. 
Using once more the relation (|5.20|) we write for it 



The first term is again local: 



jr" = 4eW„>l,W)/|^j^. (5.25) 



The second term reads as 



5.2.2 Odd contributions 

Odd contributions come exclusively from {a^(pi). Using ( |5.14| ) and denoting the current by j^^'* 
we find 

Finally, with the help of ( p. 201 ), performing the mass and external momentum expansion ac- 
cording to the method followed in the previous subsection we arrive at 
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5.2.3 Linear response and physical modes 

The full induced current is the sum of the different parts coming from (|5.19| ), (|5.22|) , (|5.25|) , 
O) and (lOSD : 

^ur, ^ ^local ^ + j(2) ^ ^ (5-29) 

wllCrG j^^'^^^ jlocal,! _|_ jlocal,2 

Here jj^^ and jj^'^ are ^-independent, while and j^^"* depend on the gauge fixing. For 
the physical characterisation of the system (for instance, damping rates) we have to find the 
independently evolving modes. In order to do this we decompose the polarisation tensor in the 
tensor basis, appropriate for finite-temperature studies | |123| |: 

P% = ^. 5.. = - + , (5.30) 

where Q = (go, q), = q'^ and uj^ = g^o - Q^.qo/Q'^ was used. 

Since II^i, of ( ^.22|) is transverse {Q^Tlfj^u = 0), it is the combination of P^ and P^ . Intro- 
ducing = Hoo and Ut = 1/2P^^W we find 



n = --^n^p^ + n^p^. (5.31) 



The Z^i, term in (|5.28|) can be written as a combination of P*^ and S. Introducing Zq = 
and Zs = {qlZc - Q'^ Zqq) / {qqo) we find 



Z = ZgP"^ + ZsS. (5.32) 

Using the completeness relation P^ + P^ + P^ = g the Fourier transform of the EOM (|5 
has the following form: 



A = 0, 

(5.33) 



{Q^-R- Ut)P^ +(q^-R+ ^Ul^ P^ + Qq=^ -R-Zg]p''-ZsS 

with R defined from m'^Afj, + + jj^^ = RA^. 

Since 5* is not a projector, but mixes the subspaces belonging to P^ and P*^, in this two- 
dimensional subspace the inverse propagator matrix still has to be diagonalized: 



Q^-R+fU, 



IQ^-R- Zg 



(5.34) 



The eigenvalues should be found with one-loop accuracy, which means that terms proportional 
to the square of one-loop corrections (i.e. n|, Zq and Z|) should be neglected. This, however. 
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implies that the eigenvalues are the diagonal entries. With appropriately rotated projectors in 
the (P*^, P^) plane, we can therefore write ( |5.33| ) 



{Q^-R- Ut)P^ 



Q^~R + ^n^^ + Qg' -r-Zg]p'' 



A = 0, 



(5.35) 



As expected, the transverse modes plus the Pl mode, which might be called longitudinal, can 
be made independent of the gauge-fixing parameter (see the discussion on the renormalization 
of R below). 



5.2.4 The infrared separation scale 

The p-integrals in all terms of j^""*^ are factorized into a radial and an angular integral. It is well 
known that the HTL term describes the dynamical screening of the electric fluctuations below 
the Debye scale, an infrared separation scale ("IR cut-off") has to be therefore introduced into 
the radial integration at p = M = Cm x eT. The effective equations one arrives at in this way 
are to be used for the modes p < M. 

Applying this cut-off to the integrals appearing in the expressions of R{Q) and H^u, one 
finds 



R = m^ 



w 



■ 3 + ^2, e'^T qo qo + q 

1 H T^e m m 

87r2 A 2tt^M q go - g. 



2 I ™2 I _2 



rr 2 go 

^2g 



go _ go + g 

g 2g2 go - g 



+ 4^^"^^ + "^^ + "^^)f^"|3f' (5-36) 



with K = 27rexp(— 7^;) and mj^ = e'^T^/3 — c^MT/tt^. The logarithmical UV divergence in R 
can be absorbed into the renormalization. With appropriate renormalization scale /j, = kT 
the ^ dependence can be made to vanish. 

The intermediate IR scale M contributes a term to the Debye mass, which depends linearly 
on M. It will be cancelled by the linear M-dependence of the self-energies, which shows up in 
the (one-loop) solution of the classical effective EOM (jSj) ^ g |lT|, |2^. 



The mass corrections in ( ^.35| ), proportional to ^'^A'^, should be interpreted as coming from 
an induced (non-local) $2^2 yeitex correction. In the subsequent classical time evolution it 
contributes to the self-energies of both fields. When (perturbatively) combined with the ~ MT 
terms coming from the tadpoles, it results in a finite contribution of order T"^, independent of 
the exact choice of the coefficient Cm in the expression of M. Therefore, these terms can by 
no means be neglected. 
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5.3 Discussion 

In this Chapter we have investigated the induced current in the effective EOM of the gauge field 
in the Abehan Higgs model. The main result is the determination of the gauge polarisation 
tensor beyond the well-known HTL expression in presence of a constant scalar background. 
The expressions appear in ( |5.35|) and ( ^.36| ). 



The most important property of the corrections is that, just as the leading HTL term, in 
the high-temperature expansion for the physical degrees of freedom they are independent of 
the gauge-fixing procedure (after appropriate renormalization). 

The result is compatible with gauge invariance if the pieces which are proportional to <l>^ are 
considered as the simplest manifestation of the nonlocal corrections to the Y(y4^(a;)$*(a;)$(x)) 
vertex appearing in Though for constant scalar background they break the gauge sym- 

metry, it is rather easy to embed them into a nonlocal gauge invariant structure. For example, 
of ( |5.26|) follows from the functional derivation with respect to the vector potential of 



the expression 

-2e^y" c/^x(D^<l>(x))*X>(D)D^$(x), (5.37) 

where is the operator of the covariant derivation and the operator £'(-D) is the obvious 
covariant generalisation of 

^'^2vr5(p^)^4^-^. (5.38) 



(27r)4 ' dpp 2p-d' 

Up to partial integrations Eq. (|5.37|) is a unique gauge invariant completion of ( ^.26|) and similar 



(s) 

type of completion can be written down for jli . For the determination of the finite temperature 
effective theory in imaginary time the same strategy was used by two of us successfully [ 124| . 



For us the only important point is that the corrections to the HTL-current are compatible 
with the gauge invariance of the full theory. The piece we explicitly picked up in the linear 
response approximation is certainly dominant for slowly varying scalar background. 

For consistency we have applied an IR cut-off M = Cm x eT to the fluctuations, and the 
effective EOM is valid below this scale. The subsequent 3D time evolution will be insensitive to 
the accurate choice of M. Partly it is cancelled by the 3D (Rayleigh- Jeans-type) divergences, 
partly it yields also non-zero, M-independent contributions, when the gauge-scalar vertex 
corrections are combined with the 3D would-be divergences. The order of magnitude of the 
tadpole contribution, remaining after the 3D divergences are canceled is (e^T/M) x (e^MT) ~ 
(9(e^T^), if the mean field masses are negligible (for example near the phase transition), relative 
to the effective cut-off. If the temperature decreases, the mean field masses might dominate 
over the momentum cut-off, and one has the estimate e^T$^/M ~ e^<l>^. That is, deeper in the 
broken phase the terms proposed in this note become more important. For e$^/r^ 
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very idea of this scheme breaks down. 

We shortly discuss here a possible strategy for the numerical implementation of the effective 
gauge field equation. The most convenient gauge fixing seems to be the Landau-gauge, where 

(3) 

the ^-dependent mode does not propagate. Then j/t can be left out of the discussion and 
the effective equations of motion are to be solved under the constraint d^A^^ = 0. It can 
be implemented, for example, by solving the equations for Ai and computing Aq from the 
constraint. 

The non-local induced currents are written in local form with the help of auxiliary fields. 
With the well-known form of the leading HTL current, we can write in this form also the 
corrections due to the non-zero scalar background: 



Jl \x) ='^D I ^ [v^Vj - ] ^ v), (5.39) 



.(2)/ X e^rri^T f rfi^t, (2). x 

where the auxiliary fields satisfy 

{do - vd)wl''^ = Fo„ (^0 - vd)wP = doA. (5.40) 

In the derivation of the expression of jj^^ we have performed a partial ^-integration in the part 
of (|5.22| ) proportional to 5'(p^), and used the fact that, to the accuracy of our calculation, 
Q^F^^{Q) = can be exploited in the expression of the induced current. The auxiliary fields 
account for the energy flow from the low frequency electromagnetic field towards the high 
frequency {k > M) modes, due to Landau damping. They represent generalisation of the 
well- understood effect in electromagnetic plasma |^ to the case of Higgs-effect. 

The equations for the scalar field receive mainly local corrections (renormalization and T- 
dependence of the couplings in the classical equations). However, for consistency also the 
^-derivative of the ^^-A^ vertex correction (|5.37|) should be introduced. 



The logics of the derivation followed in the Abelian Higgs model seem to be robust enough 
for us to attempt its generalisation to the non-Abelian case, which will be the subject of a 
future study. A numerical study of the corrected EOMs can lead us to a deeper understanding 
of the reliability of our present views on the high-temperature Higgs models, especially where 
the sphaleron rate is concerned. We will also learn in more detail of the non-equilibrium aspects 
of the cosmological onset of the Higgs regime. 
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Appendix A 

The introduction of the one-particle distribution function and the calculation of 
the damping rate 



In the relativistic kinetic theory (see for example ||125| , |126|| ) the number of world lines AN{x,p) 



representing particles with four- momenta = meffdx^{T)/dT which intersects at point x a 
surface segment A^S of a space-like three-surface in Minkowski space is given by the relation 

AN{x,p) = [ [ M{x,p)p^A'^T.^d^p, (A.l) 

J Ap JASS 

where 

M{x,p) = V / rfr5(4) [x - x,(r)] <5(^) [p - p,{t)] , (A.2) 

is a measure of the world line density. 

It is easy to construct a new function M\x^p) such that the physical constraints, the 
mass-shell condition and positive energy are factored out into the phases space measure dP = 
■^^^ 26{pq))6{p^ — mljj), in such a way that J d'^pj\f{x,p) = J dPJ\f'{x,p) holds. The construc- 
tion goes as follows. 

1 /2 

Owing to the fact that = (p^ + ^Ifj) the delta function containing p° can be factored 
out from ( |A.2| ). Then, the proper time variable r is replaced with and changing the integration 
variable from proper times Xj to the time variable x^ with the help of dri = \dx^ /dTi\~^dx^ = 
{p^/m)~^dx^, we may integrate with respect to these variables. 

N 

We obtain: J\f'{x,p) = (27r)^ ^ rf^^^^ [x - Xi(t)] 5^^^ [p - Pi(t)] , where the positions and 

i=l 

momenta of the particles are written as functions of time. Using the properties of the 6 function 
one can show that both functions satisfy the Boltzmann equation (|3.10|). 



A statistical averaging over the particle positions and momenta introduces the one particle 
distribution function: 

(AT) = 2e{po)6{p' - m%f){U') =: {27if2e{po)5{p' - m%f)n{x,p). (A.3) 
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Next we present the evaluation of three integrals appearing in Eqs. ( |3.25| ), ( |3.35| ). 

r jA 1 poo poo 

h = J -^^Mp) ^ ^ J JPo <i|p||pri(pS-i^»)(8(p„) + fi(|p„|)), (A,4) 
/, = f JpL^2nHp^ - mY-^, (A,5) 

J (277)4 /c-p 



3 



/oo roo a2 

dp, ^;|p||pp^5b^-,;2^p))(e(po)+n(bo|)).(A.6) 



(27r)4 27r2 
We use the property of the 5 function: 

5{pl-u:\^)) = -^(5{p,-uj{Y^) + 5{p, + uj{jp))), with uip) = {V^ + M'Y". (A.7) 

After the evaluation of the po integral the two terms arising from the use of ( |A.7| ) turn out 
to be identical in the case of the three integrals above, that can be easily seen using the change 
of variable p — > —p. Here we also have assumed that the 0-field consists solely of a single 
mode, characterised by the momentum 4- vector (co',k). 

The temperature-dependent part of integral (|A.4| ) is easy to evaluate in the high temperature 
limit. In dimensionless variables it can be rewritten as 

oo 

Po= 



Making use of a high temperature {l3M << 1) expansion in the above integral, namely 



F M = / dx- = ^ + 0(u^ Inu)] , A.9 

^ ^ I e"^ - 1 m2 \ 12 47r ^ ^ ' ' ^ ^ 



we obtain 

A(^) = (A.IO) 

which reproduces the thermal mass expression for the low-Zc modes correctly. 
For the integral ( |A.5| ) one obtains: 

' '^°I,nr%l/^.^ ' , , . (A.11) 



47r2 |k| Jo J_i uj{p) uj{p)ko/\k\ - \p\y + ie duj{p) 

oo 



47r2 |k| J 7_i A;o/|k|t- Vt^^y + ze 

oo 1 oo 

[ dt dys(y-^^L=]^^=-^^ [ dt^^. 
47T\k\J J \ \k\y/W^J dt 47r|k| J dt 

1 -1 1 _ 

Y^l-(fc0/|k|)2 
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where for the second hne we used the substitution t = Lj{\p\)/M and have introduced n°(t) 
l/(e^*^ — 1) Q(Mt — A). For the third hne we apphed the Landau prescription 



1 ^1 . N 

= V l-KOiz). 



(A.12) 



For the integral (|A.6|) we use the relation ddi/p^) I dp^ = ddi/p^) I dpo and perform a partial 

integration in the Pq integral: 



1 



47r2 
1 

4^ 



rf|p||pp 
d\p 



dpoSipl-uj^ip)) 



1 dn{\po\) n{\po\) Q{po) 



1 fl 



uj{p) 



n(^(p)) 



_Po d\po 

dn{uj{p)) 



pI 



duj{p) 



(A.13) 
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Equivalence of the one-loop perturbation theory and the once iterated solution of 
the generalised Boltzmann-equations 



Instead of the generalised Boltzmann-equations described in Subsection [4.1.2| we can use also 
perturbation theory to evaluate the two-point correlation functions. Here we will demonstrate 
the equivalence of the one-loop perturbation theory and the iterative solution of the generalised 
Boltzmann-equations in the case of A«. 
In the perturbation theory we write 



{M^)Mx)) = (T,0W(a;)0r(x)e-^^)^, (B.l) 



where cj)^ ' are the free fields, Sj is the part of the action which describes the interaction between 
the different field components in the presence of the background. Tc stands for the time ordering 



along a complex time path C specified in |^8[. At one loop, to linear order in $ we need only 

Si = ^ [ dxoe [ d'x [$i + 2{^k<f>b)(f>i] , (B.2) 



c 



where the integration variable xqc represents the points on the complex integration contour in 
the t plane. The field operator contractions are performed with help of the matrix propagators 



' ' ^G>,ix) iG^,{x) i I (0.(x)0,(O)) (T*0,(x)0,(O)) ' ' 



where T* denotes anti time ordering. The tree level propagators are diagonal Gab{,x) = 5abGa{x). 
Since the background depends on the real (not the contour) time we find 

(0a(x)0,(x))« = -^(0i°)(x)0f (X)5,) 

= ^ /a [^iWa, + ^aWbl + ^bWal]Sa,{x - y) , (B.4) 

where 

tSatiz) = G^(^)Gf (^) - G<iz)G<iz). (B.5) 
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The propagators can be expressed with help of the spectral function p{p) = iG^ (p) — iG^ (p) 
as 

G<{p) = n{po)p{p), G^(t,p) = e(t)p(t,p) + G<(t,p). (B.6) 
Using the {t, p) representation and finally performing time Fourier transformation leads to 

d^p dpo dp'o paipo, p)pb{p'o, P + k) 



Sab{k) 



'l + n{po) +n{pQ)). 



(B.7) 



(27r)3 2n 2n ko — po — p'q + ie 
Using free spectral functions Pa{p) = (2vr)e(po)'^(p^ — "^a) arrive at the known expression 



(see for example |pO| ) 

Sab{k) = 



d^p 1 

(27r)3 AuJatOb 

ria - rib 



l + Ua + nb 



l + Ha + nb 



ko — Ua — ujb + ie ko + uja + ujb + ie 
^ rib-Ua 



+ 



ko + Ua — uJb + ie ko — uja + ujb + ie 



(B.8) 



On the other hand introducing in ( p.7|) A^^\p) = G^{p) = (1 + n(po))p{p) = G'^{—p), next 
performing one of the po integrals, and shifting properly the p integral we obtain 



Sabik) 



d'^p 



Ai°i(p) 



+ 



^2 _|_ p2 _ 2kp — ml k"^ + p^ — 2kp — ml 



(B.9) 



The ie is assigned to ko by the Landau-prescription. Exploiting that p{p) ~ 6{p'^ — m^) one can 
write 



Sab{k) 



d'^p 



+ 



fc2 - 2kp - M2 A;2 + 2kp + 



(B.IO) 



where = m^ — m^. Finally performing a ±A;/2 shift the result is 



Sabik) 



d^p Af^{p-k/2)-A':'>{p + k/2) 



(0), 



:b.ii) 



(27r)4 2kp + M2 

Introducing for the Fourier-transform of the left hand side of Eq. ( [B.4| ) the representation 



i^a^b) 

we obtain from Eqs. ( p.4|) and 



(B.12) 



2kpA(^Xk,p) = -XaM^ik) [AW(p + A;/2) - A^^^{p-k/2)] , 



{2kp + Ml) A«(A;,p) = --^Uk) [^t\P + ^/2) - A^^P - ^/2) 



A. 



.(0), 



(0), 



(B.13) 



which exactly coincides with Eq. ( [4.25|) . 
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The dynamics of the classical 0{N) model 



An alternative approach for calculating characteristic quantities of real-time correlation func- 
tions is based on real-time dimensional reduction [|127| , p.28|| and solving the resulting classical 
effective theory relevant for the modes with high occupation numbers. On-shell as well as 
the off-shell damping rates were successfully reproduced p4|, ^ in the symmetric phase of 0^ 
theory using this approach. In this appendix we will discuss the application of the classical 
approach to the 0{N) model and compare its results with the exact one- loop dynamics. The 
Lagrangian of the classical 0{N) theory has the following form: 

1 



24 I a) 



The classical equation of motion corresponding to this Lagrangian is: 



0. 



(C.l^ 



(C.2) 



The external currents ja{x) were introduced in order to prepare the derivation of the classical 
response theory. They should not be confused with the induced currents Ja appearing in the 
effective quantum equations of motion. The classical mass mc/(A) is different from the mass 
parameter m of the quantum theory (see Eq. ( [4. 2D ). The same is true for the coupling Ac;. The 
classical mass parameter depends on the ultraviolet cutoff A. The results of the classical and 
quantum calculation could be matched by a suitable choice of mci{A) and Ac/. In particular 
we shall see below (eq.( |C.li )), that the ultraviolet divergences of the classical theory can be 

Therefore we 



eliminated if the divergent part of mc/(A) is suitably chosen |24, 47, |127| , [l2l 
will separate out the divergent part from the classical mass and write m^;(A) = + 5m^(A) 
Q Note that mx is not only the thermal mass, see Eq. ( |C.23|) . 

In the broken symmetry phase one separates out the condensate $, 



(C.3) 



^The divergent part of the classical mass parameter will be treated as interaction, similarly to quantum field 
theory. 



Ill 



and the equations of motion read 



2 / D 2 3 



In addition the following initial conditions are imposed: 

= 0,X) = F„(X), 9i$a(t,x)|4=o = Pa(x). 



(C.4) 
(C.5) 

(C.6) 



The expectation value of some quantity O (e.g. some correlation function ) is obtained by 
averaging over the initial conditions with the Boltzmann factor determined by the classical 
Hamiltonian Hci{Pa, Fa,^) corresponding to ( |C.1|) 

1 



<0> 
Z 

The explicit form of Hci{Pa, Fa, ^) is 



DF,DP,Oexpi~PH,i{Pa,Fa,^)) 



DFaDPa eM-l3Hci{Pa, Fa, ^)), 



(C.7) 
(C.8) 



H,i{Pa,Fa,^) 



d X 



^miF^r + ^{F^af + «^ + Y^')F, + \ml^' + + 6m^F^] , (C.9) 

where we have separated out the classical condensate $ since no averaging over classical con- 
densate is understood. Since the classical condensate is separated both from the dynamical 
fields and from the initial conditions the following equation holds: 



< = 0,x) >=< Fi(x) >= 0. 



(C.IO) 



Using the explicit form of Hd (Eq. (|C.9|) ) this results at one-loop level in the following equation 
for the classical condensate $: 



T 



+ (iV-l) 



T 



p2 + ml + 



0. (C.ll) 



Choosing 6m^{h) to cancel the linearly divergent piece on the left hand side this relation reduces 
in the high temperature limit {mx, ^ T) to 



and 



(C.12) 
(C.13) 
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Our procedure of solving the classical theory perturbatively closely follows Ref . |^ . Equations 
), ( |C.5| ) are rewritten in form of the following integral equations: 



where 



J V 2 6 



At) 



(2vr) 



(C.14) 
(C.15) 



is the classical retarded Green function with Ug = \J + ml, where 



ml + 



is the Higgs field mass and mf = m'^ + -^^'^ = (we have used Eq. (|C.12| )). Furthermore, 
$° are the solutions of the free equations of motion and j in the argument of refers to the 
functional dependence on j 
functions 

x') = (C.16) 



(ji,jj). Following Ref. one introduces linear response 

5*a(a;,j) 



Using the integral equation (|C.14| ) one can derive a coupled set of integral equations also for 
the linear response functions H'^ by functional differentiation of eq. ( |C.14|) (see Ref. for 
details). These integral equations can be solved iteratively in the weak coupling limit. When 
solving the equations, it is important to exploit the fact that only diagonal components of 
have terms 0{X^). The classical retarded response function is the ensemble average of with 
respect the initial conditions : 



X 



--< H%\x ~x')> . 



(C.17) 



By eq. (|C.14|) one easily finds that it satisfies a Dyson-Schwinger equation of the general form 



G^Kx - x') 



DUx - x')Sat + J dSd^DUx - y)SadIli{y - y')G^tiy' - x'), (C.18) 



where Habiv ~ v') classical self-energy. 

In the perturbative expansion the average is done with the free Hamiltonian and therefore 
all thermal n-point functions are expressed as products of the two-point function of the free 
fields (solutions of the free equation of motion). The two point function of the free fields reads 



PI, M as 



A:'°(x - x') =< <l>°(x)„$0(x')a >°= T 



d^q 
(2^ 



^iq(x-x') 



(C.19) 



This classical two point function is analogous to the free two point function of the quantum the- 
ory Aia''(x, x') Using eq. ( |(J.14| ) one gets the following self-energies for the Higgs and Goldstone 
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fields at leading order in the coupling constant Xd- 



a 

— 2 

(^) I dtd'ye^-'-^^y {At\y)D'niy) + ^f'\y)DUy)) , (C.20) 



Using the explicit form of A^''°(0,0) and Eq. ( |C.13 ) one can easily verify that all divergencies 
present in the above expression cancel. Furthermore, in the high temperature limit and at 
leading order in the coupling constant only the last terms contribute in the expression of 
nf]^(c<j,k) and Ilfj^{uj,k.). One has to evaluate integrals of the following intrinsic form 

jdtj d'ye-^~^-yDUy)^t%) = IJ 2^ (^"(^o) + n^\p'o)) , (C.21) 

where uja = a/p^ + ml, ooh = a/(p + k)^ + and /p = / (^~^- The above expression coincides 
with the result of the quantum calculation (|B.7|) , except the fact the there is no T = con- 
tribution and the Bose-Einstein factors are replaced by the classical distribution: T/pq. Then 
it is easy to write down the explicit expression for the classical self-energies using the formal 
analogy with the result of the one-loop quantum calculations. For example for the classical 
on-shell damping rate of the Goldstone modes one easily gets the following expression: 

— 11 I 

, Imng(cu= |k|,k) A2;$2 /-Tw+W ^ 
rf (k) = " oiui = " o. ,2 dpn'^'ip) 




(C.22) 

Now let us discuss the correspondence between classical and quantum calculations. It was 
shown in Refs. [24, 47, 127 , 128|| that the result of classical calculations can reproduce the high 
temperature limit of the corresponding quantum results if the parameters m^;(A) and Xd are 
fixed to the values determined by static dimensional reduction. For our case this implies: 

^^z(A) = m^ + ^(iV + 2) (--—), A,, = A, (C.23) 

where m and A are the mass and the coupling constant of the corresponding quantum theory. 
From the above equations it is easy to see that the divergent part in m^i coincides with 
determined from the ultraviolet finiteness of the classical result (cf. Eqs. ( |Cll| ),( |Cl3| ) and 
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( p.20| )) and mi is the high temperature hmit of Mi (see Eqs. ( |C12| ) and ( [4.12| )). For small 
values of |k| (|k| ^ mi) the logarithm in this expression can be expanded and one obtains 
the result of Eq. ( |4.39| ), which fails to reproduce the result of the quantum calculation. The 
same is true for the whole imaginary part of the Goldstone self-energy. The high temperature 
limit of the imaginary part of the Higgs self-energy (see Eq. ( [4.29| )), on the other hand, is well 
reproduced by the classical theory. 

This result can be easily understood by looking at the effective quantum equation of motion 
( [4.10| ). In the equation for the Goldstone fields the induced current Jj has a non-local contri- 
bution from loop momenta p ~ Mf/|k| >> T (cf. (|4.32|) ). However, no such term is present 
in the corresponding classical equation of motion on one hand and the classical theory cannot 
describe fluctuation with wave length much smaller than on the other hand. The induced 
current Ji in the effective equation of motion for the Higgs fields receives non-local contribution 
only from the loop momenta around p ~ Mi which can be described in the framework of the 
classical theory. 
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Appendix D 



Evaluation of some integrals appearing in Subsection 4.1.3 



In this Appendix we illustrate the steps of evaluation of the integrals in Eq. (|4.28|) . 



Imaginary parts. As example of the evaluation of a relevant integral we discuss in detail 

d^p A(o) (p - A;/2) - AW {p + k /2) 



Imi?i(A;,M) = Im j 



(27r)4 pk ' ^^'^^ 

where the Wigner-functions A*^"^ can be either of type 11 or ii. In order to implement the 
Landau-prescription we transform away any k dependence from A'^''^ by shifting the p integral 
by ±.k/2. Then we use 

(D.2) 



lim Im- 



— = — 7re(a)5(x), 
a^o x + ia V / w 

and finally shift the integrals back. We write the 4D integration measure as 

oo oo 1 

1 



dp 
(2^ 



87r3 



dpp^ / dpo I dx, 



:d.3) 







-1 



where x = pk stands for the cosine of the angle between the spatial momenta. The x integration 
is trivial because it appears in the Dirac-delta arising from the application of the principal value 
theorem 



/ dx6{poko - p\k\x) = -l—Q{p\k\ - \poko\). 



(D.4) 



Using the explicit form of the Wigner-functions A*^°) (see Eq. ( |4.24| )) and the identity 0(u;) + 
n(|u;|) = e(to')(l -|- n^u)) we find 

oo a 

Iff k k 

-lmRi{k,M) = -^-^^ j dppj dpot{pQ--^)e{po + -^) X 

-a 



5{pl-S' 



(D.5) 
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where a = p|k|/|A;o| and S"^ = + — The po integration over the Dirac-delta gives a 

constraint for the p integration of the form 

S<a =^ p2-2< — -Ml (D.6) 

Kg 4 

This can be fulfilled only for k'^ > 4M^ (above the two-particle threshold), or for A;^ < 
(Landau damping). After elementary algebra one can establish the value of the sign functions 
and one arrives at the formula appearing in Eq. (f4.32| ). 



A similar analysis can be performed for Imi?j, however, in this case it proved to be more 
convenient to start from the equivalent form 

M) = Im / Ai?(p-t)-Air(ri 

' ^ J (27r)4 2pk-k^ + M^ ^ ' 

After implementing carefully the Landau prescription and performing the x integration as 
described before, we arrive at 

-Imi?i(/c,M) = j dppj dpo e{po) e{po-ko) 5{pl-p'^) [n{po-ko) - n{po)] , (D.8) 

6_ 

where 

e-M^ iki , , 

The Po integration over the mass-shell delta-function again restricts the domain of integration 
in the p integral: 6_ < ±p < 6+, which, however, does not restrict the possible values for fc^. 
After the solution of these linear inequalities and the analysis of the sign functions we arrive 
at the result appearing in Eq. (|4.32| ). 

Real parts. The relevant integrals in Rei?i are 

J (27r)4 pk±k^/2' ^ ^ 

where the Wigner-f unctions A'^''^ can be either of type 11 or ii. With their help we find 
Rei?i = — . The real part comes from the principal value integration. We decompose 



the integration measure as in the calculation of the imaginary parts. When we use Eq. ( 4.24 ) 
for A^^\ the value of po is fixed by the mass-shell delta function. The x integration can be 
performed as 

1 

-1 



118 



where arth(x) = l/21n|(l + x)/(l — a;)|. Using the properties of the absolute value we find 



/+(A;o,k) = -/-(-A;o,k). 



(D.12) 



Then we directly arrive at Eq. ( [4.30|) . 
We write Rei?, in the form 



ReRi = Re 



(27r)4 2pk + k^ + M^ 



Re 



(27r)4 2pfc-P + M2' 



(D.13) 



Then the previous scheme of calculation goes through directly. 

The only problem still to be discussed is the zero temperature contribution, which diverges 
logarithmically. The regularization and renormalization prescriptions are better formalized in 
the language of the propagators and couplings. Using the results of Appendix ^ to go over to 
the perturbation theory, we have to evaluate 

d^p 



R 



J [G''ip)G''ip-k)-G<ip)G<ip-k 



(D.14) 



The first term at T = is the usual time ordered product; with finite four- dimensional cutoff 
one finds 



d'^p 



G^{p)G''{p-k) 



-1 
16^ 



1 + J dx \n 





\k'^x{l — x) — mix — m^il — x) 

A2 



(D.15) 



The divergence is canceled by the coupling constant counterterm. In modified minimal sub- 
traction (MS) scheme it reads 



, , Ik'^xil — x) — m^x — 171^(1 — x) 
dxln- ^ i — 



(D.16) 



IGtt^ Jo a* 
In the second term of Eq. (p.l4|) at T = we can use G^{p) = 0(— po)(2vr)5(p^ — 



Because of the delta functions this piece yields finite result (as is expected by the arguments of 
the renormalizability) 



d^v 1 
^ G<{p)G<{p-k) 



(2vr) 



in 



1 — 0(A;^ — 4m^), if mi = 1712 = m 
(1 — ^) 6(/c^ — m^), if mi = 0, m2 = m. 



(D.17) 
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Summary of main results 



In my thesis I discuss the real time non-equihbrium dynamics of systems in the broken symmetry 
phase. Beside reviewing the fields in which this kind of investigation are relevant, I presented 
the equation derived to be used in the broken phase of scalar models and in the Abelian Higgs 
model. 

The main results of the analytical and numerical investigations are as follows: 

1. A classical kinetic theory was constructed for particles whose mass square depends quadrati- 
cally on the local amplitude of a real, low frequency scalar wave. It was shown that the damping 
rate of the waves due to the scattering of the particles agrees with the result of the quantum 
perturbation theory in the broken phase of the one-component <I>-theory. 

2. The effective equation of motion of the order parameter of the real theory extracted from 
numerical simulations follows the Maxwell construction when a non-equilibrium first order 
phase transition occurs in the 2 + 1 classical theory. The statistical features of the decay of the 
false vacuum agrees with the results obtained by expanding around the critical bubble. 

3. The damping rate of the quantum Goldstonc modes was determined to leading order in the 
0{N) model with perturbation theory for low wave number k and in presence of small explicit 
symmetry breaking (h). It was shown that the rate of relaxation vanishes non-analytically 
when |k| — ^ and /i ^ 0, as a direct consequence of Goldstone's theorem. The discrepancy 
between the classical an quantum results was discussed. 

4. The analytically determined large time asymptotic power law of the relaxation of an arbi- 
trary field configuration has been checked in a numerical simulation of the 2-1-1 dimensional 
0(2) symmetric classical field theory. In the linear regime of the relaxation the numerically 
determined damping rate of the Goldstone mode was found to be in agreement with the one 
calculated in the classical finear response theory. The real time evolution of the symmetry 
breaking ground state into the symmetrical one was checked. The suppression of the order ex- 
pected on the basis of the Mermin- Wagner theorem was checked to proceed through a diffusive 
mechanism. 

5. The current density of the Maxwell equation has been determined to O(e^) accuracy in scalar 
electrodynamics. This is induced by the presence of a mean gauge field and also influenced by 
a symmetry breaking homogeneous scalar background The result of quantum perturbation 
theory displays a contribution proportional to in addition to the well-known "hard thermal 
loop" (HTL) expression. The gauge-fixing parameter independent degrees of freedom were 
constructed. 
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A fobb eredmenyek osszefoglalasa 



A doktori ertekezesben a sertett szimmetriaj li fazis valosideju nem-egyensiilyi dinamikajat 
vizsgalom. A valos idejii nem-egyensiilyi dinamika temakorehez kapcsolodo irodalom ismer- 
tetese mellett bemutatom az altalam levezetctt egyenleteket, melyek a skalartereket illetve az 
abeli mertektcrckct tartalmazo elmeletek scrtctt szimmetriaj u fazisaban hasznalhatok. Anali- 
tikus vizsgalatokat skalar terelmelctckben cs az abeli Higgs modellben vegeztem. Numerikus 
szimulaciokat a skalar terelmelet kcrctei kozott folytattam. 

A fobb eredmenyeket a kovetkezo pontokban foglalom ossze: 

1. Megkonstrualtam egy klasszikus kinetikus elmeletet az onkolcsonhato egyensulyi skalartcr 
nagyfrekvencias modusait reprezentalo gazra. E gaz relativisztikus skalar rcszccskcinek tomege 
kvadratikusan fiigg a hosszuhullamu hatter amplitudojatol. Megmutattam, hogy a hosszuhul- 
lamii modusok csillapftasi rataja, amely a reszecskeken valo szorodas kovetkezmenye, a sertett 
fazisban megegyezik a kvantum egy-hurok szamolas eredmenyevel. 

2. 2+1 dimenzios skalar elmelet nem-egyensiilyi elsorendii fazisatalakulasa numerikus szimu- 
laciojanak adataival egy effektfv mozgasegyenletet konstrualtam a rendparamcterre. Megmu- 
tattam, hogy a rendszer dinamikajat leiro effektfv mozgasegyenlet tiikrozi a Maxwcll-konstruk- 
ciot: a rendparameter a stabil ertekek kozott eromentes mozgast vegez. A mctastabil — > sta- 
bil atmcnct idostatisztikajabol adodo atmeneti rata osszhangban van a nukleacios elmeletben 
szamolt rataval. 

3. Az 0(N) szimmetriajii modell scrtctt fazisaban kiszamoltam a pcrturbacioszamitas elso 
rendjeben a Goldstone-modus csillapftasi ratajat a |k| hullamszam, illctvc az explicit szimme- 
triasertest jellemzo h kiilso ter fiiggvenyeben. Megmutattam, hogy a Goldstone-tetel ertelme- 
ben, a |k| ^ modusok csillapftasi rataja h = 0-ra ncm-analitikus modon tfinik el. Kimutattam 
az elterest a kvantumos cs a klasszikus clmelctbol szamolt csillapftasi rata kozott. 

4. 2 + 1 dimenzios numerikus szimulacioval leellenoriztem az 0(2) szimmetriajii skalar ter ana- 
litikusan mcgliatarozott, aszimptotikus idokrc vonatkozo rclaxacios torvenycnck teljesiileset. 
A relaxacio linearis szakaszaban a Goldstonc jcllcgi'i szal:)adsagi fok mcrt csillapftasi ratajat 
osszhangl:)an talaltam a klasszikus linearis valaszclmclctbol szamolt rataval. Kimutattam a 
szimmctriascrto allapotnak a valosidcjii atalakulasat a szimmctrikus alapallapotban. A Mer- 
min-Wagner-tetel dinamikai megvalosulasara diffiizios jcllcmzest adtam. 

5. Skalar elektrodinamikaban O(e^) rendfi pontossaggal meghataroztam a hossziihullamii mer- 
tekterek Maxwell egyenleteben szereplo indukalt aramot. Az aramot a mertek atlagter in- 
dukalja cs tartalmazza a szimmetriasertcst tiikrozo homogen skalar hatter hatasat. Ennck $^ 
rendfi figyelembevetele tiilmutat az irodalomban szereplo, un. "hard thermal loop" (HTL) 
egyenleteken, hatasa fontos lehet a fazisatalakulas lefrasakor. Megmutattam, hogy a fizikai 
szabadsagi fokok fiiggetlenek a mertekrogzfto parametertol. 
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